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CHAPTER  I 


INTRODUCTION 


Current  evidence  shows  that  deoxyribonucleic  acid(DNA) 
is  the  carrier  of  genetic  information  within  the  nucleus  of 
the  cell.  In  protein  synthesis,  the  DNA  is  transcribed  into 
a  complementary  chain  of  messenger  ribonucleic  acid(mRNA). 
The  mRNA  chain  is  "read"  in  the  cytoplasm.  The  reading 
process  involves  the  incorporation  of  either  an  amino  acid 
or  the  termination  of  the  amino  acid  chain,  depending  on  the 
specific  three  nucleotide  sequence  of  the  mRNA  chain. 

The  genetic  code  which  translates  the  nucleotide 
triplet  into  a  specific  amino  acid  (or  a  chain  termination 
signal)  has  been  found  to  be  universal  across  all  species  of 
life  so  far  studied  (Epstein,  1 966) (Mel char ,  1970).  The 
implication  is  that  all  forms  of  life  have  evolved  after  the 
genetic  code  became  fixed  in  its  present  form  some  time  in 
the  distant  past  (Conrad,  1970) (Fox,  1 97 4) ( Har tman , 
1975)  (Jukes,  1973)(Orgel,  1972).  The  theoretical  problem  of 
the  origin  of  life  as  known  on  this  planet  can  then  be 
divided  into  two  component  parts.  First,  the  problem  is  to 
understand  the  chemistry  and  exobiology  responsible  for  the 
development  of  protocells  and  the  fixation  of  the  genetic 
code  represented  by  the  universality  of  the  genetic  code 


(Chernavskn,  1  97  5)  (Keosian,  1974)  (Lacey,  1  975)(Mazin, 
1975)  (Mikelsaar,  1975).  The  second  problem  is  to  understand 
the  mechanisms  at  work  after  the  fixation  of  the  genetic 
code.  Within  this  second  area  fall  the  mechanisms 
responsible  for  the  range  and  diversity  of  life  as  observed 
today  (Conrad,  1 97 0) (Geraci tano ,  1971).  The  range  and 
diversity  is  represented  by  the  great  number  and  complexity 
of  species  as  well  as  the  vast  amount  of  difference  observed 
within  each  species  (Crow,  1 97 0) (Geracitano,  1971) 

The  evolutionary  theories  of  Darwin  and  the  genetic 
systems  of  Mendel  have  been  instrumental  in  providing  a 
framework  for  categorizing  and  understanding  many  of  the 
processes  that  have  occured  since  the  fixation  of  the  code. 
The  theories  of  Darwin  and  Mendel  formed  a  school  of  thought 
which  examined  phenotype  to  form  generalizations  about 
genotype.  This  methodology  has  been  continued  by  the  modern 
population  biologists  (Jungck,  1 97 1)  (Lewontin,  1966)  (Ohto, 
1971).  The  main  focus  is  on  changes  in  population 
statistics  as  observed  and  applied  using  neo-Darwinian 
theories  which  emphasize  selective  advantage,  or 
disadvantage  (Epstein,  1 96 7 ) (K imur a ,  1 96 9) ( Papen t in , 
1973)  (Richmond,  1  970) . 

After  the  identification  and  characterization  of  the 


genetic  code,  many  researchers  began  to  look  at  the 
phenominon  of  evolution  from  the  genotypic  level  (as 
represented  by  amino  acid  sequence  within  protein)  towards 


the  phenotypic  level.  The  redundant  nature  of  the  genetic 
code  along  with  known  characteristics  of  DNA  enabled  certain 
statistical  projections  about  amino  acid  sequence  and 
nucleotide  sequence  to  be  made.  Random  point  mutations  of 
the  DNA  chain  could  yield  statistical  evidence  for  changes 
in  the  amino  acid  structure  of  protein  over  a  time  period 
long  enough  to  reflect  evolutionary  processes.  Such 
projections  lead  to  surprisingly  accurate  correlation 
between  predicted  time  of  seperation  of  species  and 
archaeological  data  (Acher,  1974) (Clarke,  1970).  Amino  acid 
differences  within  species  (known  as  polymorphic  proteins) 
could  be  partially  explained  by  random  point  mutations  of 
neutral  selection  value  (Kimura,  1 9 7 1 )  (  S u e o k a , 
1961)  (Zuckerkandl,  1972). 

The  overall  objective  of  this  thesis  will  be  to  provide 
a  framework  to  examine  the  random  mutation  hypothesis.  From 
this  framework  a  better  understanding  of  the  possibilities 
and  limitations  imposed  by  random  mutation  can  be  evaluated. 
Implications  about  the  evolutionary  potential  of  the 
information  carrying  system  and  the  present  equilibrium  of 
that  system  may  be  possible. 

This  thesis  is  written  in  five  chapters.  This,  the 
first  chapter,  has  presented  a  brief  introduction  to  several 
of  the  ideas  of  current  interest  in  evolutionary  biology. 
This  section  identifies  the  organization  of  the  remainder  of 
the  paper.  The  second  chapter  reports  the  relevant 
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information  gathered  during  the  course  of  the  research  into 
the  evolutionary  potential  of  bio-organic  information 
transmitting  compounds,  the  characteristics  of  mutations 
occurring  within  the  system,  and  information  on  hemoglobin's 
evolutionary  and  mutational  behavior.  The  third  chapter 
addresses  the  methodology  to  be  employed  in  this  thesis. 
This  includes  the  development  of  a  modeling  approach  to  the 
point  mutation  of  the  DNA  sequence  responsible  for  producing 
viable  beta  hemoglobin  chains  as  well  as  the  techniques  used 
to  evaluate  the  model's  performance  and  the  results  on 
hypothesis  testing.  The  fourth  chapter  lists  the  results  of 
analysis  I  performed  on  the  simulation  results  obtained. 
The  fifth  chapter  reflects  the  conclusions  that  I  have  been 
able  to  gather  from  the  modeling  and  analysis  results 
obtained  to  date. 


CHAPTER  II 


BACKGROUND 

The  first  section  will  be  a  quick  tutorial  of  different 
portions  of  the  information  transmission  system  of  cells. 
This  will  consist  of  brief  outlines  of  the  compounds  and 
simplified  descriptions  of  the  functioning  of  the  system. 
This  section  will  also  include  a  review  of  genetic 
mutations.  The  second  section  will  contain  a  bibliographic 
review  of  the  concept  of  evolutionary  forces  at  work  in  the 
mutation  of  DNA.  Different  theories  as  to  the  operate 
mechanism  will  be  discussed  as  well  as  some  of  the 
implications  of  these  theories. 

The  chromosomes  of  organisms  are  composed  of 
nucleoproteins.  Proteins  such  as  histones  and  nucleic 
acids  comprise  the  nucleoproteins.  The  nucleic  acids  have 
been  associated  with  the  transmission  of  genetic  information 
while  the  histones  play  a  role  in  the  stabilization  of  the 
compounds  and  mechanisms  invol ved(Stansf ield,  1969)  (Stryer, 
1975)  . 

Deoxyribonucleic  acid  is  the  repository  of  genetic 
information  in  all  eukaryotes.  DNA  is  stabilized  into 
chains  thousands  of  nucleotides  long  that  form  a  double 
stranded  helix.  The  backbone  of  each  chain  consists  of  five 
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carbon  sugar  compounds  oxygenated  at  the  2'  position  linked 
by  phospate  to  form  5'=>3'  phosphodiester  linkages.  The 
helix  is  held  together  by  the  interaction  of  the  bases 
attached  to  the  1'  carbon  of  the  sugar.  The  paired  organic 
bases  found  in  DNA  are  the  purines  adenine (A)  and  guanine (G) 
and  the  pyrimidines  thymine(T)  and  cytosine(C).  In  normal 
DNA  with  unmodified  bases,  cytosine  always  pairs  with  three 
hydrogen  bonds  with  guanine  while  adenine  always  pairs  with 
two  hydrogen  bonds  with  thymine.  The  nucleotide  unit 
consists  of  the  sugar,  the  phosphate,  and  the 
base (Stansf  ield,  1969)  (Stryer,  1975). 

Ribonucleic  acid(RNA)  is  another  class  of  nucleic  acid. 
RNA  differs  from  DNA  in  that  the  base  sugar,  again  a 
pentose,  does  not  have  an  oxygen  at  the  2'  position.  RNA  is 
thought  of  as  primarily  a  single  stranded  nucleotide, 
however,  the  ability  of  RNA  to  stabilize  itself  by  forming 
double  and  even  four  stranded  chains  is  thought  to  play  a 
key  role  in  several  of  the  mechanisms  and  the  results  of 
several  thermodynamic  considerations.  RNA  constains  the 
pyrimidine  uracil(U)  instead  of  thymine  within  its 
structure.  RNA  is  typically  much  shorter  than  DNA  and  is 
found  both  within  and  outside  of  the  cell 
nucleus(Stansf ield,  1969) , (Stryer,  1975). 

RNA  plays  a  predominate  role  in  protein  synthesis. 
There  are  three  common  types  of  RNA  designated  as  mRNA, 
tRNA,  and  rRNA.  Messenger  RNA,  or  mRNA,  is  transcribed  from 


DNA  within  the  nucleus  of  the  cell  and  is  brought  to  the 
cytoplasm.  Transfer  RNA,  or  tRNA,  is  found  in  the  cytoplasm 
and  is  active  in  selectively  positioning  the  appropriate 
amino  acid  into  the  growing  protein  chain  based  on  the  codon 
of  the  mRNA  carrying  the  template  for  protein  synthesis. 
Ribosomal  RNA,  or  rRNA,  is  actually  a  class  of  different 
molecular  weight  components  with  protein  components  which 
serve  as  the  binding  site  for  mRNA  and  tRNA  in  the 
translation  process (Stansfield,  1969)  (Stryer,  1975). 

The  individual  hydrogen  bonds  linking  the  bases  of  the 
DNA  chains  are  not  strong  compared  to  typical  covalent 
chemical  bonds  such  as  those  that  connect  the  sugar  to  the 
phosphate  or  to  the  base,  however,  the  large  number  of  bonds 
form  a  very  powerful  and  a  very  stable  chain.  DNA  is 
replicated  with  the  assistance  of  enzymes  known  as  DNA- 
polymerases  which  assist  in  the  separation  of  the  DNA  helix. 
The  two  single  DNA  strands  each  form  a  new  double  stranded 
DNA  helix  by  the  addition  of  nucleotides  in  a  new  second  DNA 
strand  which  forms  in  the  5,=>3'  direction.  The  replication 
of  DNA  occurs  extremely  fast.  As  a  practical  consideration, 
single  stranded  DNA  is  very  rarely  found(Stansf ield, 
1969)  (Stryer,  1975). 

The  nucleotide  sequence  of  DNA  is  the  blueprint  for 
protein  synthesis.  Proteins  form  the  major  structural 
components  of  the  cell  as  well  as  the  enzymes  necessary  to 
catalyze  the  biological  reactions  which  are  essential  for 
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life  as  we  know  it.  Proteins  are  composed  of  a  large  number 
of  amino  acid  residues  which  forms  the  primary  structure. 
Proteins  form  stable  secondary  and  tertiary  structures  based 
on  the  thermodynamic  equil librium  of  the  chain  in  the  ionic 
environment  in  which  it  exists.  While  the  primary  structure 
of  the  protein  may  be  responsible  for  the  production  of  an 
active  site  in  an  enzyme,  it  is  the  secondary  structure 
which  is  reflected  within  the  environment  which  ultimately 
determines  the  catalytic  ability  of  the  protein.  If  the 
ionic  or  temperature  conditions  cause  the  secondary 
structure  to  change,  the  catalytic  ability  of  the  enzyme  is 
usually  lost.  Proteins  are  comprised  of  20  different  amino 
acids  as  reflected  in  table  2.1 (Stansf ield,  1969) (Stryer, 
1975)  . 


Table  2.1 

Amino  Acids  By  Decreasing  Frequency  of  Occurance 


1 

G 

Gly 

Glycine 

11 

I 

He 

Ileucine 

2 

A 

Ala 

Alanine 

12 

R 

Arg 

Arginine 

3 

L 

Leu 

Leucine 

13 

N 

Asn 

Asparagine 

4 

S 

Ser 

Serine 

14 

Q 

Gin 

Glutamine 

5 

V 

Val 

Valine 

15 

F 

Phe 

Phenylalanine 

6 

K 

Lys 

Lysine 

16 

Y 

Tyr 

Tyrosine 

7 

T 

Thr 

Threonine 

17 

C 

Cys 

Cysteine 

8 

E 

Glu 

Glutamic  Acid 

18 

H 

His 

Histidine 

9 

P 

Pro 

Proline 

19 

M 

Met 

Methionine 

10 

D 

Asp 

Aspartic  Acid 

20 

W 

Trp 

Tryptophan 

The  mRNA  is  translated  by  tRNA  while  being  bound  by 
rRNA  to  form  proteins.  Each  amino  acid  has  an  amino  group 
at  one  end  and  a  carboxyl  group  at  the  other.  Proteins  are 


synthesized  by  the  addition  of  new  amino  acids  at  the 
carboxyl  terminal  of  the  growing  protein  chain  based  on  the 


■reading"  of  the  mRNA  in  the  5  *  =  > 3  *  direction.  The 

■reading"  of  the  RNA  message  is  conducted  in  accordance  with 

the  genetic  code.  The  genetic  code  is  a  catalog  of  sixty- 

four  three  base  codons  which  form  a  degenerate  reference  set 

for  the  twenty  available  amino  acids  from  the  four 

Table  2.2 
Genetic  Code 

5J.CM1L  51  RNA  11  Amino  Acid  51  Dm  31  51  RNA  3'  Amino  Acid 


AAA 

OUU 

Phe 

GAA 

OOC 

Phe 

TAA 

UUA 

Leu 

CAA 

DOG 

Leu 

AGA 

tJCU 

Ser 

GGA 

UCC 

Ser 

TGA 

UCA 

Ser 

CGA 

OCG 

Ser 

ATA 

UAO 

Tyr 

GTA 

OAC 

Tyr 

TTA 

OAA 

Stop 

CTA 

OAG 

Stop 

ACA 

DGO 

Cys 

GCA 

OGC 

Cys 

TCA 

DGA 

Stop 

CCA 

OGG 

Trp 

AAG 

CUD 

Leu 

GAG 

cue 

Leu 

TAG 

CDA 

Leu 

CAG 

CUG 

Leu 

AGG 

ecu 

Pro 

GGG 

CCC 

Pro 

TGG 

CCA 

Pro 

CGG 

CCG 

Pro 

ATG 

CAU 

His 

GTG 

CAC 

His 

TTG 

CAA 

Gin 

CTG 

CAG 

Gin 

ACG 

CGU 

Arg 

GCG 

CGC 

Arg 

TCG 

CCG 

Arg 

CCG 

CGG 

Arg 

AAT 

AOD 

He 

GAT 

AOC 

lie 

TAT 

AOA 

lie 

CAT 

AUG 

Met 

AGT 

ACU 

Thr 

GGT 

ACC 

Thr 

TGT 

ACA 

Thr 

CGT 

ACG 

Thr 

ATT 

AAD 

Asn 

GTT 

AAC 

Asn 

TTT 

AAA 

Lys 

CTT 

AAG 

Lys 

ACT 

AGU 

Ser 

GCT 

AGC 

Ser 

TCT 

AGA 

Arg 

CCT 

AGG 

Arg 

AAC 

GUU 

Val 

GAC 

GOC 

Val 

TAC 

GDA 

Val 

CAC 

GOG 

Val 

AGC 

GCU 

Ala 

GGC 

GCC 

Ala 

TGC 

GCA 

Ala 

CGC 

GCG 

Ala 

ATC 

GAO 

Asp 

GTC 

GAC 

Asp 

TTC 

GAA 

Glu 

CTC 

GAG 

Gly 

ACC 

GGO 

Gly 

GCC 

GGC 

Gly 

TCC 

GGA 

Gly 

CCC 

GGG 

Gly 
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nucleotide  bases  available.  The  code  can  be  seen  in  Table 

2.2. 

The  attachment  of  an  amino  acid  to  its  tRNA  is  mediated 
by  specific  enzymes  in  a  process  called  activation.  The 
tRNA  sequence  contains  a  three  nucleotide  anticodon  which  is 
complementary  to  the  mRNA  codon.  The  tRNA  is  positioned  at 
the  translation  site  where  the  rRNA  is  "reading"  mRNA  in  the 
5,ss>3'  direction.  The  amino  group  of  the  tRNA  bound  amino 
acid  is  brought  into  a  proximal  position  to  the  carboxyl 
group  of  the  growing  peptide  chain.  Enzymes  are  active  in 
the  formation  of  the  peptide  bond  formed  at  this  site  which 
joins  the  the  amino  acid  to  the  protein.  The  tRNA  is 
released  from  both  its  amino  acid  and  the  mRNA  chain.  When 
the  ribosome  meets  the  end  of  the  mRNA  or  encounters  an  in- 
phase  termination  codon,  the  protein's  primary  structure, 
consisting  of  the  sequence  of  amino  acids  is  complete.  The 
secondary  structure,  or  the  three  dimensional  structure  of 
the  protein,  is  probably  also  complete  at  this  point  or 
shortly  after  release  from  the  ribosome.  The  secondary 
structure  is  created  by  the  formation  of  three  dimensional 
structures  such  as  the  alpha-helix  or  beta-pleat  sheet. 
Quarternary  structures  that  stabilize  the  hydrophobic  and 
hydrophillic  regions  of  the  protein  as  well  as  the  formation 
of  internal  covalent  bonds  (such  as  between  cysteine 
residues).  The  tertiary  structure  of  multiple  proteins  may 
form  a  quaternary  structure.  For  instance,  two  alpha  and 


two  beta  hemoglobin  chains  combine  to  form  the  functional 


complex.  Any  alteration  of  the  final  configuration  may 
inactivate  the  protein(Stansf ield,  1969)  (Stryer,  1975). 

Each  tRNA  consists  of  75  to  80  nucleotides  in  a 
specific  order.  The  31  end  of  all  tRNA's  contain  the 
sequence  CCA.  The  5*  end  terminates  in  a  guanine  residue. 
Pseudouradilic  acid,  inosinic  acid,  and  a  few  other  modified 
bases  are  found  in  tRNAs.  The  tRNAs  have  a  characteristic 
three  dimensional  four  arm  stucture  with  helical  sections 
and  single  stranded  loops.  One  of  the  unpaired  loops 
contains  the  anticodon(Stansf  ield,  1969)  (Stryer,  1975). 
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V  <-  t-RNA 
PROTEIN 

Figure  2.1 

Simplified  Schematic  Of  Protein  Synthesis 


More  than  30  different  proteins  have  been  associated 
with  rRNA  in  ribosomes.  Three  sites  conceptually  exist  on 
functional  ribosomes.  A  decoding  site  binds  the  activated 
tRNA  to  the  mRNA  codon.  A  condensing  site  joins  the  amino 
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acid  to  the  growing  polypeptide  chain.  An  exit  site  at  which 
|  the  tRNA  becomes  separated  from  the  protein,  mRNA,  and 

ribosome.  The  tRNA  is  processed  in  turn  thru  all  three 
'  sites  of  the  ribosome(Stansf ield,  1969)  (Stryer,  1975). 

j  Most  genes  are  relatively  stable  and  mutation  is  a  very 

rare  event.  The  average  mutation  rate  has  been  estimated  at 
one  accepted  point  mutation  occurring  in  million  genes. 
|  Other  genes  have  mutation  rates  estimated  at  one  mutation 

per  occurring  in-  ten  genes.  Mutation  at  the  genotypic  level 
can  only  be  examined  directly  very  recently.  Present 
j  estimates  are  directly  tied  to  observations  of  the  proteins 

produced.  When  the  proteins  that  are  produced  cause 
differences  observable  in  the  organism's  behavior,  a 
|  mutation  could  be  identified  which  caused  the 

change (Stansfield,  1969)  (Stryer,  1975). 

The  majority  of  mutations  are  lethal  to  the  organism 
and  result  in  the  death  of  the  organism.  The  mutant  gene, 
even  in  the  case  of  lethal  effects,  are  kept  at  a  low 
frequency  with  in  the  population.  The  effects  on  the  target 
population  and  the  frequency  of  the  allele  can  be  predicted 
by  modeling  techniques  used  in  population  genetics. 
Population  genetics  have  been  used  to  study  the  transfer  of 
advantagious  and  neutral  alleles  within  a  population. 
Mutant  types  produce  organism  which  compete  within  the 
ecosystem.  Natural  selection  appears  to  play  an  important 
role  in  the  level  of  an  allele  in  a  population.  To  date. 
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all  of  the  population  genetic  studies  that  have  been 
conducted  rely  on  phenotypic  manif estions.  Any  genotypic 
change  that  occurs  in  the  proteins  of  an  organism  can  only 
be  evaluated  using  present  techniques  when  the  genotypic 
change  causes  a  phenotypic  change.  The  hiearchy  of 
information  transmission  is  presented  in  figure 
2.2 (Stansf ield,  1969)  (Stryer,  1975). 

. - >  DNA  molecule 

I  I 

I  V 

j  .->  RNA  molecule 
I  I  I 

I  I  V 

' — ' - Protein  <-. 

I  I 

V  I 

. - >  Cell - 

I  I 

I  V 

' - Organism  < - . 

I  I 

V  I 

.->  Population  - 

I  I 

I  V 

' - Ecosystem  <-. 

I  I 

V  I 

Universe  - 

Figure  2.2 

Hieracrhy  Of  Biological  Information  Transfer 

There  a  several  types  of  mutations  possible.  Each  type 
may  be  identified  by  several  different  sets  of  terminology 
which  can  be  used  to  describe  non-exclusive  mutation  events. 
When  discussing  the  size  of  a  mutation,  the  two 
classifications  most  used  are  point  and  gross.  Point 
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mutations  are  usually  associated  with  changes  at  a  single 
nucleotide  pair.  Gross  mutations  involve  large  areas  of 
genetic  information.  This  could  describe  mutations  to 
multiple  nucleotides,  an  entire  gene,  or  even  entire 
chromosomes (Stansfield,  1 96 9)  (Stryer ,  1975). 

Mutations  can  also  be  classified  as  to  quality. 
Structural  mutations  affect  the  nucleotide  content  of  the 
gene.  Substitution  mutations,  deletion  mutations,  and 
insertion  mutations,  and  insertion  mutations  belong  in  this 
category.  Transition  mutations  occur  when  one  purine 
nucleotide  is  substituted  for  another  purine  or  when  one 
pyrimidine  is  substituted  for  another  pyrimidine.  Deletion 
mutations  occur  when  some  part  of  the  gene  is  lost  from  the 
gene.  Insertion  mutations  occur  when  nucleotides  are  added 
to  the  gene(Stansf ield,  1969). 

The  origin  of  the  mutation  is  also  a  source  of 
classification.  A  mutation  is  called  spontaneous  when  the 
origin  is  unknown.  This  is  the  category  where  the  random 
mutations  thought  to  be  important  in  evolutionary  drift. 
Genetic  control  mutations  are  thought  to  occur  when  the 
influence  of  a  specific  gene  affects  one  or  many  other 
genes.  Induced  mutations  occur  when  ionizing,  noionizing 
radiation,  and  chemical  mutagens  which  are  not  normally 
included  in  the  environment  affects  the  gene (Stansf ield, 
196  9)  . 

The  magnitude  of  the  phenotypic  effect  is  also  a 
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classification  scheme  for  mutation.  Some  alleles  can  be 
distinguished  by  the  frequency  with  which  they  mutate. 
Isoalleles  which  produce  seemingly  identical  pheotypes  when 
in  homozygous  or  heterozygous  combinations,  but  in 
combination  with  other  alleles  are  identifiable.  Mutants 
which  affect  the  viability  or  the  reproductive  vigor  of  the 
organism  can  be  identif ied(Stansf ield,  1969). 

Mutations  can  be  identified  by  the  direction  of  the 
phenotypic  response.  A  mutation  is  considered  forward  when 
the  normal  phenotypic  characteristic  suddenly  becomes  an 
abnormal  phenotype.  A  back  mutation  is  the  return  to  normal 
phenotype  after  an  abnormal  phenotype  was  observed.  Back 
mutations  have  been  observed  to  occur  either  by  a  change  at 
the  initial  location  within  the  gene  or  at  a  remote  site  in 
a  different  gene  (Stansf ield,  1969). 

Mutations  can  also  be  cataloged  as  to  the  type  of  cell. 
A  somatic  mutation  occurs  in  the  non-reproducti ve  cells  of 
the  body.  This  phenominon  is  seen  when  only  a  part  of  the 
organism  demonstrates  the  effect  of  the  mutation.  Gametic 
mutation  occurs  in  the  reproductive  cells  of  the  organism 
and  produces  an  heritable  change  seen  in  the  following 
generations  of  the  organism(Stansf ield,  1969). 

To  systematically  understand  the  relationships  involved 
in  the  mutation  of  an  organisms  DNA,  it  is  important  to 
understand  the  control  relationships  that  exist.  Figure  2.3 
examplifies  this  type  of  approach. 
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Figure  2.3 

Functional  Environmental  Relationships 


In  the  1960's  it  was  established  that  proteins 
reflected  the  genotypic  content  of  the  gene.  This 
reflection  was  found  to  be  common  across  all  forms  of  life 
according  to  the  genetic  code.  The  amino  acid  sequence  was 
the  result  of  a  specific  and  ordered  sequence  of  nucleotides 
in  DNA.  It  is  clear  that  these  proteins  represent  the  end 
products  of  molecular  evolution.  The  study  of  amino  acid 
sequences  has  enabled  a  new  perspective  on  the  study  of 
evolution.  These  sequences  represent  a  more  detailed  view 


of  the  phenotype  of  the  organism  than  is  possible  from  even 
the  most  detailed  observation  of  the  characteristics  of  the 
organism.  Atlas  Q±  Protein  Sequence  and  Stuictuie  edited  by 
Dayhoff  and  published  in  1972  provided  a  great  deal  of 
information  on  the  amino  acid  sequences  of  many  different 
proteins  and  many  different  species.  Because  a  large 
qunatity  of  purified  protein  and  the  massive  labor  efforts 
involved  in  doing  amino  acid  sequences,  the  amount  of 
information  which  can  be  found  is  still  pretty  sparse. 
Cytochrome  C  and  hemoglobin  are  the  two  proteins  which  have 
received  the  most  attention.  Sequences  from  many  species 
have  been  collected  for  these  two  proteins. 

The  second  generation  of  data  on  information 
macromolecules  has  already  begun.  New  techniques  for  the 
rapid  and  automatic  sequencing  of  very  small  quantities  of 
nucleic  acid  are  now  available.  Sequencing  of  mRNA  for  the 
hemoglobin  chains  is  now  available  for  use  in  this  study. 
The  sequencing  of  tRNA  have  been  completed.  Using  the 
information  from  the  structure  of  informational  molecules  is 
better  than  using  either  observable  characteristics  or  even 
amino  acid  sequencing.  Proteins  reflect  the  informational 
structure  of  the  gene  far  better  than  any  observable 
characteristic  of  the  organism.  The  chemistry  of  the  cell 
is  the  result  in  large  measure  to  the  functional  capability 
of  the  available  proteins.  Amino  acid  sequences  directly 
reflect  the  content  of  the  DNA.  A  single  nucleotide  change 
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can  result  in  the  placement  of  a  new  amino  acid  into  a 
protein.  This  amino  acid  replacement  may  or  may  not  affect 


any  observable  characteristic  of  the  protein.  A  single 
amino  acid  replacement  in  hemoglobin  causes  sickle  cell 
anemia.  However,  since  the  genetic  code  is  reduntant  with 
some  amino  acids  being  coded  by  as  many  as  6  codons,  many 
changes  in  nucleic  acid  result  in  no  change  to  the  primary 
protein  structure.  In  fact,  studies  show  that  the  structure 
of  the  genetic  code  is  such  that  most  of  the  changes  in 
nucleotide  content  either  produce  the  same,  or  a  chemically 
related  amino  acid.  Another  factor  to  be  considered  is  the 
existance  of  multiple  genetic  loci  coding  for  the  same 
protein.  This  can  result  in  the  production  of  multiple 
proteins  with  the  same  basic  structure  and  function.  The 
study  of  the  nucleic  acid  sequence  of  functional  proteins 
allows  a  more  in-depth  view  of  the  phenominon  at  work  at  the 
genetic  level  while  keeping  the  product  tied  to  the 
functional  structural  requirements  placed  on  the  protein 
chain. 

Use  of  a  thorough  understanding  of  the  structure- 
function  relationships  of  proteins  together  with  information 
available  from  analysis  can  result  in  an  understanding  of 
the  evolutionary  mechanism  that  is  far  better  than  what  is 
commonly  known  to  date.  Even  the  work  using  these 
techniques  that  have  been  accomplished  using  the  amino  acid 
sequence  information  has  the  increased  difficulty  associated 
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with  the  redundancy  of  the  code.  Understanding  evolutionary 
mechanics  at  the  level  of  the  protein  chain  does  not 
directly  reflect  the  process  at  the  level  of  nucleic  acid. 
For  instance,  the  colinearity  of  DNA  and  protein  has  been 
presumed  for  some  time  now.  Most  of  the  information  was 
derived  from  the  translation  of  mRNA  into  proteins.  When 
the  DNA  sequence  of  the  hemoglobin  gene  was  identified,  it 
was  found  to  have  a  large  region  of  non-coded  intervening 
sequences.  Also,  recently  the  discovery  of  large  portions 
of  repeated  sequences  of  DNA  identified  in  mammalian  DNA 
have  spurred  new  questions.  Study  of  the  evolutionary 
impact  of  random  mutation  and  the  effect  of  processes 
occurring  at  the  level  of  the  nucleic  acids  must  be 
conducted  using  relationships  at  that  level.  Using 
relationships  at  the  amino  acid  level  given  the  unique 
information  transmission  mechanism  of  the  genetic  code  is 
likely  to  give  a  different  perspective  of  reality. 

From  an  evolutionary  viewpoint,  the  primary  structure 
of  a  protein  is  in  large  part  the  end  point  of  a  biological 
phenomenon.  The  secondary  structure  and  binding  to 
cofactors  or  in  the  case  of  hemoglobin,  the  binding  of  the 
heme,  to  produce  a  functionally  proficient  molecule  within 
the  cell  is  the  direct  result  of  the  primary  sequence.  Any 
changes  in  the  primary  structure  which  alters  the 
function  of  the  protein  will  have  very  large  consequences  on 
the  survivability  of  the  cell,  and  the  organism.  On  the 
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other  hand,  many  changes  at  different  parts  of  the  protein 
result  in  little  or  no  change  in  the  capability  of  the 
protein  to  perform  its  biochemical  function. 

For  hemoglobin,  X-ray  crystallographic  determinations 
on  the  structure  of  hemoglobin  chains  from  many  species  show 
nearly  the  same  tertiary  and  quartinary  structure. 
Biochemical  differences  in  the  oxygen  transport  capability 
of  the  molecule  are  very  small.  Since  many  of  the  amino 
acid  positions  have  been  observed  to  have  different  residues 
that  result  in  changes  that  do  not  affect  function,  there 
has  been  considerable  debate  on  the  effect  of  neutral 
mutations  and  the  role  that  random  genetic  drift  may  have  on 
the  evolutionary  process.  Changes  that  are  not  the  result 
of  some  selective  advantage  appear  to  be  occuring  throughout 
the  proteins  of  the  cells.  Even  a  protein  with  the  survival 
requirements  that  are  tied  to  the  functional  demands  of 
hemoglobin  appear  to  have  a  large  number  of  residues  that 
may  be  changed  among  several  amino  acids.  Figure  3.2 
reflects  the  variety  of  residues  that  a  review  of  the 
literature  suggests  as  possible.  Of  note  at  this  time  is 
the  use  of  Figure  3.2.  That  these  amino  acids  have  been 
found  in  some  of  the  beta  chains  of  different  species  is  not 
to  say  by  any  stretch  of  the  imagination  that  all 
replacements  are  neutral.  In  fact,  many  of  these  changes 
have  been  associated  with  decreased  function  of  the 
hemoglobin.  Further,  many  of  these  residues  are  only 
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acceptable  when  other  residues  on  the  chain  have  already 
been  substituted.  Other  changes  are  undoubtedly  not 
acceptable  on  a  survival  basis  when  combined  with  other 
changes.  Because  of  the  limited  nature  of  detailed 
functional  and  survival  information  on  each  individual 
possible  change  the  simplifing  assumtion  made  in  this  study 
was  that  all  possible  combinations  would  be  allowed  in  this 
modeling  effort. 
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CHAPTER  III 


METHODOLOGY 

After  completing  the  literature  review,  it  was  obvious 
that  there  exists  strong  evidence  that  point  mutations  in 
the  amino  acid  sequences  of  several  different  proteins 
appear  to  have  significant  evolutionary  patterns  that  can 
be  studied.  The  two  most  commonly  characterized,  studied, 
and  analyzed  protein  chains  were  cytochrome  c  and 
hemoglobin.  Both  consist  of  protein  chains  longer  than  100 
amino  acids.  Both  have  had  extensive  studies  of  polymorphic 
variants  and  of  the  secondary,  and  tertiary  structures. 
Many  investigators  have  sequenced  both  compounds  from  many 
different  species.  Recent  work  on  sequencing  the  mRNA  of 
the  different  hemoglobin  chains  in  several  species  have 
enabled  this  study. 

The  target  sequence  of  the  DNA  represention  of  the  beta 
chain  of  hemoglobin  was  chosen.  The  unique  availability  of 
this  information  complemented  the  existing  data  on  point 
mutations  and  amino  acid  differences  among  and  between 
different  species.  Work  by  Dayhoff  and  by  Fitch  to 
construct  hemoglobin  ancestrial  trees  based  on  similarity 
studies  of  the  amino  acid  differences  between  species  and 
between  different  chains  of  hemoglobin  provided  insightful 
guidance. 
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For  this  study,  two  known  mRNA  sequences  of  beta 
hemoglobin  were  chosen,  that  of  human  and  that  of  rabbit.  A 
third  sequence  based  on  the  amino  acid  sequence  of  an 
hypothetical  ancestral  beta  hemoglobin  chain  which  should  be 
evolutionary  common.  Approximately  equal  numbers  of  point 
mutations  were  accepted  into  the  human  and  the  rabbit 
chains.  The  ancestral  beta  hemoglobin  chain  was 
characterized  and  included  in  the  modeling  effort.  The 
ancestral  beta  hemoglobin  I  selected  to  use  for  this  effort 
was  identified  as  locus  28  in  Dayhoff's  Atlas  of  Protein 
Sequence  and  Structure. 

Table  3.1  presents  the  DNA  representation  of  the  beta 
hemoglobin  chain  that  were  chosen  for  the  present  study. 
Very  few  point  mutations  exist  in  the  data  between  the  human 
and  the  rabbit  DNA.  Likewise,  there  are  few  point  mutations 
between  the  ancestral (ear ly)  DNA  and  either  the  human  or  the 
rabbit  DNA.  In  formulating  the  ancestral  nucleotide  string, 
minimizing  the  differences  between  existing  nucleotides  in 
the  human  and  rabit  DNA  was  a  goal.  The  ancestral  protein 
string  had  very  few  changes  in  amino  acids.  If  the  human 
and  rabbit  codon  at  any  particular  location  was  identical 
and  if  the  ancestral  amino  acid  was  that  same  codon,  then 
the  ancestral  amino  acid  was  assumed  to  have  the  same  codon 
as  the  other  two  DNA  chains.  If  the  human  and  rabbit  codons 
were  different,  then  the  codon  chosen  to  represent  the 
ancestral  codon  was  chosen  to  have  the  minimum  number  of 


nucleotide  differences  with  both  chains  and  still  code  for 


the  predicted  amino  acid.  If  both  the  human  and  the  rabbit 
codons  were  equidistant,  then  the  nearest  neighbor 
nucleotide  transition  rate  was  used.  The  nucleotide  chosen 
for  the  ancestral  codon  was  that  nucleotide  that  had  the 
highest  transition  probability  to  the  needed  nucleotide.  No 


equivalent  probabilities  were  encountered. 

Table  3.1 

Beta  Hemoglobin  DNA 


Amino  5,=>3'  Amino  5  *  =>3  *  Amino  5'=>3' 

Acid  Acid  Acid 


1 

Han  Rabbit 

Early 

1 

Han  Rabbit 

Early 

1  Man  Rabbit 

Early 

0 

CAT 

CAT 

CAT 

50 

AGT 

AGA 

AGA 

100 

AGG 

AGG 

AGG 

1 

CAC 

CAC 

CAC 

51 

AGG 

TGC 

AGC 

101 

CTC 

CTC 

CTC 

2 

GTG 

ATG 

ATG 

52 

ATC 

ATT 

ATC 

102 

GTT 

GTT 

GTT 

3 

CAG 

CAC 

CAC 

53 

TGC 

TGC 

TGC 

103 

GAA 

GAA 

GAA 

4 

AGT 

GGA 

GGT 

54 

AAC 

AAC 

AAC 

104 

CCT 

CCT 

CTT 

5 

AGG 

ACT 

AGC 

55 

CAT 

CAT 

CAT 

105 

GAG 

GAG 

GAG 

6 

CTC 

CTC 

CTC 

56 

GCC 

GTT 

GCC 

106 

CAG 

CAG 

CAG 

7 

CTC 

CTC 

CTC 

57 

GTT 

ATT 

ATT 

107 

GCC 

GCC 

GCC 

8 

CTT 

CTT 

CTT 

58 

AGG 

AGG 

AGG 

108 

GTT 

GTT 

GTT 

9 

AGA 

AGA 

AGC 

59 

CTT 

CTT 

CTT 

109 

CAC 

CAC 

CAC 

10 

CGC 

CGC 

CGC 

60 

CAC 

CAC 

CAC 

110 

CAG 

CAG 

CAG 

11 

AAC 

GAC 

AAT 

61 

CTT 

CTT 

CTT 

111 

GAC 

AAC 

AAC 

12 

AGT 

AGT 

AGT 

62 

AGC 

AGC 

AGC 

112 

ACA 

AAT 

AAT 

13 

GGC 

GGC 

GGA 

63 

ATG 

ATG 

ATG 

113 

CAC 

CAC 

CAC 

14 

CAG 

CAG 

CAG 

64 

GCC 

GCC 

GCC 

114 

CAG 

CAG 

CAG 

15 

CCA 

CAC 

CAC 

65 

CTT 

CTT 

CTT 

115 

GGC 

AGA 

AGC 

16 

GCC 

GCC 

GCC 

66 

TTT 

CTT 

TTT 

116 

ATG 

ATG 

AGG 

17 

CTT 

CTT 

CTT 

67 

CAC 

CAC 

CAC 

117 

GTG 

ATG 

ATG 

18 

CAC 

CAC 

CAC 

68 

GAG 

CAG 

CAG 

118 

AAA 

AAA 

AAA 

19 

GTT 

ATT 

ATT 

69 

ACC 

AGC 

ACC 

119 

GCC 

GCC 

GCC 

20 

CAC 

CAC 

CAC 

70 

GGC 

GGC 

GGC 

120 

TTT 

TTT 

TTT 

21 

ATC 

TTC 

TTC 

71 

AAA 

GAA 

AAA 

121 

TTC 

TTC 

TTC 

22 

TTC 

TTC 

TTC 

72 

ACT 

ACT 

ACC 

122 

GAA 

GAA 

GAA 

23 

AAC 

AAC 

AAC 

73 

ATG 

CTC 

ATC 

123 

GGT 

AGT 

AGT 

24 

ACC 

ACC 

ACC 

74 

GCC 

ACC 

ACC 

124 

TGG 

AGG 

AGG 

25 

ACC 

ACC 

ACC 

75 

CAG 

CAG 

CAG 

125 

TGG 

CTG 

TTC 

26 

CTC 

CTC 

CTC 

76 

AGC 

ACT 

TTT 

126 

CAC 

CAC 

CAC 

27 

GGC 

GGC 

GGC 

77 

GTG 

GTG 

GTG 

127 

CTG 

CTG 

CTG 

28 

CAG 

CAG 

CAG 

78 

CAG 

CAG 

CAG 

128 

AGC 

AGC 

AGC 

29 

GCC 

GCC 

GCC 

79 

GTC 

GTC 

GTC 

129 

GGC 

GGC 

GGC 

30 

CAT 

CCT 

CAT 

80 

GTT 

GTT 

GTT 

130 

ATA 

ATA 

CCA 

31 

CAG 

CAG 

CAG 

81 

GAG 

GAG 

GAG 

131 

CTG 

CTG 

CTG 

Table  3  .1  Continued 


Amino  5  *  => 3  *  Amino  5  *  =>3  *  Amino  5  *  => 3 ' 

Acid  Acid  Acid 


± 

Mil  Rabbit 

Early 

± 

Man  Rabbit 

Early 

i  Man  Rabbit 

Enxly 

32 

CAG 

CAG 

CAG 

82 

CTT 

TTT 

TTT 

132 

TTT 

CTT 

TTT 

33 

CAC 

AAC 

AAC 

83 

GCC 

GCC 

GCC 

133 

CAC 

CAC 

CAC 

34 

GAC 

GAC 

GAC 

84 

GGT 

GGT 

GTT 

134 

CAC 

CAC 

CAC 

35 

GTA 

GTA 

GTA 

85 

AAA 

AAA 

AAA 

135 

AGC 

AGC 

AGC 

36 

AGG 

TGG 

AGG 

86 

GGC 

AGC 

AGC 

136 

ACC 

ACC 

ACC 

37 

CCA 

CCA 

CCA 

87 

TGT 

CTT 

TTT 

137 

CAC 

CAC 

CAC 

38 

GGT 

GGT 

GGT 

88 

CAG 

CAG 

CAG 

138 

AGC 

GGC 

AGC 

39 

CTG 

CTG 

CTG 

89 

ACT 

ACT 

ACT 

139 

ATT 

ATT 

ATT 

40 

CCT 

CCT 

CCT 

90 

CTC 

TTC 

TTC 

140 

GGC 

GGC 

GGC 

41 

GAA 

GAA 

GAA 

91 

CAG 

CAG 

CAG 

141 

CAG 

CAG 

CAG 

42 

AAA 

GAA 

AAA 

92 

GTG 

GTG 

GTG 

142 

GGC 

AGC 

AGC 

43 

CTC 

CTC 

GTC 

93 

ACA 

ACA 

ACA 

143 

GTG 

GTG 

GTG 

44 

GGA 

GGA 

GGA 

94 

GTC 

GTC 

GTC 

144 

CTT 

TTT 

TTT 

45 

AAA 

AAA 

AAA 

95 

CTT 

CTT 

CTT 

145 

ATA 

GTA 

ATA 

46 

CCC 

CCC 

CCC 

96 

CAG 

CAG 

CAG 

146 

GTG 

GTG 

GTG 

47 

ATC 

GTC 

ATC 

97 

GTG 

GTG 

GTG 

147 

TTA 

TTA 

TTA 

48 

CAG 

CAG 

CAG 

98 

CAC 

CAC 

CAC 

49 

GGA 

GGA 

GGA 

99 

ATC 

ATC 

ATC ( BARRALLE , 

1977A 

s 1093) 

The  model  was  set  up  so  that  all  three  DNA  sequences 
were  run  with  the  same  random  number  stream  generating  point 
mutations  for  each  run  of  the  model.  The  complete  model  can 
be  found  in  Appendix  A  of  this  report. 

Potential  point  mutations  were  assessed  by  comparison 
of  the  random  number  with  the  assigned  probability  of  a 
mutation  occurring.  The  probability  of  a  mutation  occurring 
at  any  nucleotide  was  dependent  upon  the  structure  of  the 
model  and  the  selection  of  the  data  coefficient  chosen. 
Three  structures  were  tested  for  creating  point  mutation. 
The  first  technique  is  called  the  random  method.  If  the 
random  number  representing  the  probability  of  having  a 
point  mutation  is  below  the  threshold  established  for  having 


then  the  mutation  was  tested  for 


a  mutation , 

validity/reasonableness.  Each  nucleotide  was  assigned  an 
equal  probability  of  becoming  changed  to  any  of  the  other 
three  nucleotides  and  this  was  the  point  mutation  tested. 
Point  mutation  testing  consisted  of  evaluating  whether  or 
not  the  amino  acid  that  would  result  from  the  new  codon  was 
possible  at  that  position.  Possibility  was  arbitrarily 
estabilished  in  this  effort  as  having  the  amino  acid  in 
question  at  that  position  recorded  either  as  a  known  point 
mutation,  or  protein  polymorph,  or  even  an  amino  acid  found 
in  one  of  the  other  species  for  which  beta  hemoglobin 
protein  sequenceing  had  been  conducted.  See  Table  3.2  for  a 
compiled  list  of  acceptable  amino  acids  at  each  position  on 
the  chain. 


Table  3.2 

Allowable  Amino  Acid  by  Position 
Beta 

Hemoglobin  Amino 


Beta 

Hemoglobin  Amino 

Eflsition  Acid 

1  T,V,M,P,G 

2  H,L,Q,E,D,N,Y,R 

3  L,K,F,W 

4  T  S  A 

5  p'dJs,G,E,A,K,N 

6  E,A,D,G,K,Q,P,V,S 

7  E,A,G,D,N,Q,K 

8  K ,  D ,  S 

9  S,A,N,T,V,D,C,H,G,E 

10  A,T,H,L,E,Q,S 

11  V,G, I 

12  T,S,N,A,L,D,R,I,V,K 

13  A,G,T,S,C,N 

14  L,S,F,R 

15  W,F ,G,S 


Positi.Qn  Acid 

74  G,A,D 

75  L,V,M, I ,T 

76  A,N,S,T,H,K,G,Q,E 

77  H,N,E,Q,D 

78  L,  P 
D,N,E 
N,D,K,S 
L,I 
K 

G,A,N 
T,A,H 
F,Y 

86  A,  S 

87  T,S,K,E,A,Q,N,H 

88  L,R,P 


79 

80 
81 
82 

83 

84 

85 


Beta 

Hemoglobin  Amino 
Position  Acid 

17  K,I,E,H 

18  V,K,I,H 

19  N,K,H,G, 

20  V,h,S,E, 

21  D,EfH,A, 

22  EfDfRfAj 

23  V,DfAfT, 

24  G,L,R,V 

25  Gf K,Af R 

26  EfQ,YfK 

27  A,H,T 

28  Lf  P 

29  Gf  A 

30  RfS 

31  L 

32  L 

33  V,L,I 

34  V 

35  Y/F 

36  P 

37  W,S 

38  T 

39  EfR,Q,S 

40  R 


Table  3.2  Continued 
Beta 

Hemoglobin 


*  F  * 

FfS 

EfSfDfQ/T/AfR 
SfHfTf A, Ef Q 

FfL 
Gf  E 
DfA,N 
L 

SfG 

TfSfNfD,K 

PfA 

Df  Nf  Sf  HfGfKfE,i 

A 

Vfl 

MfLfC 

GfNfSfOfAfH 

N,D,A 

PfA,R 

K,E,TfM,Q 

V 

Rf LfNfEfR 


n  Amino 
Acid 
EfQ,R 
LfY,P 
H,  Y 

C, S 

D, V,N 

R  f  A  f  E  f  Q 
L 

HfE,R 

V,M 

D,N,H,Y 

P 

Ef  Qf  A 

NfR,TfD 

F 

RfR,N 
Lf  R 
L 
G 

Nf  D 

VfM,S,I,A 

LflfFfV 

VfDfNfA,I 

CfLflfSfVfTfH 

V,CfEfA 

L 

AfGfS 

HfDfN,RfIfE 

H,N,R 


RfD,NfSfHf E 
EfDfQfNfEfQ,K 

F 

T,Df  N,S 


V,AfL,M,T,EfC 

QfErL 

AfV,H,S 

AfSfD,QfEfGfN 

Y,FfH,Wf LfD 
Qf  Ef  L 
RfAfQ,E 
VfM,LfH 
V ,  F 

AfTfCfSfR 

G,A,D,V 


Table  3.2  Continued 


Beta 

Hemoglobin  Amino 

Beta 

Hemoglobin  Amino 

Position 

Acid 

Position 

Acid 

64 

G 

137 

V 

65 

K ,  S  r  A,  E 

138 

A,G,S 

66 

K,E 

139 

N,A,D,T,S,H 

67 

V,A,D,E 

140 

A 

68 

LrI 

141 

L,R 

69  P,R 

,G,T,A,S,D,N,H,V,E,Q 

142 

A,S,G 

70 

A,T,S 

143 

H,S,K,R 

71 

FrL,I,S 

144 

K,R,A 

72 

S/C/G/AfD/K 

145 

y,c,h 

73 

D,E,Q,N 

146 

H,D 

The  second  structure  tested  relied  on  data  compiled 
from  three  different  data  sources  to  address  the  frequency 
of  any  nucleotide  being  replaced  by  mutation  to  the  other 
three.  These  estimate  were  derived  from  different 
collections  of  proteins.  This  structure  results  in  three 
assessments  of  nucleotide  mutation  potential  and  amino  acid 
possibility  assessment.  The  transition  probabilities  are 
listed  in  Table  3.3. 


Table  3.3 

DNA  Nucleotide  Change  Frequency 


Transition 

Dayhof f 

Zuckerkandl 

Fitch 

frequency 

frequency 

fx.eguency 

A=>C 

.0244 

.0610 

.0113 

A**>G 

.0731 

.0360 

.0000 

A»>T 

.0488 

.0690 

.0602 

C=>A 

.0122 

.0850 

.0458 

C=>G 

.0854 

.0  980 

.0989 

C=>T 

.2927 

.1590 

.0960 

G=>A 

.0610 

.02  80 

.0277 

G=>C 

.0976 

.0920 

.0635 

G=>T 

.0732 

.0890 

.1995 

T=>  A 

.0366 

.0580 

.1140 

T=>C 

.1098 

.1210 

.1237 

T=>G 

.0854 

.0940 

.1694 
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The  third  structure  conformed  to  a  first  order  Ising 
model  in  which  the  ordered  nucleotide  composition  of  the  two 
nearest  neighbors  determined  what  the  transition  probability 
at  any  particular  nucleotide  site  on  the  genome.  Again, 
amino  acid  possibility  assessment  was  made  before  any 
nucleotide  change  was  allowed  to  occur.  The  transition 
probabilities  are  listed  in  Table  3.4. 


Table  3  .4 

DNA  Nearest  Neighbor  Nucleotide  Change  Frequency 


h=±&  h=>.Q 

A_A  .9946  .00016 

A_C  .9901  .00198 

A_G  .9957  .00217 

A_T  .9872  .00230 

C_A  .9948  .00035 

C_C  .9901  .00194 

C_G  .9947  .00198 

C_T  .9874  .00342 

G_A  .9946  .00135 

G_C  .9900  .00198 

G_G  .9947  .00217 

G_T  .9872  .00199 

T_A  .9947  .00035 

T_C  .9901  .00194 

T_G  .9947  .00198 

T_T  .9872  .00259 

5  1  ->3  *  G=>A  S=2£ 
.00151  .01389 
A_C  .00692  .00053 
A_G  .00185  .00238 
/UT  .00176  .00821 
C_A  .00295  .01263 
C_C  .00640  .00049 
C_G  .00167  .00215 
C_T  .00165  .00660 
G_A  .00205  .01231 
G_C  .00692  .00053 
G_G  .00185  .00228 
G_T  .00183  .00  87  9 
T_A  .00215  .01354 
T_C  .00640  .00049 
T_G  .00167  .00215 
T_T  .00227  .00765 


A=>G  C=>A 

00175  .00349  .00018 
00713  .00079  .00061 
00120  .00193  .00094 
00722  .00328  .00125 
00247  .00247  .00072 
00699  .00097  .00063 
00110  .00221  .00093 
00315  .00603  .00032 
00135  .00270  .00000 
00713  .00079  .00061 
00120  .00193  .00094 
00796  .00284  .00125 
00247  .00247  .00032 
00699  .00097  .00063 
00110  .00221  .00093 
00484  .00536  .00034 
fi=2».G  ,G=>T  T=>-A 

9840  .00060  .00311 

9867  .00585  .00101 

9926  .00317  .00152 

9871  .00293  .00137 

9840  .00042  .01333 

9867  .00640  .00053 

9926  .00358  .00157 

9871  .00465  .00072 

9840  .00164  .00467 

9867  .00585  .00101 

9926  .00317  .00152 

9871  .00228  .00123 

9840  .00031  .01333 

9867  .00650  .00101 

9926  .00358  .00231 

9871  .00298  .00072 


C=>C  C=>G  C=>T 
9978  .00180  .00024 

9935  .00427  .00163 

9913  .00126  .00650 

9840  .00847  .00627 

9976  .00168  .00000 

9935  .00440  .00147 

9913  .00124  .00652 

9933  .00260  .00378 

9936  .00260  .00378 

9935  .00427  .00163 

9913  .00126  .00650 

9840  .00847  .00627 

9976  .00208  .00000 

9650  .00440  .00147 

9913  .00124  .00652 

9913  .00394  .00442 

T=>C  T=>G 
00206  .00033  .9945 
00730  .00579  .9859 
00532  .00196  .9912 
01164  .00479  .9822 
00000  .00263  .9840 
00904  .00393  .9865 
00876  .00207  .9876 
00417  .00251  .9926 
00050  .00033  .9945 
00730  .00579  .9859 
00532  .00196  .9912 
01043  .00614  .9822 
00000  .00267  .9840 
00504  .00746  .9865 
00707  .00303  .9876 
00417  .00251  .9926 
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When  the  possibility  assessment  was  favorable,  the 
transition  is  recorded  and  the  new  nucleotide  will  be 
substituted  into  the  chain.  When  the  possibility  assessment 
did  not  reveal  evidence  for  the  particular  amino  acid 
appearing  at  that  position  of  the  hemoglobin  chain,  the 
nucleotide  will  be  left  without  change  and  the  failure  to  be 
able  to  change  the  nucleotide  is  recorded. 

Many  types  of  information  can  be  recorded  from  the 
model.  In  addition  to  writting  out  checks  to  validate  the 
input  data,  the  model  can  be  used  to  write  out  information 
to  trace  the  point  mutations  that  occurred  sequentially  from 
the  start  of  the  experiment.  The  information  on  transition 
rates  for  both  amino  acids  and  nucleic  acids  is  output  as 
well  as  the  number  of  times  that  possibility  assessment  did 
not  allow  the  transition.  The  number  of  times  that  the 
mutation  of  an  nucleic  acid  or  an  amino  acid  back  mutates  to 
an  earlier  configuration  at  a  location  along  the  chain. 
Other  recorded  data  includes  the  simularity  of  the  resulting 
amino  acid  sequence  with  the  starting  protein  chains  of 
human,  rabbit,  and  ancestrial  beta  hemoglobin.  The  amino 
acid  and  nucleotide  sequences  entered  at  each  location  was 
also  put  into  data  files  for  storage.  This  information  was 
stored  into  intermediate  files  to  be  used  by  the  Statistical 
Package  for  Social  Scientists (SPSS)  for  data  analysis  and 
computations . 
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Statistical  tests  were  used  to  evaluate  the  results  of 
the  simulation.  Linear  regression  and  analysis  of  variance 
were  used  to  help  understand  the  relationships  that  were 
revealed  during  the  experiment.  Review  and  categorization 
of  the  results  was  used  to  examine  relationships  and  to 
determine  the  applicability  of  different  techniques  and  the 
implications  of  the  results  of  this  simulation. 


CHAPTER  IV 


RESULTS 

The  simulation  model  used  for  this  study  was 
implemented  and  tested  on  the  Cyber  computer  system  used  by 
the  School  of  Engineering  at  the  Air  Force  Institute  of 
Technology.  Model  verification  was  conducted.  Data 
collection  consisted  of  1275  runs  of  the  model  on  the  IBM 
370/155  operated  by  the  Human  Engineering  Division  of  the 
Air  Force  Aerospace  Medical  Research  Laboratory.  425  runs 
of  the  model  were  made  for  each  of  the  three  initial 
hemoglobin  DNA  sequences  to  be  examined.  The  number  of 
iterations  selected  were  randomly  determined  out  of  three 
groups.  A  stream  of  random  digits  between  0  and  9  were  used 
to  establish  the  number  of  iterations (probability  assessment 
of  a  mutation  at  each  position  along  the  entire  chain)  to  be 
used.  The  three  groups  were  formed  by  sequentially 
multiplying  the  sum  of  the  chosen  random  digit  and  1  by  1, 
10,  and  1000.  This  resulted  in  a  stream  of  iterations  which 
varied  between  1  and  1000.  Random  number  streams  were 
generated  within  the  model  for  tl—  probabal is  tic 
determination  of  a  mutation  occuring  at  any  location. 

Statistical  testing  of  the  random  digit  stream  revealed 
no  trends  or  bias.  Verification  included  of  examining  and 


comparing  the  resultant  nucleotide  transition  frequencies 
with  the  transition  tables  for  each  of  the  methods  used  in 
the  model.  There  was  no  observable  abnormality  and  the  Chi 
Square  test  revealed  no  significant  difference  at  the  0.05 
level.  The  pattern  of  changes  in  amino  acids  produced  in 
the  model  appeared  to  be  consistent  with  the  pattern  of 
mutations  observed  by  phenotypic  studies  of  abnormal 
hemoglobin  in  vivo.  The  amino  acid  mutations  produced  by 
the  model  resulted  in  primary  sequences  indistinguishible 
from  primary  sequences  observed  in  different  species. 


Table  4.1 

Nucleotide  Transition  Anova 


street 

igmiicance 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Hemoglobin 

<0.01 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

Interations 

<0.001 

Method 

<0.001 

The  results  of  the  anova  used  to  examine  the  changes  in 
nucleotides  observed  in  the  model  can  be  found  in  table  4.1. 

The  anova  revealed  that  both  the  method  being  used  to 
generate  transition  probabilities  and  the  number  of 
iterations  of  the  model  has  a  consistently  significan  effect 
on  the  transition  of  a  nucleotide  to  one  of  the  other  three 
possible  nucleotides.  Hemoglobin  starting  chain  was  only 
found  significant  for  the  transition  of  guanine  to  adenine. 

The  results  of  the  anova  conducted  to  examine  the 
transitions  between  amino  acids  can  be  found  in  table  4.2. 


Table  4.2 

Amino  Acid  Transition  Anova 


I 


I 


Transition 

Effects 

(Significance  Level) 

Interations  Method 

Hemoglobi 

G=>A 

<0.001 

<0.001 

G=>S 

<0.001 

<0.001 

G=>V 

<0.001 

<0.001 

G=>E 

<0.001 

<0.001 

<0.001 

G=>D 

<0.001 

<0.001 

<0.05 

G=>R 

<0.001 

<0.001 

G=>C 

<0.001 

<0.001 

G=>W 

<0.001 

A=>G 

<0.001 

<0.001 

A=>S 

<0.001 

<0.001 

A=>V 

<0.001 

<0.001 

A=>K 

<0.001 

<0.001 

A=>T 

<0.001 

<0.001 

A=>E 

<0.001 

<0.001 

<0.05 

A=>P 

<0.001 

<0.001 

<0.001 

A=>D 

<0.001 

<0.001 

<0.001 

A=>R 

<0.001 

<0.001 

<0.001 

A=>M 

<0.001 

<0.001 

<0.001 

L=>S 

<0.001 

<0.001 

L=>V 

<0.001 

<0.001 

L=>P 

<0.001 

<0.001 

<0.001 

L=>I 

<0.001 

<0.001 

L=>R 

<0.001 

<0.001 

L=>Q 

<0  .001 

<0.001 

L=>F 

<0.001 

<0.001 

L*>H 

<0.001 

<0.001 
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Table  4.2  Continued 


Transition 

Effects 

(Significance  Level) 

Interations 

Hemoglobi 

L=>M 

<0.001 

<0.001 

<0.001 

L=>W 

<0.001 

<0.001 

S=>G 

<0.001 

<0.001 

S=>A 

<0.001 

<0.001 

S=>L 

<0.001 

<0.001 

S=>T 

<0.001 

<0.001 

S=>P 

<0.001 

<0.01 

<0.001 

S  =  >I 

<0.001 

<0.001 

S=>R 

<0.001 

<0.001 

S=>N 

<0.001 

<0.001 

<0.001 

S=>F 

<0.001 

<0.001 

S=>Y 

<0  .001 

<0.00. 

S*>C 

<0.001 

<0.001 

s=>w 

<0  .001 

<0.001 

V=>G 

<0.001 

<0.001 

V=>A 

<0.001 

<0.001 

V=>L 

<0.001 

<0.001 

V»>E 

<0.001 

<0.001 

<0.05 

V=>D 

<0.001 

<0.001 

V=>I 

<0.001 

<0.001 

<0  .001 

V=  >F 

<0.001 

<0.001 

V=>M 

<0.001 

<0.01 

K=>A 

<0.001 

<0.001 

K=>T 

<0.001 

<0.001 

K=>E 

<0.001 

<0.001 

K=>I 

<0.001 

<0.001 

<0.05 

K=>R 

<0.001 

<0.001 

K=>N 

<0.001 

<0.001 

K=>Q 

<0.001 

<0.001 

K=>M 

<0.001 

<0.001 

T*>A 

<0.001 

<0.001 

T=>S 

<0.001 

<0.001 

<0.001 

T=>K 

<0.001 

<0.001 

<0.001 

T*>P 

<0.001 

<0.001 

<0.05 

T=>I 

<0.001 

<0.001 

T=>R 

<0.001 

<0.001 

T=>N 

<0.001 

<0.001 

F=>G 

<0  .001 

<0.001 

F*>A 

<0  .001 

<0.001 

F*>V 

<0 .001 

<0.001 

F*>K 

<0.001 

<0.001 

F*>D 

<0.001 

<0.001 

<0.001 

F*>Q 

<0.001 

<0.001 

P*>A 

<0.001 

<0.001 

<0.001 

P*>L 

<0.001 

<0.001 

P=>S 

<0.001 

<0.001 

<0.001 

P«>T 

<0.001 

<0.001 

P*>R 

<0.001 

<0.001 
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Table  4.2  Continued 


Transition 

Effects 

(Significance  Level) 

Interations  Method 

Hemoglobi 

P=>Q 

<0.001 

<0.001 

P*>H 

<0  .001 

<0.001 

D=>G 

<0.001 

<0.001 

D=>A 

<0.001 

<0.001 

D=>V 

<0.001 

<0.001 

D=>E 

<0.001 

<0.001 

D«>N 

<0.001 

<0.001 

D*>Y 

<0.001 

<0.001 

D»>H 

<0.001 

<0.001 

I*»>L 

<0.001 

<0.001 

I»>S 

<0.001 

<0.001 

<0.05 

I*>V 

<0  .001 

<0.001 

1->K 

<0.001 

<0.001 

<0.001 

I*>T 

<0 .001 

<0.001 

I*>N 

<0.001 

<0.001 

I->F 

<0.001 

<0.001 

I«>M 

<0.001 

<0.01 

R«>G 

<0.001 

<0.001 

<0.05 

R«>A 

<0.001 

<0.001 

R»>L 

<0.001 

<0.001 

R*>S 

<0.001 

<0.001 

R«>K 

<0.001 

<0.001 

<0.05 

R»>T 

<0.001 

<0.001 

R«>P 

<0.001 

<0 .001 

R»>I 

<0.001 

<0.01 

R«>Q 

<0.001 

<0.001 

<0  .01 

R«>C 

<0.001 

<0.001 

<0.05 

R=>H 

<0  .001 

<0.001 

R*>M 

<0.001 

<0.001 

R=>W 

<0.001 

<0.001 

N*>S 

<0.001 

<0.001 

N=>K 

<0  .001 

<0.001 

Nss>T 

<0.001 

<0.001 

<0.05 

N»>D 

<0.001 

<0.001 

N«>I 

<0.001 

<0.001 

N»>Y 

<0.001 

<0.05 

N=>H 

<0.001 

<0.001 

Q»>L 

<0  .001 

<0.001 

Q=>K 

<0.001 

<0.001 

Q->E 

<0.001 

<0.001 

Q=>P 

<0.001 

<0.001 

<0.001 

Q=>R 

<0.001 

<0.001 

Q=>H 

<0.001 

<0.001 

F  =  >L 

<0.001 

<0.001 

F=>S 

<0.001 

<0.001 

F =>V 

<0.001 

<0.001 

F=>I 

<0.001 

<0.001 

F  =  >Y 

<0  .001 

<0.001 
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Table  4.2  Continued 


Transition 

Effects 

(Significance 

:  Level) 

Interatior 

L&  Method 

imoclobi 

F=>C 

<0.001 

<0.001 

<0.05 

Y=>S 

<0.001 

<0.001 

Y=>D 

<0.001 

<0.001 

Y=>N 

<0.001 

<0.001 

Y=>F 

<0.001 

<0.05 

Y=>C 

<0.001 

<0.001 

<0.001 

Y=>H 

<0.001 

<0.001 

C=>G 

<0.001 

<0.001 

C=>S 

<0.001 

<0.001 

<0.001 

C=>R 

<0.001 

C=>F 

<0.001 

<0.001 

C=>Y 

<0.001 

<0.001 

H=>L 

<0.001 

<0.001 

H=>P 

<0.001 

<0.001 

H=>D 

<0.001 

<0.001 

H=>R 

<0.001 

<0.001 

H=>N 

<0.001 

<0.001 

H=>Q 

<0.001 

<0.001 

H=>Y 

<0.001 

<0.001 

M=>A 

<0.001 

<0.001 

M=>L 

<0.001 

<0.001 

M=>V 

<0.001 

<0.001 

M=>K 

<0.001 

<0.001 

M=>I 

<0.001 

M=>R 

<0.001 

<0.01 

<0.05 

W=>G 

<0.001 

W=>L 

<0.001 

<0.001 

W=>S 

<0.001 

<0.001 

W=>R 

<0.001 

<0.001 

As  is  seen  in  examining  the  results  of  the  nucleotide 
anova,  both  the  number  of  iterations  of  the  chain  and  the 
method  of  assessing  the  transition  probabilities  resulted  in 
overwhelmingly  significant  effects.  Hemoglobin  was  found  to 
be  a  significant  effect  for  several  of  the  amino  acid 
transitions. 

Table  4.3  presents  information  on  the  frequency  of  the 
transitons  observed  in  the  position  and  iteration  averaged 
values  of  the  nucleotide. 


Table  4.3 


Average  Nucleotide  Transition  Frequency 

Nearest- 

Random  Dayhoff  Zuckerkandl  Fitch  Neighbor 


Allowed 

Change 

49.43 

117.90 

108.75 

111.55 

43.94 

Repeated 

Change 

6.70 

38.45 

36.27 

38.03 

6.34 

Rejected 

Change 

211.58 

542.46 

526.98 

530.59 

197.91 

Table  4.4  contains  the  results  from  the  regression  of 
the  value  of  the  positioned  averaged  transition  of 
nucleotides. 


Table  4.4 

Average  Nucleotide  Transition  Regression 


Intercept/ 

Slope 

Random 

Dayhoff 

Zuckerkandl  Fitch 

Nearest- 

Neighbor 

Allowed 

1.311 

3.144 

5.296 

3.434 

2.236 

Change 

.2254 

.5375 

.4846 

.5064 

.1953 

Repeated 

-1.515 

-7.924 

-6.721 

-6.073 

-1.454 

Change 

.0385 

.2172 

.2013 

.2108 

.0365 

Rejected 

-2.160 

0.000 

0.396 

-2.531 

-2.792 

Change 

1.160 

2.944 

2.858 

2.894 

1.089 

In  these  tables,  the  variable  that  tracks  allowed 
change  is  incremented  every  time  that  a  nucleotide  change 
produces  a  reasonable  amino  acid  at  the  particular  amino 
acid  position  coded  for  by  the  nucleotide.  Repeated  change 
refers  to  any  nucleotide  change  for  which  the  nucleotide 
previously  had  the  same  base  at  some  earlier  time  during  the 
run  of  the  model.  Rejected  change  is  incremented  whenever 
the  proposed  nucleotide  transition  would  result  in  the 
coding  of  an  amino  acid  residue  that  has  not  been  previously 
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associated  with  any  of  the  polymorphs,  abnormal  hemoglobin 
variants,  or  species  with  sequenced  hemoglobin  in  Atlas  of 


Protein  Sequence  and  Structure  at  that  particular  location. 

Table  4.5  presents  information  on  the  frequency  of  the 
transitons  observed  in  the  position  and  iteration  averaged 
values  of  the  amino  acid  residues. 


Table  4.5 

Average  Amino  Acid  Transition  Frequency 

Nearest- 

Random  Dayhoff  Zuckerkandl  Fitch  Neighbor 

Accepted  42.38  105.40  97.97  100.45  41.12 

Change 

Repeated  49.42  117.87  108.75  111.55  43.94 

Amino  Acid 


Table  4.6  contains  the  results  from  the  regression  of 
the  value  of  the  positioned  averaged  transition  of  amino 
acid  residues. 


Table  4.6 

Average  Amino  Acid  Transition  Regression 


Intercept/ 

Slope  Random  Dayhoff  Zuckerkandl  Fitch 
Accepted  1.837  1.972  4.527  2.772 
Change  .2177  .5553  .5017  .5244 
Repeated  1.311  3.144  5.256  3.434 
Amino  Acid  .2254  .5375  .4846  .5064 


Nearest- 
Neighbor 
2.065 
.20  97 
2.236 
.1953 


In  the  two  preceeding  tables,  accepted  changes  are  the 


result  of  nucleotide  changes  that  were  compatible  with 
coding  for  accepted  amino  acids  at  the  appropriate  position. 
The  repeated  amino  acid  line  of  the  tables  presents  that 
subset  of  the  accepted  amino  acid  changes  that  result  in  an 
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amino  acid  that  was  previously  found  during  a  particular 
simulation  run  at  the  same  position  of  the  protein  chain. 

Table  4.7  contains  a  listing  of  the  regression  of 
transition  frequency  as  a  function  of  interation  listed  by 
nucleotide  position. 


Table  4.7 

Nucleotide  Transition  Regression 


Intercept/ 

*  « 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

»  "  H 

1 

None  Seen 

No  Repeats 

None  Rejected 

v. 

2 

0  .0240 

-0.0160 

0.0679 

k 

0.0015 

0.0004 

0.0056 

3 

-0  .0291 

-0  .0287 

0.1318 

/\\" 

0.0016 

0.0007 

0.0055 

4 

-0.0094 

-0.0033 

-0.0235 

0.0001 

0.0000 

0.0025 

i 

5 

-0.0410 

-0.0259 

0.1093 

0.0015 

0.0006 

0.0058 

6 

0  .0438 

-0.0035 

-0.2064 

0.0018 

0.0002 

0.0054 

7 

None  Seen 

No  Repeats 

0.1686 

•  \ 

0.0022 

■ 

8 

-0.0214 

-0.0014 

0.0231 

0.0008 

0.0000 

0.0059 

9 

-0.0260 

-0.0125 

0.0624 

.V 

0.0004 

0.0002 

0.0060 

10 

-0.0033 

-0.0004 

0  .0386 

*3""* 

0  .0000 

0 .0000 

0  .0034 

11 

0  .0463 

-0.0016 

0  .0294 

0p 

0.0009 

0.0000 

0.0081 

,*• 

12 

0.0039 

0.0016 

-0.0365 

»•, 

0.0002 

0 .0001 

0.0076 

>•- 

13 

-0.0000 

No  Repeats 

-0.0456 

0.0000 

0.0015 

14 

0.0608 

-0.0176 

-0.1250 

0.0011 

0.0005 

0.0066 

15 

0.0999 

0.0016 

-0.0358 

»  * 

0.0012 

0.0001 

0.0056 

►  \ 

16 

0  .0885 

-0.0000 

-0.0540 

■  ", 

0.0008 

0.0000 

0.0028 

17 

-0  .0087 

-0  .0074 

-0 .0406 

W 

±'m 

0.0004 

0.0002 

0.0050 

40 

$ 

L '■V' v*w!*-v-v  v  .  ...  .  . 

J  J 
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Intercept/ 

Table  4.7 

Slope 

Allowed 

Position 

Change 

18 

0.0653 

0.0005 

19 

-0.0018 

0.0010 

20 

0.1507 

0.0011 

21 

-0.0056 

0.0018 

22 

-0.1228 

0.0030 

23 

-0.0120 

0.0034 

24 

-0.0360 

0.0008 

25 

-0.0072 

0.0001 

26 

0.0145 

0.0014 

27 

-0.0066 

0.0009 

28 

None  Seen 

29 

0.0470 

0.0010 

30 

0.0580 

0.0004 

31 

None  Seen 

32 

0.0382 

0.0006 

33 

-0.0041 

0.0006 

34 

None  Seen 

35 

None  Seen 

36 

None  Seen 

37 

None  Seen 

38 

0.0010 

0.0001 

39 

-0.0454 

0.0029 

40 

None  Seen 

41 

-0.0200 

0.0007 

Continued 


Repeated 

Rejected 

Change 

Change 

0.0018 

-0.1654 

0.0001 

0.0068 

-0  .0100 

-0  .0092 

0.0005 

0.0011 

-0.0000 

-0.0907 

0.0000 

0  .0068 

-0.0084 

0.1703 

0.0010 

0.0066 

-0.0770 

0.0410 

0.0015 

0.0018 

-0.0  818 

0.1544 

0.0017 

0.0076 

-0.0111 

-0.0565 

0.0003 

0.0098 

-0.0012 

-0.0216 

0.0000 

0.0008 

-0.0117 

-0.0622 

0.0006 

0  .0046 

-0.0056 

0.0551 

0.0002 

0  .0041 

No  Repeats 

-0.0043 

0.0022 

-0.0080 

-0.1331 

0.0002 

0.0062 

No  Repeats 

-0.1106 

0.0064 

No  Repeats 

-0.0449 

0.0022 

0.0004 

0.0196 

0.0000 

0.0062 

-0.0087 

-0.0257 

0.0003 

0.0065 

No  Repeats 

-0.0902 

0.0021 

No  Repeats 

0.0538 

0.0056 

No  Repeats 

-0  .1616 
0.0077 

No  Repeats 

-0.0440 

0.0009 

-0.0006 

-0.1044 

0.0000 

0.0066 

-0.0504 

0.0850 

0  .0013 

0.0045 

No  Repeats 

-0.0331 

0.0007 

-0.0086 

-0.0003 

0.0002 

0  .0045 
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Table  4.7  Continued 


Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

42 

-0.0136 

-0  .0366 

-0.0745 

0.0025 

0.0012 

0.0054 

43 

None  Seen 

No  Repeats 

-0.0008 

44 

0.0485 

-0.0188 

0.0059 

0.1758 

0.0018 

0.0007 

0.0067 

45 

-0.0143 

-0.0099 

0.0086 

0.0009 

0.0003 

0.0062 

46 

-0.0230 

-0.0179 

0.0039 

0.0008 

0.0004 

0.0013 

47 

0.0643 

-0  .0200 

0.0545 

0.0021 

0.0010 

0.0069 

48 

-0.0003 

-0.0084 

0.1020 

0.0014 

0.0008 

0  .0052 

49 

None  Seen 

No  Repeats 

-0.0427 

50 

-0.0459 

-0.0284 

0  .0047 
0.1272 

0.0014 

0.0006 

0.0062 

51 

0.0556 

No  Repeats 

-0.1116 

52 

0.0011 

-0.0207 

-0.0015 

0.0058 
-0  .0292 

0.0007 

0.0002 

0.0006 

53 

-0.0170 

-0.0166 

-0.0175 

0.0020 

0.0003 

0.0059 

54 

-0.0414 

-0.0234 

-0.0246 

0.0011 

0.0005 

0  .0080 

55 

0.0093 

-  0.0008 

-0.0321 

0.0002 

0  .0000 

0.0026 

56 

-0.0248 

-0.0234 

-0.1041 

0.0015 

0.0008 

0.0060 

57 

0.0978 

-0.0309 

-0.0850 

0.0030 

0.0009 

0.0065 

58 

0.0073 

-0  .0000 

-0  .0464 

0.0001 

0.0000 

0.0021 

59 

0.0567 

-0.0062 

0  .0215 

0.0006 

0.0001 

0.0041 

60 

-0.0192 

-0.0311 

0.0449 

0.0018 

0.0008 

0.0058 

61 

None  Seen 

No  Repeats 

0.0108 

62 

0.0709 

-0.0555 

0.0019 

-0.0782 

0.0031 

0.0013 

0.0054 

63 

-0.0170 

-0.0309 

-0.1159 

0.0016 

0.0008 

0.0073 

64 

None  Seen 

No  Repeats 

-0.0460 

65 

-0.0080 

-0.0037 

0.0032 
-0  .1562 

0.0002 

0.0001 

0.0093 

42 


Table  4.7  Continued 


Intercept/ 

Slope 

Allowed 

Reper  ted 

Rejecced 

Position 

Change 

Change 

Change 

66 

-0.0012 

No  Repeats 

-0.0388 

0.0000 

0.0066 

67 

-0  .0048 

-0.0020 

-0.0175 

0.0009 

0.0003 

0.0007 

68 

-0.0118 

-0.0123 

0.0030 

0.0007 

0.0002 

0  .0034 

69 

0.0984 

-0.1056 

-0.0030 

0.0071 

0  .0043 

0  .0023 

70 

-0  .0008 

0  .0036 

-0.0262 

0.0004 

0.0001 

0.0014 

71 

0.1193 

0.0130 

0.1506 

0.0019 

0.0008 

0.0057 

72 

-0.0215 

-0.0206 

0.0155 

0.0014 

0.0008 

0  .0071 

73 

-0.0243 

-0.0087 

0.0186 

0.0009 

0.0002 

0.0030 

74 

0.0513 

-0.0567 

0.0511 

0.0030 

0.0012 

0  .0055 

75 

-0.0412 

-0.0480 

0.0079 

0.0024 

0.0009 

0  .0045 

76 

-0  .0247 

-0.0045 

0.0776 

0.0015 

0.0004 

0.0030 

77 

0.0231 

-0.0635 

0.0330 

0.0030 

0.0015 

0.0057 

78 

-0.0161 

-0.0832 

-0.1682 

0.0040 

0.0022 

0  .0040 

79 

-0.0025 

-0.0012 

0.0136 

0.0000 

0.0000 

0.0014 

80 

-0.0425 

-0.0814 

0.0052 

0.0037 

0.0019 

0.0043 

81 

-0.0317 

-0.0180 

0.0401 

0.0017 

0.0007 

0.0066 

82 

-0.0375 

-0.0117 

-0.0416 

0.0011 

0.0002 

0.0008 

83 

0.0044 

-0.0135 

-0.0184 

0.0015 

0.0002 

0.0065 

84 

-0  .0378 

-0.0258 

0.0981 

0.0010 

0.0004 

0.0071 

85 

-0.0293 

-0.0279 

-0.0482 

0.0019 

0.0009 

0.0023 

86 

0.1172 

-0.0174 

-0.1215 

0.0022 

0.0010 

0.0068 

87 

-0.0129 

-0.0441 

-0.0985 

0.0025 

0.0011 

0.0068 

88 

None  Seen 

No  Repeats 

0.1238 

0.0023 

89 

0.0008 

-0.0185 

-0.0172 

0  .0020 

0.0004 

0.0078 

Table  4.7  Continued 


Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

90 

0.0103 

-0.0301 

0.1865 

0  .0022 

0  .0010 

0  .0045 

91 

-0.0019 

-0.0006 

0.0537 

0.0000 

0.0000 

0.0011 

92 

-0  .0115 

-0.0161 

0  .0829 

0.0009 

0.0004 

0  .0036 

93 

-0.0110 

-0.0129 

0.1076 

0.0002 

0  .0002 

0  .0061 

94 

None  Seen 

No  Repeats 

-0.0673 

95 

-0.0244 

-0  .0132 

0  .0021 

0  .0683 

0.0017 

0.0003 

0  .0054 

96 

-0.0507 

-0.0299 

-0  .0461 

0.0011 

0  .0005 

0.0063 

97 

-0.0031 

-0.0012 

-0.0704 

0.0001 

0  .0000 

0.0017 

98 

-0  .0032 

-0.0161 

-0  .0433 

0.0013 

0.0005 

0.0045 

99 

-0.0260 

-0.0134 

-0.0892 

0.0010 

0.0005 

0  .0061 

100 

None  Seen 

No  Repeats 

-0.0051 

101 

0.0060 

-0.0186 

0.0002 

-0.1118 

0.0020 

0.0006 

0.0046 

102 

0.0193 

-0.0160 

0.0385 

0.0014 

0.0004 

0.0078 

103 

None  Seen 

No  Repeats 

0.0145 

104 

0.0115 

0.0066 

0.0048 
-0 .0791 

0  .0002 

0  .0001 

0  .0069 

105 

0  .0018 

-0.0000 

0.0456 

0.0001 

0.0000 

0.0050 

106 

-0  .0428 

-0  .0257 

-0  .0138 

0.0015 

0.0006 

0.0004 

107 

0  .0480 

-0.0164 

0  .0581 

0.0020 

0.0008 

0.0041 

108 

0.1660 

-0.0615 

0  .0822 

0.0043 

0.0018 

0.0056 

109 

None  Seen 

No  Repeats 

-0.0484 

110 

-0.0518 

-0.0382 

0.0039 

0  .0015 

0.0015 

0  .0008 

0.0061 

111 

0.0393 

-0.0025 

-0.0255 

0.0013 

0  .0000 

0  .0046 

112 

0.0048 

No  Repeats 

-0.0139 

113 

0.0001 

0 .0685 

0.0114 

0  .0015 
0.0073 

0.0017 

0.0003 

0  .0044 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

114 

0.0153 

-0  .00562 

0.0671 

0.0004 

0.0003 

0.0104 

115 

None  Seen 

No  Repeats 

0.1773 

116 

0.1958 

0.0077 

0.0018 

0.0145 

0  .0033 

0.0015 

0.0044 

117 

-0.0411 

-0  .0272 

0.0012 

0.0012 

0.0007 

0.0077 

118 

0.0024 

No  Repeats 

-0.0328 

119 

0.0000 

0.0034 

-0.0031 

0.0027 

0.0030 

0.0006 

0.0001 

0  .0095 

120 

0  .0483 

-0.0113 

0  .0680 

0.0013 

0  .0007 

0  .0052 

121 

-0.0025 

-0.0012 

-0.0068 

0.0000 

0.0000 

0.0031 

122 

-0.0298 

-0.0171 

-0.0125 

0.0008 

0.0004 

0.0036 

123 

0.0313 

-0.0221 

-0.0510 

0.0019 

0.0009 

0  .0057 

124 

None  Seen 

No  Repeats 

0.0522 

125 

None  Seen 

No  Repeats 

0.0025 

0  .0283 

126 

0  .0588 

-0.0720 

0.0037 

0.1128 

0.0041 

0.0020 

0.0038 

127 

0.0136 

-0.0017 

0.0274 

0.0003 

0.0000 

0.0023 

128 

0.0323 

-0.0015 

-0.0116 

0.0016 

0  .0002 

0.0063 

129 

-0.0127 

-0.0066 

-0.0463 

0.0006 

0.0002 

0  .0062 

130 

None  Seen 

No  Repeats 

-0.0294 

131 

-0.0121 

No  Repeats 

0.0043 

0.0221 

132 

0.0012 

-0.0031 

No  Repeats 

0.0058 

-0.0117 

133 

0.0001 

None  Seen 

No  Repeats 

0.0039 

0.0003 

134 

None  Seen 

No  Repeats 

0.0040 

0.0621 

135 

None  Seen 

No  Repeats 

0.0065 

0  .0969 

136 

-0  .0979 

-0.0439 

0.0058 

0.0279 

0.0030 

0.0012 

0.0010 

137 

0  .0893 

0.0101 

-0  .1334 

0.0015 

0.0005 

0  .0069 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

13  8 

-0.0012 

No  Repeats 

-0.0873 

13  9 

0.0000 

None  Seen 

No  Repeats 

0.0092 

-0.0189 

140 

-0.0158 

-0.0131 

0.0033 
-0  .1156 

0.0009 

0.0005 

0  .0086 

141 

0.0316 

-0.0031 

-0  .0570 

0.0007 

0.0001 

0.0069 

142 

-0.0161 

No  Repeats 

0  .0626 

143 

0.0004 

-0.0155 

-0.0006 

0.0027 

0.1275 

0.0004 

0.0000 

0.0058 

144 

-0.0050 

No  Repeats 

0  .0674 

145 

0.0003 

None  Seen 

No  Repeats 

0.0064 

-0.0049 

146 

0.0024 

-0.0000 

0.0001 

-0.1568 

0.0000 

0.0000 

0  .0083 

147 

0.0146 

-0.0000 

-0.0180 

0.0001 

0.0000 

0.0069 

148 

None  Seen 

No  Repeats 

0.0566 

149 

a. 0000 

No  Repeats 

0.0015 

-0.1345 

150 

0.0000 

None  Seen 

No  Repeats 

0.0084 

-0.1076 

151 

None  Seen 

No  Repeats 

0.0064 

-0.0176 

152 

0.0123 

-0.0008 

0.0010 

0.0670 

0.0003 

0.0001 

0.0051 

153 

0.0403 

-0.0188 

0  .0262 

0.0013 

0.0003 

0.0058 

154 

-0.1360 

-0.1065 

0  .0043 

0.0042 

0.0020 

0.0010 

155 

0.0227 

-0.0039 

-0.0897 

0.0004 

0.0001 

0.0080 

156 

0.0131 

-0.0023 

-0.0663 

0.0003 

0.0001 

0.0067 

157 

-0  .0483 

-0.0279 

0.1042 

0.0018 

0.0008 

0.0005 

158 

-0.0027 

-0.0017 

-0.0748 

0.0001 

0.0000 

0  .0079 

159 

0.0774 

-0.0260 

0.1247 

0.0025 

0.0016 

0.0049 

160 

0.0030 

-0.0000 

0.0218 

0.0001 

0.0000 

0.0014 

161 

0.0267 

No  Repeats 

0.0038 

0.0001 

0.0064 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

162 

None  Seen 

No  Repeats 

-0.1115 

163 

None  Seen 

No  Repeats 

0.0080 

0.1657 

164 

None  Seen 

No  Repeats 

0.0081 

0  .0346 

165 

None  Seen 

No  Repeats 

0.0050 

0.0589 

166 

-0.0004 

No  Repeats 

0  .0082 
-0.0295 

167 

0.0000 

0.0034 

-0.0623 

0.0023 

0.0002 

0.0029 

0.0017 

0.0042 

168 

0.0102 

-0.0217 

-0.0905 

0.0009 

0.0003 

0.0072 

169 

None  Seen 

No  Repeats 

-0.0571 

170 

None  Seen 

No  Repeats 

0  .0047 
0.0160 

171 

None  Seen 

No  Repeats 

0.0080 

-0.0512 

172 

None  Seen 

No  Repeats 

0.0091 

-0.0643 

173 

-0.0638 

-0.0636 

0.0020 
-0  .0292 

0.0022 

0.0016 

0.0075 

174 

0.0895 

No  Repeats 

-0.0263 

175 

0.0010 

0.1249 

No  Repeats 

0.0072 

-0.1250 

176 

0.0008 

0  .0502 

-0.0190 

0  .0042 
-0.1460 

0.0009 

0.0003 

0.0076 

177 

None  Seen 

No  Repeats 

-0.0490 

178 

-0.0098 

-0.0061 

0.0078 

0.0098 

0.0001 

0.0001 

0.0025 

179 

-0.0241 

-0.0037 

-0.0352 

0.0006 

0.0001 

0.0061 

180 

0.0324 

No  Repeats 

0.0907 

181 

0.0007 

None  Seen 

No  Repeats 

0.0078 
-0  .0440 

182 

0.0275 

-0.0019 

0.0028 

0.0128 

0.0010 

0.0000 

0.0071 

183 

-0.0119 

-0.0050 

0.0982 

0.0014 

0.0005 

0.0061 

184 

-0  .0883 

-0.0407 

0.0233 

0.0027 

0.0009 

0.0014 

185 

0.0773 

-0.0019 

-0  .1211 

0.0012 

0  .0000 

0.0069 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

186 

-0.0589 

-0.0611 

-0.0389 

0.0023 

0.0015 

0  .0051 

187 

None  Seen 

No  Repeats 

0.0635 

188 

-0.0569 

-0.0272 

0.0027 

0.1970 

0.0015 

0.0007 

0.0068 

189 

0.0085 

-0.0000 

-0.1098 

0.0012 

0.0000 

0.0086 

190 

None  Seen 

No  Repeats 

-0.0767 

191 

0.0103 

-0  .0487 

0.0017 

0.0366 

0.0021 

0.0010 

0.0079 

192 

0.1047 

-0.0025 

-0.0567 

0.0010 

0.0001 

0.0076 

193 

None  Seen 

No  Repeats 

-0.0023 

194 

0.0596 

0.0001 

0.0042 
-0  .0832 

0.0020 

0.0010 

0.0064 

195 

0.1308 

0.0005 

0  .0830 

0.0028 

0.0015 

0.0043 

196 

None  Seen 

No  Repeats 

-0.0677 

197 

0.0279 

No  Repeats 

0.0031 

-0.0426 

198 

0.0007 

0.0033 

No  Repeats 

0.0062 

0.0568 

199 

0.0000 

None  Seen 

No  Repeats 

0.0073 

0  .0939 

200 

0.0695 

-0.0160 

0.0016 

0.1742 

0.0010 

0.0003 

0.0051 

201 

-0.0487 

-0.0597 

-0.0813 

0.0027 

0.0012 

0.0065 

202 

None  Seen 

No  Repeats 

0.0721 

203 

None  Seen 

No  Repeats 

0  .0021 
0.1042 

204 

None  Seen 

No  Repeats 

0.0060 

0.0227 

205 

None  Seen 

No  Repeats 

0.0063 

-0.0614 

206 

-0.0604 

-0.0418 

0.0012 

-0.1538 

0.0012 

0  .0007 

0.0056 

207 

0.0591 

0  .0038 

-0  .0358 

0.0010 

0.0001 

0.0057 

208 

-0.0050 

-0.0000 

0  .1785 

0  .0001 

0.0000 

0.0016 

209 

-0  .0353 

-0.0025 

-0  .1131 

0.0013 

0.0000 

0.0065 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

210 

-0.0037 

-0.0000 

0.1737 

0.0001 

0.0000 

0.0105 

211 

0.0587 

No  Repeats 

-0.0441 

212 

0.0008 

0.0036 

-0.0033 

0.0026 

0.1166 

0.0020 

0.0005 

0.0039 

213 

0.0227 

-0.0055 

0.0342 

0.0005 

0.0001 

0.0054 

214 

None  Seen 

No  Repeats 

-0.0617 

215 

-0.0044 

-0.0041 

0.0029 

-0.0722 

0  .0003 

0  .0001 

0.0079 

216 

0.0032 

0.0006 

-0.0016 

0.0001 

0.0000 

0.0092 

217 

0.0673 

No  Repeats 

-0.0636 

218 

0.0004 

0.0352 

-0.0260 

0.0024 

0.0754 

0.0023 

0.0011 

0.0050 

219 

None  Seen 

No  Repeats 

-0  .0539 

220 

-0.0006 

No  Repeats 

0.0075 

0.0300 

221 

0.0000 

0.0853 

-0.0021 

0.0010 

-0.0740 

0.0013 

0.0001 

0.0064 

222 

None  Seen 

No  Repeats 

0.0138 

223 

-0.0218 

-0.0043 

0.0089 

0.0336 

0.0004 

0.0001 

0.0005 

224 

-0.0173 

-0.0012 

0.0278 

0.0010 

0.0000 

0.0052 

225 

-0.0203 

-0.0043 

0.1336 

0.0004 

0.0001 

0.0061 

226 

-0.0385 

-0.0183 

-0.0539 

0.0014 

0.0004 

0.0019 

227 

0.0575 

-0.0188 

-0.0251 

0  .0020 

0.0007 

0.0051 

228 

0.0086 

-0.0319 

0.0011 

0.0017 

C  .0009 

0.0053 

229 

0.1118 

No  Repeats 

-0  .0748 

230 

0.0008 

-0.0426 

No  Repeats 

0.0024 

-0.0756 

231 

0.0021 

-0.0610 

-0.0240 

0.0058 

0.0754 

0.0010 

0.0004 

0.0053 

232 

None  Seen 

No  Repeats 

-0 .0244 

233 

-0.0052 

-0  .0272 

0.0005 

0.0196 

0.0012 

0.0007 

0.0054 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

234 

0.0018 

-0.0000 

0.1055 

0.0001 

0.0000 

0  .0074 

235 

0.1115 

No  Repeats 

-0.0544 

236 

0.0008 

0.0212 

-0  .0388 

0.0020 

-0.0008 

0.0021 

0.0008 

0.0057 

237 

-0.0128 

-0.0138 

-0.0029 

0.0015 

0.0006 

0.0068 

23  8 

None  Seen 

No  Repeats 

-0.0862 

239 

0.0503 

No  Repeats 

0.0036 

0.0482 

240 

0.0003 

0.0567 

No  Repeats 

0.0059 

-0.0165 

241 

0.0009 

None  Seen 

No  Repeats 

0.0089 

0.0661 

242 

-0.0192 

-0.0092 

0.0012 

-0.0165 

0.0005 

0.0002 

0.0063 

243 

-0.0069 

-0.0366 

0.0291 

0.0028 

0.0015 

0.0060 

244 

0.1090 

-0.0106 

-0.0427 

0.0012 

0.0002 

0  .0043 

245 

-0.0006 

-0.0031 

0.0361 

0.0001 

0.0001 

0.0077 

246 

0.0151 

-0.0364 

0.0180 

0.0021 

0.0011 

0.0070 

2  47 

None  Seen 

No  Repeats 

-0.0540 

248 

0.0278 

-0.0532 

0.0024 
-0  .0539 

0.0023 

0.0014 

0.0061 

249 

0.1042 

-0.0125 

0.0667 

0.0016 

0.0009 

0.0068 

250 

0.0433 

-0.0056 

0.0712 

0.0009 

0.0002 

0.0017 

251 

-0.0128 

-0.0062 

0.0619 

0.0007 

0.0003 

0  .0085 

252 

0.0044 

-0.1361 

0.1065 

0.0066 

0.0043 

0  .0047 

253 

None  Seen 

No  Repeats 

-0.0007 

254 

0.0263 

No  Repeats 

0.0021 

0.0280 

255 

0.0008 

-0.0025 

No  Repeats 

0  .0064 

0  .0929 

256 

0.0000 

-0.0590 

-0.0272 

0.0055 

0.0659 

0.0024 

0.0007 

0.0014 

257 

0.0016 

-0.0710 

-0  .0398 

0  .0033 

0.0021 

0.0078 
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Intercept/ 

^  Slope 

Allowed 

Repeated 

Rejected 
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Change 

Change 

Change 

258 

0.0244 

-0.0030 

0.1344 

P 

0.0003 

0.0002 

0.0077 

259 

0.0031 

-0.0028 

-0.0296 

i'. 

0.0012 

0.0004 

0.0011 

-  260 

0.0635 

-0.0269 

-0.0279 

1 

0.0010 

0.0007 

0.0050 

261 

0.0181 

-0.0016 

0.0407 

0.0014 

0  .0008 

0  .0055 

262 

0.0100 

No  Repeats 

0.0595 

0.0002 

0.0014 

%  263 

-0.0146 

-0.0044 

0.0470 

I 

0.0007 

0.0002 

0  .0049 

v  264 

-0.0138 

-0.0037 

0.1057 

0.0006 

0.0001 

0  .0053 

265 

None  Seen 

No  Repeats 

-0  .0393 

0.0018 

266 

-0.0394 

-0.0296 

0.0148 

M 

0.0016 

0.0007 

0.0074 

267 

0.0940 

No  Repeats 

-0.0816 

0.0011 

0.0060 

268 

-0.0196 

-0.0075 

0.0810 

K* 

0.0008 

0.0003 

0.0023 

■  269 

0.0746 

-0.0157 

0.0216 

1 

0.0019 

0.0006 

0.0083 

270 

-0.0587 

-0.0326 

0.0085 

0.0024 

0.0009 

0.0037 

>  271 

0.0018 

-0.0000 

0.0003 

' 

0.0001 

0.0000 

0.0016 

|  272 

-0.0365 

-0.0173 

-0  .0072 

■ 

0.0010 

0.0004 

0.0061 

273 

-0.0032 

0.0018 

0.1048 

0.0002 

0.0000 

0.0067 

274 

-0.0055 

-0.0025 

0.1639 

0.0001 

0.0000 

0.0062 

Z  275 

-0.0684 

-0.0382 

-0.0214 

R 

0.0016 

0.0008 

0.0067 

276 

0.0464 

-0.0467 

-0.17  43 

0.0028 

0.0009 

0.0064 

277 

-0.0031 

No  Repeats 

0.0102 

>  ■, 
i* 

0.0001 

0.0015 

kj  *  278 

0.0307 

-0.0049 

-0.1415 

P* 

0.0009 

0.0001 

0.0075 

279 

0.0174 

-0.0031 

0  .0851 

■  ' 

0.0007 

0.0001 

0.0064 

‘  280 

0.0052 

-0.0032 

-0.0254 

0.0004 

0.0003 

0.0043 

jj  281 

0  .0054 

-0.0013 

-0.0052 

0.0002 

0.0001 

0.0073 

J 
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Intercept/ 

Slope 

Allowed 

Position 

Change 

282 

0.0263 

0.0003 

283 

-0.0917 

0.0027 

284 

-0.00  96 
0.0004 

285 

0.0521 

0.0011 

286 

-0  .0333 
0.0006 

287 

0.0239 

0.0006 

288 

0.0012 

0.0006 

289 

0.0029 

0.0002 

290 

0.0871 

0.0012 

291 

0.0277 

0.0021 

292 

None  Seen 

293 

-0.0031 

0.0001 

294 

None  Seen 

295 

-0.0694 

0.0017 

296 

0.0285 

0  .0009 

297 

0.0152 

0.0018 

298 

-0.0752 

0.0018 

299 

0.0287 

0.0038 

100 

-0.0024 

0.0000 

$01 

0.0152 

0.0001 

$02 

-0  .0028 
0.0026 

103 

0.0338 

0.0002 

$04 

None  Seen 

$05 

-0.0605 

0.0028 

Repeated 

Rejected 

Change 

Change 

-0.0001 

0.0268 

0.0001 

0  .0086 

-0.0519 

-0.0296 

0.0013 

0  .0026 

-0.0111 

-0.0727 

0.0002 

0  .0078 

-0.0019 

0.0460 

0.0000 

0.0062 

-0.0142 

-0.0176 

0.0003 

0.0008 

No  Repeats 

0.0028 

0.0050 

-0.0002 

0.0568 

0.0003 

0.0068 

0.0349 

0.0015 

-0.0063 

-0  .0435 

0.0006 

0.0072 

-0.0204 

-0.1257 

0.0005 

0.0064 

No  Repeats 

None  Rejected 

*o  Repeats 

-0.0365 

0.0062 

No  Repeats 

-0.0480 

0.0077 

-0.0306 

-0.0664 

0.0006 

0.0025 

-0.0032 

0.0524 

0.0003 

0.0059 

-0.0111 

-0.0036 

0.0003 

0.0051 

-0.0253 

0.0278 

0.0005 

0.0014 

-0  .0833 

0.0422 

0.0028 

0.0071 

No  Repeats 

-0.0675 

0.0097 

0.0073 

-0  .0368 

0.0000 

0.0036 

-0  .0418 

-0  .0783 

0.0008 

0.0047 

No  Repeats 

0.1787 

0 .0068 

No  Repeats 

0.1483 

0.0037 

-0.0408 

0.2168 

0.0012 

0  .0051 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

306 

0.0735 

-0.0012 

-0.0665 

0.0011 

0.0000 

0.0072 

307 

None  Seen 

No  Repeats 

-0.0701 

308 

None  Seen 

No  Repeats 

0.0033 

0.0412 

309 

None  Seen 

No  Repeats 

0.0045 

-0.0844 

310 

None  Seen 

No  Repeats 

0.0050 

-0.0115 

311 

0.0236 

-0.0007 

0.0034 

0.1183 

0.0005 

0.0001 

0.0076 

312 

0.0518 

-0.0058 

0  .0973 

0.0009 

0  .0001 

0.0059 

313 

None  Seen 

No  Repeats 

-0.0492 

314 

-0.0051 

No  Repeats 

0.0021 

-0.1122 

315 

0.0006 

-0.0265 

-0.0129 

0.0056 

-0.0200 

0.0005 

0.0002 

0.0049 

316 

-0.0051 

No  Repeats 

-0  .0074 

317 

0.0002 

-0.0140 

-0.0049 

0.0017 

0.0091 

0.0009 

0.0001 

0.0044 

318 

-0.0087 

-0.0019 

-0.0980 

0.0003 

0.0000 

0.0046 

319 

None  Seen 

No  Repeats 

-0.0302 

320 

-0.0631 

-0.0482 

0.0014 

0.1247 

0.0020 

0.0012 

0  .0090 

321 

0.0503 

No  Repeats 

0.0285 

322 

0.0016 

-0.0061 

-0.0025 

0.0060 

0.0015 

0.0001 

0.0000 

0.0053 

323 

-0.0392 

-0.0222 

0  .0051 

0.0012 

0.0004 

0  .0093 

324 

0  .0892 

-0.0199 

-0.0600 

0.0027 

0.0006 

0.0071 

325 

None  Seen 

No  Repeats 

-0.0241 

326 

0.0308 

0.0024 

0.0013 

0.0139 

0.0000 

-0.0000 

0.0061 

327 

0.0440 

0.0018 

-0.0036 

0.0013 

0.0001 

0.0048 

328 

0.0158 

0.0037 

0  .1169 

0.0001 

0.0000 

0.0035 

329 

0.0060 

-0.0006 

0.0770 

0.0008 

0.0000 

0.0058 

Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

330 

-0.0332 

-0.0190 

-0.0158 

0.0006 

0.0003 

0.0044 

331 

-0  .0004 

No  Repeats 

-0.0160 

332 

0.0000 
-0.00  47 

-0.0008 

0.0002 

0.0086 

0.0001 

0.0000 

0  .0051 

333 

-0.0039 

-0.0010 

-0.1022 

0.0001 

0.0000 

0.0088 

334 

None  Seen 

No  Repeats 

0.0168 

335 

-0.0013 

-0.0025 

0.0026 
-0  .0469 

0.0001 

0.0001 

0  .0054 

336 

0.0610 

0.0024 

-0.0624 

0.0008 

0.0000 

0.0070 

337 

None  Seen 

No  Repeats 

-0.0726 

338 

0.0040 

No  Repeats 

0.0015 

-0.1484 

339 

0.0000 

0.0285 

No  Repeats 

0.0048 

-0.0120 

3  40 

0.0012 

-0.0062 

-0.0025 

0.0065 

0.0007 

0.0001 

0.0001 

0.0007 

3  41 

-0 .0320 

-0.0135 

0.0629 

0.0007 

0.0002 

0.0043 

3  42 

0.1032 

No  Repeats 

-0  .0902 

3  43 

0.0011 

None  Seen 

No  Repeats 

0.0069 
-0.03  82 

344 

0.0295 

0.0073 

0.0012 

0.0270 

0.0004 

0.0000 

0.0058 

345 

0.0053 

0.0012 

-0.0122 

0.0000 

0.0000 

0.0065 

3  46 

None  Seen 

No  Repeats 

-0.0716 

3  47 

None  Seen 

No  Repeats 

0.0022 

-0.1222 

348 

0.0070 

No  Repeats 

0.0070 
-0  .0690 

349 

0.0000 

-0.0102 

-0.0023 

0.0075 

0.1705 

0.0001 

0.0000 

0.0035 

350 

0.0407 

No  Repeats 

0  .0660 

351 

0.0011 

-0.0078 

-0.0025 

0.0050 

0.0692 

0.0001 

0.0000 

0.0076 

352 

None  Seen 

No  Repeats 

-0.0065 

353 

0.0207 

-0.0037 

0.0015 

-0.0244 

0.0005 

0.0000 

0  .0062 

Table  4.7  Continued 


Intercept/ 

Slope 

Position 

354 

355 

356 
3  57 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 


Allowed 
Change 
0.0140 
0.0010 
None  Seen 

-0.0327 
0.0008 
0.0201 
0.0005 
None  Seen 

0.0027 
0.0005 
None  Seen 

None  Seen 

0.0669 
0.0014 
-0.0593 
0.0015 
-0.03  93 
0.0006 
-0  .06  81 
0.0014 
0.0102 
0.0006 
-0.0572 
0.0020 
0.0600 
0.0024 
0.0006 
0.0000 
None  Seen 

0.0256 
0.0017 
0.1330 
0.0020 
-0.0046 
0.0001 
0.0559 
0.0019 
0.0299 
0.0004 
0.0000 
0.0000 
-0  .0047 
0.0008 


Repeated 
Change 
-0.0074 
0.0001 
No  Repeats 

-0.0210 
0.0004 
No  Repeats 

No  Repeats 

-0.0012 
0.0000 
No  Repeats 

No  Repeats 

-0.0049 
0  .0001 
-0.0334 
0  .0009 
-0.0141 
0.0002 
-0.0412 
0.0007 
No  Repeats 

-0.0222 
0.0006 
-0.0098 
0.0005 
No  Repeats 

No  Repeats 

-0.0317 
0.0008 
-0.0068 
0.0004 
-0.0045 
0.0001 
-0.0218 
0.0007 
-0.0010 
0.0000 
No  Repeats 

0.0012 

0.0001 


Rejected 
Change 
-0.0329 
0  .0078 
-0.0504 
0.0018 
-0.0633 
0.0077 
-0  .0385 
0.0065 
-0.0317 
0.0041 
0  .0626 
0.0063 
-0.0661 
0.0067 
0  .0702 
0.0017 
0.0364 
0.0080 
.  -0.0735 
0.0061 
0.0122 
0.0022 
0.1755 
0.0047 
-0  .0432 
0.0081 
0.0192 
0.0018 
-0  .0194 
0.0040 
-0.0320 
0.0082 
-0  .0357 
0.0012 
0.0276 
0.0059 
-0.0480 
0.0068 
0.0115 
0.0016 
-0.1633 
0.0047 
0.0766 
0.0065 
-0.0319 
0.0019 
-0  .0238 
0.0046 


Table  4.7  Continued 


Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

378 

0.0054 

0.0048 

0  .0630 

379 

0.0002 

0.0001 

0  .0074 

0.0446 

No  Repeats 

-0.0548 

3  80 

0.0008 

0.0073 

-0.0426 

0  .0020 
-0.0435 

3  81 

0.0019 

0.0010 

0.0066 

-0.0236 

-0.0265 

0  .0968 

3  82 

0.0012 

0.0006 

0.0063 

-0.0037 

No  Repeats 

-0.0182 

383 

0.0001 

0.0069 

-0.0012 

0.0030 

0.0723 

3  84 

0.0001 

0  .0000 

0  .0079 

0.0839 

-0  .0783 

0.0225 

385 

0.0034 

0.0017 

0.0053 

-0.0066 

-0.0012 

-0  .0946 

3  86 

0.0001 

0.0000 

0.0019 

0.0582 

-0.0149 

-0.0679 

0.0024 

0.0004 

0.0045 

3  87 

-0.0417 

-0.0348 

0.0773 

388 

0.0014 

0.0008 

0.0077 

None  Seen 

No  Repeats 

0.0324 

389 

-0.0061 

No  Repeats 

0  .0050 
-0.1540 

3  90 

0.0001 

None  Seen 

No  Repeats 

0.0078 

0.0833 

3  91 

None  Seen 

No  Repeats 

0.0075 

-0.1154 

3  92 

-0.0245 

-0.0246 

0.0041 

0.1098 

3  93 

0.0016 

0.0008 

0  .0077 

0.0738 

-0.0006 

0.1109 

394 

0.0013 

0.0000 

0.0068 

-0.0166 

-0.0151 

0.0813 

395 

0.0004 

0  .0002 

0  .0068 

-0.0133 

-0.0078 

-0.1425 

3  96 

0.0010 

0.0001 

0.0086 

-0.0448 

-0.0279 

-0.0011 

3  97 

0.0011 

0.0005 

0.0061 

None  Seen 

No  Repeats 

-0.0155 

398 

0.0556 

-0.0006 

0.0018 

-0.0729 

399 

0.0009 

0.0001 

0.0058 

0.0006 

-0.0000 

-0  .1685 

400 

0.0001 

0.0000 

0.0069 

None  Seen 

No  Repeats 

-0.0020 

J01 

None  Seen 

No  Repeats 

0.0003 
-0 .0000 

0.0056 

56 


Table  4.7  Continued 


Intercept/ 

Slope 

Allowed 

Repeated 

Rejected 

Position 

Change 

Change 

Change 

402 

None  Seen 

No  Repeats 

-0.2218 

403 

-0.0093 

-0.0043 

0.0102 

-0.0291 

0.0002 

0.0001 

0.0007 

404 

-0.0219 

-0.0068 

0.0720 

0.0007 

0.0001 

0.0065 

405 

0.0561 

-0.0012 

0.0662 

0.0011 

0.0000 

0.0044 

406 

None  Seen 

No  Repeats 

-0.0081 

407 

0.0619 

0.0036 

0.0005 

-0.1284 

0.0007 

0.0000 

0.0058 

408 

-0.0309 

-0.0124 

0.0727 

0.0010 

0.0004 

0.0071 

409 

-0.0500 

-0.0117 

-0.0215 

0.0011 

0.0002 

0.0006 

410 

-0.0215 

-0.0735 

-0.0077 

0.0011 

0.0002 

0.0006 

411 

-0.0946 

-0.0704 

-0.1072 

0.0030 

0.0017 

0.0052 

412 

-0.0364 

-0.0148 

0.0179 

0.0008 

0.0004 

0.0088 

413 

-0.0569 

-0.0250 

-0.0182 

0.0021 

0.0008 

0.0006 

414 

0.0473 

-0.0557 

-0.0364 

0.0033 

0.0017 

0.0077 

415 

-0.0221 

-0.0135 

-0.0686 

0.0003 

0.0002 

0.0034 

416 

0.0290 

0.0024 

-0.0758 

0.0000 

0.0000 

0.0073 

417 

0.0952 

-0.0000 

0.0352 

0.0013 

0.0001 

0.0058 

418 

-0.0184 

-0.0074 

-0.0277 

0.0003 

0.0001 

0.0022 

419 

0.0102 

-0.0222 

-0.0072 

0.0018 

0.0005 

0.0059 

420 

-0.0181 

-0.0307 

-0.1104 

0.0015 

0.0005 

0.0078 

421 

None  Seen 

No  Repeats 

-0  .0450 

422 

0.0582 

-0  .0037 

0.0029 

0.0813 

0.0007 

0.0001 

0.0069 

423 

0.0048 

No  Repeats 

0  .0086 

424 

0.0001 

None  Seen 

No  Repeats 

0.0087 

-0.0207 

425 

-0.0118 

-0.0089 

0.0010 

-0.0347 

0.0004 

0.0002 

0.0076 

57 
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Intercept/ 

Slope 

Allowed 

Repeated 

Rej  ected 

Position 

Change 

Change 

Change 

426 

None  Seen 

No  Repeats 

-0.0663 

427 

0.0180 

No  Repeats 

0.0067 

-0.1183 

428 

0.0001 

-0.0213 

-0.0123 

0.0039 
0.07  93 

0.0003 

0.0002 

0.0072 

429 

0.0042 

No  Repeats 

-0  .0209 

430 

0.0001 

None  Seen 

No  Repeats 

0.0055 
-0  .0503 

431 

-0.0312 

-0.0472 

0.0034 

0  .0322 

0.0024 

0  .0012 

0 .0056 

432 

0.1037 

No  Repeats 

0  .1043 

433 

0.0006 

None  Seen 

No  Repeats 

0.0073 

0  .0396 

434 

0.0256 

-0.0154 

0.0040 

-0.0015 

0.0014 

0.0003 

0.0063 

435 

0.0067 

-0.0166 

-0  .0364 

0.0013 

0.0002 

0  .0044 

436 

None  Seen 

No  Repeats 

0.0054 

437 

-0  .0265 

-0.0315 

0.0027 

-0.1455 

0.0017 

0.0007 

0.0071 

438 

-0.0147 

-0.0098 

0.0123 

0.0002 

0.0001 

0  .0081 

Table  4.8  contains  the  regression  analysis  by  position 
for  the  transition  of  amino  acid  residues  as  a  function  of 
iteration. 


Table  4.8 

Amino  Acid  Transition  Regression 


Intercept/ 

Slope  Accepted  Reverse 


Amino 

Ancestral 

Human 

Rabbit 

Acid 

Change 

Change 

Change 

Change 

1 

0.0425 

0.0425 

0.0425 

-0.0051 

0.0009 

0.0009 

0.0009 

0  .0031 

2 

0.0526 

0 .0526 

0  .0526 

-0  .0066 

0.0012 

0.0012 

0  .0012 

o 

• 

o 

3 

0 .6671 

0  .3180 

0.6671 

-0  .0474 

0.0003 

0.0004 

0.0003 

0.0012 

Table  4.8  Continued 


Intercept/ 

Slope 

Accepted 

Reverse 

Amino 

Ancestral 

Human 

Rabbit 

Acid 

Change 

Change 

Change 

Change 

4 

0  .3793 

0.3793 

0  .6724 

0  .0470 

0.0008 

0.0008 

0.0000 

0.0012 

5 

0.7315 

0.7919 

0  .7503 

0.1607 

0.0005 

0.0004 

0.0004 

0.0023 

6 

0.1434 

0.1434 

0  .1434 

0  .1450 

0.0012 

0.0012 

0.0012 

0.0017 

7 

0.1525 

0.1525 

0  .1525 

0.1433 

0.0010 

0.0010 

0.0010 

0.0039 

8 

0  .0870 

0.0870 

0  .0870 

-0.1708 

0.0012 

0.0012 

0.0012 

0.0072 

9 

0.6741 

0  .3623 

0  .3623 

0  .0008 

0 .0002 

0.0001 

0.0001 

0.0024 

10 

0.1155 

0.1155 

0  .1155 

0.1051 

0.0011 

0.0011 

0.0011 

0.0014 

11 

0.6670 

0.3694 

0  .3694 

0.0341 

0.0001 

0.0005 

0.0005 

0.0013 

12 

0.0000 

0.0000 

0.0000 

Not  Repeated 

0.0000 

0.0000 

0.0000 

13 

0.6688 

0  .3580 

0  .3580 

-0.0440 

0.0004 

0.0002 

0.0002 

0.0029 

14 

0  .0356 

0.0356 

0.0356 

-0.0335 

0.0005 

0.0005 

0.0005 

0.0032 

15 

0.0635 

0.0635 

0.0635 

0.0342 

0.0008 

0.0008 

0.0008 

0.0027 

16 

0.0888 

0.0888 

0.0888 

0  .0411 

0.0008 

0.0008 

0.0008 

0  .0043 

17 

0.0556 

0.0556 

0.0556 

0  .0097 

0.0011 

0.0011 

0.0011 

0.0026 

18 

0.0170 

0.0170 

0.0170 

*'0.0791 

0.0012 

0.0012 

0.0012 

0  .  1038 

19 

0.1349 

0.1349 

0  .1349 

0.0823 

0.0013 

0.0013 

0.0013 

0.0047 

20 

0.0805 

0.0805 

0.0805 

0.0447 

0.0010 

0.0010 

0.0010 

0.0025 

21 

0.4308 

0.7365 

0.4308 

0.0539 

0.0008 

0.0004 

0.0008 

0.0046 

22 

-0.0031 

-0.0031 

-0.0031 

-0  .0093 

0.0001 

0.0001 

0.0001 

0.0002 

23 

0.1762 

0.1762 

0  .1762 

0  .0818 

0.0010 

0.0010 

0.0010 

0.0086 

24 

0.1565 

0.1565 

0.1565 

0.0969 

0.0007 

0.0007 

0.0007 

0  .0037 

25 

0.1565 

0.1565 

0  .1565 

-0.0142 

0.0012 

0.0012 

0.0012 

0.0062 

26 

0  .1970 

0.1970 

0.1970 

-0.0177 

0.0010 

0.0010 

0.0010 

0  .0085 

Table  4.8  Continued 


Intercept/ 

Slope 

Accepted 

Rabbit 

Reverse 

Amino 

Ancestral 

Human 

Change 

Acid 

Change 

Change 

Change 

27 

0.0762 

0.0762 

0.0762 

-0  .0766 

0.0011 

0.0011 

0.0011 

0.0054 

28 

0.0491 

0.0491 

0.0491 

-0.0709 

0.0009 

0.0009 

0  .0009 

0.0036 

29 

0.1765 

0.1765 

0.1765 

0.0750 

0 .0011 

0.0011 

0.0011 

0.0066 

30 

0.7279 

0  .3740 

0  .3740 

0.01110 

0.0004 

0.0008 

0.0008 

0  .0042 

31 

0.0170 

0.0170 

0.0170 

-0  .0243 

0.0003 

0.0003 

0.0003 

0.0013 

32 

0.0190 

0.0190 

0.0190 

-0  .0752 

0.0007 

0.0007 

0.0007 

0.0027 

33 

0.0170 

0.0170 

0 .0170 

-0.0323 

0.0004 

0 .0004 

0.0004 

0.0023 

34 

0.0884 

0  .0884 

0  .0884 

0.0253 

0.0014 

0.0014 

0.0014 

0.0034 

35 

0.0020 

0.0020 

0.0020 

0.0133 

0.0001 

0.0001 

0.0001 

0.0003 

36 

0.2160 

0.2160 

0.2160 

0.1712 

0.0012 

0.0012 

0.0012 

0.0078 

37 

0.0442 

0.0442 

0.0442 

-0.0125 

0.0012 

0.0012 

0.0012 

0.0028 

38 

0.0719 

0.0719 

0.0719 

0.0886 

0.0010 

0.0010 

0.0010 

0.0021 

39 

0.1794 

0.1794 

0.1794 

0.1547 

0.0008 

0  .0008 

0.0008 

0.0044 

40 

0 .0766 

0.0766 

0.0766 

0  .0541 

0.0008 

0.0008 

0.0008 

0.0021 

41 

0 .0504 

0.0504 

0.0504 

-0.0010 

0 .0007 

0.0007 

0.0007 

0.0028 

42 

0.1457 

0.1457 

0.1457 

0  .0588 

0.0011 

0.0011 

0.0011 

0.0041 

43 

0.6605 

9.6146 

9.6146 

0.0331 

0.0004 

0.0007 

0.0007 

0.0025 

44 

-0.0121 

-0.0121 

-0.0121 

-0.0152 

0 .0012 

0.0012 

0.0012 

0.0013 

45 

0.0000 

0  .0000 

0.0000 

Not  Repeated 

46 

0 .0000 
0.0714 

0.0000 

0.0714 

0.0000 

0.0714 

-0.0097 

0  .0009 

0.0009 

0.0009 

0.0045 

47 

0.0372 

0.0372 

0.0372 

0.0159 

0.0006 

0.0006 

0.0006 

0.0016 

48 

-0.0121 

-0.0121 

-0.0121 

-0.0366 

0.0012 

0.0012 

0.0012 

0  .0011 

49 

0.0146 

0.0146 

0  .0146 

0.0170 

0  .0000 

0  .0000 

0.0000 

0.0001 

Table  4.8  Continued 


Intercept/ 


Slope 

Accepted 

Reverse 

Amino 

Ancestral 

Human 

Rabbit 

Acid 

Change 

Change 

Change 

Change 

50 

0  .3333 

0.6667 

0.3333 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

51 

0.3753 

0.7054 

0  .3753 

0.0525 

0.0002 

0.0003 

0.0002 

0.0016 

52 

0.3634 

0.3634 

0.7116 

-0.1003 

0.0002 

0.0002 

0.0004 

0.0048 

53 

0  .0806 

0.0806 

0  .0806 

0.0264 

0  .0004 

0.0004 

0.0004 

0  .0044 

54 

0.0267 

0.0267 

0.0267 

0  .0297 

0.0001 

0.0001 

0.0001 

0.0002 

55 

0.0000 

0.0000 

0  .0000 

Not  Repeated 

0.0000 

0.0000 

0.0000 

56 

0.4010 

0.4010 

0  .6562 

0.0132 

0.0007 

0.0007 

0.0001 

0.0038 

57 

0.0000 

0.0000 

0.0000 

Not  Repeated 

0.0000 

0.0000 

0.0000 

58 

0.0937 

0.0937 

0  .0937 

0.0257 

0.0010 

0.0010 

0.0010 

0.0033 

59 

0.2034 

0.2034 

0.2034 

0.1751 

0.0012 

0.0012 

0.0012 

0.0017 

60 

0.0282 

0.0282 

0.0282 

-0.0015 

0.0008 

0.0008 

0.0008 

0  .0015 

61 

0.0275 

0.0275 

0.0275 

0.0156 

0.0009 

0.0009 

0.0009 

0.0024 

62 

0.0823 

0.0823 

0  .0823 

0.0700 

0.0011 

0.0011 

0.0011 

0.0062 

63 

0.0085 

0.0085 

0  .0085 

-0.0484 

0.0012 

0.0012 

0.0012 

0.0027 

64 

0.1735 

0.1735 

0.1735 

0.1149 

0.0012 

0.0012 

0.0012 

0.0031 

65 

0.1442 

0.1442 

0  .1442 

0.1903 

0.0008 

0.0008 

0.0008 

0.0048 

66 

0.0299 

0  .0299 

0.0299 

0.0312 

0.0006 

0.0006 

0  .0006 

0  .0007 

67 

0.1097 

0  .1097 

0.1097 

0.0208 

0.0012 

0.0012 

0.0012 

0.0036 

68 

0  .0000 

0.0000 

0  .0000 

Not  Repeated 

0  .0000 

0.0000 

0.0000 

69 

0.3827 

0.3827 

0  .6649 

-0.0013 

0.0008 

0.0008 

0.0001 

0.0022 

70 

-0.0290 

-0 .0290 

-0  .0290 

-0.0440 

0.0012 

0.0012 

0.0012 

0.0016 

71 

0.1141 

0.1141 

0  .1141 

0  .0850 

0.0013 

0.0013 

0.0013 

0.0033 

72 

0.6660 

0  .3345 

0  .3345 

-0.0012 

0.0002 

0.0001 

0.0001 

0.0004 
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i 

Intercept/ 

Slope 

Table  4.8 

Accepted 

Continued 

Reverse 

Amino 

Ancestral 

Human 

Rabbit 

Acid 

Change 

Change 

Change 

Change 

73 

0.4418 

0.4418 

0  .6212 

0.1025 

> 

0.0008 

0.0008 

0.0001 

0.0027 

74 

0.0860 

0  .0860 

0.0860 

0  .0847 

i 

0.0012 

0.0012 

0.0012 

0.0013 

75 

-0.0112 

-0.0112 

-0 .0112 

-0.0594 

%  • 

0.0010 

0.0010 

0.0010 

0.0019 

76 

0  .7060 

0  .7009 

0  .7093 

0  .0275 

>  * 

0.0040 

0.0001 

0.0004 

0.0051 

77 

0.1225 

0.1225 

0 .1225 

0  .0383 

l 

0.0013 

0.0013 

0.0013 

0.0039 

78 

0.0343 

0.0343 

0  .0343 

-0.0034 

y 

0.0004 

0.0004 

0.0004 

0.0015 

79 

0.1927 

0.1927 

0.1927 

0.1199 

0.0012 

0.0012 

0.0012 

0.0044 

80 

0.1069 

0.1069 

0  .1069 

0.1069 

a 

0.0012 

0.0012 

0.0012 

0.0012 

81 

0.0460 

0  .0460 

0  .0460 

-0.0261 

0  .0006 

0.0006 

0.0006 

0.0033 

82 

0.1869 

0.1869 

0.1869 

0.1234 

0.0012 

0.0012 

0.0012 

0.0035 

" 

83 

0.1624 

0.1624 

0.1624 

0.13  21 

I 

0.0008 

0.0008 

0.0008 

0.0039 

84 

0.1999 

0.1999 

0.1999 

0.0035 

0  .0009 

0.0009 

0.0009 

0  .0082 

V 

85 

0.0263 

0  .0263 

0.0263 

0  .0239 

0.0008 

0.0008 

0.0008 

0.0008 

■i 

86 

0.1196 

0.1196 

0.1196 

-0.0329 

1 

0.0005 

0.0005 

0.0005 

0  .0060 

•- 

87 

0.4162 

0.6755 

0.4162 

0.0847 

0.0007 

0.0001 

0.0007 

0.0047 

■ . 

88 

0.0121 

0.0121 

0.0121 

-0.0185 

0.0007 

0.0007 

0.0007 

0.0015 

1 

89 

0.0958 

0.0011 

0.0958 

0.0011 

0  .0958 
0.0011 

0.0546 

0.0028 

[: 

90 

0.0853 

0  .0853 

0  .0853 

-0.0037 

I., 

0.0012 

0.0012 

0.0012 

0  .0051 

’ 

91 

-0.0142 

-0.0142 

-0.0142 

-0.0379 

0.0004 

0.0004 

0.0004 

0.0012 

92 

0.1009 

0.1009 

0.1009 

-0.0276 

i 

0.0012 

0.0012 

0 ,0012 

0.0045 

93 

0.0604 

0.0604 

0.0604 

-0.0450 

0.0013 

0.0013 

0.0013 

0.0017 

r*’ 

94 

0  .0227 

0.0227 

0  .0227 

0.0369 

0.0002 

0.0002 

0.0002 

0  .00081 

r; 

95 

0.0598 

0.0598 

0  .0598 

-0  .0491 

a 

r. 

r, 

r. 

*•_ 

«*, 

J 

0.0013 

0.0013 

0.0013 
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0.0043 
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Intercept/ 

Slope 

Accepted 

Reverse 

Amino 

Ancestral 

Human 

Rabbit 

Acid 

Change 

Change 

Change 

Change 

96 

0.0239 

0.0239 

0  .0239 

-0.0081 

0.0006 

0.0006 

0.0006 

0.0018 

97 

0.1635 

0.1635 

0.1635 

0.1170 

0.0012 

0.0012 

0.0012 

0.00351 

98 

-0.0031 

-0.0031 

-0.0031 

-0.0031 

0.0001 

0.0001 

0.0001 

0.00018 

99 

0.0698 

0  .0699 

0.0699 

-0.0258 

0.0012 

0.0012 

0.0012 

0  .0044 

100 

0.1185 

0.1185 

0.1185 

-0.0490 

0.0004 

0  .0004 

0.0004 

0.0056 

101 

0.1067 

0.1067 

0.1067 

0.0462 

0.0013 

0.0013 

0.0013 

0.0029 

102 

0.0747 

0 .0747 

0.07  47 

0.0130 

0.0011 

0.0011 

0.0011 

0.0039 

103 

0.0000 

0.0000 

0  .0000 

Not  Repeated 

104 

0.0000 

0.7410 

0.0000 

0  .3444 

0.0000 

0  .3444 

0  .0753 

0.0004 

0.0005 

0.0005 

0.0015 

105 

-0.0045 

-0.0045 

-0.0045 

-0.03160 

0.0006 

0.0006 

0.0006 

0.0011 

106 

-0.0142 

-0.0142 

-0.0142 

-0.0308 

0.0012 

0.0012 

0.0012 

0.0014 

107 

0.0721 

0.0721 

0.0721 

-0.0129 

0.0012 

0.0012 

0.0012 

0.0036 

108 

0.1210 

0.1210 

0.1210 

0.0439 

0.0013 

0.0013 

0.0013 

0.0039 

109 

0.0681 

0  .0681 

0.0681 

0.0748 

0.0013 

0.0013 

0.0013 

0.0013 

110 

0.0090 

0.0090 

0.0090 

-0.0115 

0.0007 

0.0007 

0.0007 

0.0015 

111 

-0.0031 

-0.0031 

-0.0031 

-0.0091 

0.0001 

0.0001 

0.0001 

0.0002 

112 

0.3921 

0 .6656 

0  .3921 

0.0597 

0.0008 

0.0000 

0.0008 

0.0009 

113 

0.0322 

0.0322 

0.0322 

0.0325 

0.0012 

0.0012 

0.0012 

0.0013 

114 

0.1032 

0.1032 

0  .1032 

0.0650 

0.0011 

0.0011 

0.0011 

0.0019 

115 

0  .3327 

0  .3327 

0  .6867 

0.0347 

0.0001 

0.0001 

0.0003 

0.0005 

116 

0  .3263 

0  .3333 

0.3333 

0.0070 

0.0000 

0.0000 

0.0000 

0.0000 

117 

0.0406 

0.0406 

0.0406 

0  .0227 

0.0011 

0.0011 

0.0011 

0.0013 

118 

0.0467 

0  .0467 

0  .0467 

0  .0347 

0.0011 

0.0011 

0.0011 

0.0015 

Table  4.8  Continued 


Intercept/ 


Slope 

Amino 

Ancestral 

Accepted 

Human 

Rabbit 

Reverse 

Acid 

Change 

Change 

Change 

Change 

119 

0.0201 

0.0201 

0.0201 

-0.0127 

120 

0.0005 

0.0005 

0.0005 

0.0013 

0.0052 

0.0052 

0.0052 

0.0027 

121 

0.0005 

0.0005 

0.0005 

0.0005 

0.0817 

0.0817 

0.0  817 

0.0076 

122 

0.0011 

0.0011 

0.0011 

0.0028 

0.0102 

0.0102 

0.0102 

-0  .0972 

123 

0.0006 

0.0006 

0.0006 

0.0027 

0.0776 

0  .0776 

0  .0776 

0  .0034 

124 

0.0012 

0.0012 

0.0012 

0  .0044 

0.1779 

0.1779 

0.1779 

0.1586 

125 

0.0012 

0.0012 

0.0012 

0.0037 

0.7112 

0  .7073 

0.6611 

0.0811 

126 

0.0004 

0.0004 

0.0001 

0.0024 

-0.0095 

-0.0095 

“0.0095 

0.0006 

127 

0.0006 

0.0006 

0.0006 

0.0010 

0.1023 

0.1023 

0.1023 

0  .0283 

128 

0.0013 

0.0013 

0.0013 

0  .0039 

0.1918 

0.1918 

0.1918 

0  .0872 

129 

0.0010 

0.0010 

0.0010 

0.0036 

0.0812 

0.0812 

0.0812 

0.0099 

130 

0.0013 

0.0013 

0.0013 

0.0038 

0.6646 

0.3292 

0.3292 

-0.0061 

131 

0.0000 

0.0001 

0.0001 

0.001 

0.0887 

0.0887 

0  .0887 

0.0493 

132 

0.0011 

0.0011 

0.0011 

0 .0028 

0.0036 

0.0036 

0.0036 

-0.0747 

133 

0.0009 

0.0009 

0.0009 

0.0025 

0.0563 

0.0563 

0  .0563 

0  .0562 

134 

0.0007 

0.0007 

0.0007 

0.0010 

0.0000 

0.0000 

0.0000 

No  Repeats 

135 

0.0000 

0.0611 

0.0000 

0.0611 

0.0000 

0.0611 

0.0250 

136 

0.0011 

0.0011 

0.0011 

0.0020 

0.0552 

0.0552 

0  .0552 

0.0310 

137 

0.0007 

0.0007 

0.0007 

0.0017 

-0.0034 

-0.0034 

1.0034 

-0  .1809 

138 

0.0007 

0.0007 

0.0007 

0.0049 

0.1073 

0.1073 

0.1073 

-0  .0088 

139 

0.0006 

0.0006 

0.0006 

0.0054 

0.0983 

0.0983 

0  .0983 

0.1021 

140 

0.0012 

0.0012 

0.0012 

0.0016 

0.0708 

0.0708 

0.0708 

-0.0263 

141 

0.0012 

0.0012 

0.0012 

0  .0035 

0.0703 

0.0703 

0.0703 

0.0630 

0.0007 

0.0007 

0.0007 

0.0008 

Table  4.8  Continued 


Intercept/ 

Slope  Accepted  Reverse 


Amino 

Ancestral 

Human 

Rabbit 

Acid 

Chance 

Change 

Change 

Change 

142 

-0.0006 

-0.0006 

-0.0006 

-0.0118 

0.0001 

0.0001 

0.0001 

0.0004 

143 

0.0222 

0.0222 

0.0222 

0.0010 

0.0002 

0.0002 

0.0002 

0.0005 

144 

0.1261 

0.1261 

0.1261 

0.0725 

0.0012 

0.0012 

0.0012 

0.0030 

145 

0.0794 

0.0794 

0.0794 

0  .0323 

0.0013 

0.0013 

0.0013 

0.0027 

146 

0.0071 

0.0071 

0  .0071 

-0.0413 

0.0009 

0.0009 

0.0009 

0.0020 

Table  4.9  shows  the  number  of  amino  acid  residues  that 
remain  the  same  as  the  starting  hemoglobin. 


Table  4.9 

Hemoglobin  Similarity 


Starting 

Hemoglobin 

Ancestra 

Rabbit 

Human 

Combined 


Number  of  amino  acids 
Ancestral  Rabbit 

Intercept  Slope  Intercept  Slope 


135.7737  -.1239 
118.1745  -.1075 
119.1841  -.1083 
124.3775  -.1132 


118.4503  -.1057 
136.0249  -.1208 
123.5119  -.1077 
125.9957  -.1114 


similar  to: 

Human 

Intercept  Slope 
119.3220  -.1071 
123.2636  -.1078 
136.0913  -.1202 
126.2256  -.1117 


Table  4.10  takes  the  process  one  step  further  by 
presenting  the  hemoglobin  singularity  regression  for  each  of 
the  methods  used  in  this  study. 
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Table  4.10 

Hemoglobin  Similarity  Regression 


Intercept/ 

Nearest- 

Slope 

Random 

Dayhof f 

Zuckerkandl 

Fitch 

Neighbor 

Ancestral 

Starting  Hemoglobin 

Ancestral 

141.09 

131.85 

131.95 

133.12 

140.86 

-.1178 

-.1392 

-.1268 

-.1334 

-.1031 

Rabbit 

122.85 

114.99 

115.28 

116.08 

123  .05 

-.1002 

-.1181 

-.1098 

-.1137 

-.0869 

Human 

123.83 

115.78 

116.23 

116.90 

123.88 

-.1015 

-.1189 

-.1103 

-.1155 

-.0893 

Rabbit  Starting  Hemoglobin 

Ancestral 

122.81 

114.55 

115.03 

115.76 

122.72 

-.1023 

-.1194 

-.1114 

-.1147 

-.0894 

Rabbit 

141.36 

132.04 

132.95 

133.49 

140.79 

-.1139 

-.1340 

-.1256 

-.1137 

-.1014 

Human 

127.79 

119.66 

120.08 

121.08 

127.69 

-.1016 

-.1204 

-.1120 

-.1159 

-.0892 

Human  Starting  Hemoglobin 

Ancestral 

123 .86 

115.55 

115.96 

116.80 

123.75 

-.1032 

-.1210 

-.1110 

-.1167 

-.0894 

Rabbit 

127.97 

119.90 

120.43 

121.43 

127.83 

-.1008 

-.1213 

-.1119 

-.1165 

-.0881 

Human 

141.27 

132.10 

132.04 

133.56 

141.15 

-.1146 

-.1341 

-.1235 

-.1292 

-.0991 

Table  4.11  presents  the  final  tabular  results  examined 
in  this  effort.  The  amino  acid  residues  similarity  to  the 
starting  hemoglobins  is  used  in  a  regression  by  iterations 
for  each  position  of  each  of  the  three  starting  hemoglobins. 

Table  4.11 

Amino  Acid  Similarity  Regression 
Amino  Ancestral  Human  Rabbit 

Acid  Intercept  Slope  Intercept  Slope  Intercept  Slope 


Table  4.11  Continued 


Amino  Ancestral  Human  Rabbit 


Acid 

Intercept  Slope 

Intercept  Slope 

Intercept  Slope 

7 

.8475 

-.0010 

.8475 

-.0010 

.8475 

-.0010 

8 

.9130 

-.0012 

.9130 

-.0012 

.9130 

-.0012 

9 

.3259 

-.0002 

.6377 

-.0001 

.6377 

-.0001 

10 

.8845 

-.0011 

.8845 

-.0011 

.8845 

-.0011 

11 

.3330 

-.0001 

.6306 

-.0005 

.6306 

-.0005 

12 

1.0000 

0.0000 

1.0000 

0  .0000 

1.0000 

0.0000 

13 

.3312 

-.0004 

.6420 

-.0002 

.6420 

-.0002 

14 

.9644 

-.0005 

.9644 

-.0005 

.9644 

-.0005 

15 

.9365 

-.0008 

.9365 

-.0008 

.9365 

-.0008 

16 

.9112 

-.0008 

.9112 

-.0008 

.9112 

-.0008 

17 

.9444 

-.0011 

.9444 

-.0011 

.9444 

-.0011 

18 

.9830 

-.0012 

.9830 

-.0012 

.9830 

-.0012 

19 

.8651 

-.0013 

.8651 

-.0013 

.8651 

-.0013 

20 

.9195 

-.0010 

.9195 

-.0010 

.9195 

-.0010 

21 

.5692 

-.0008 

.2635 

-.0004 

.5692 

-.0008 

22 

1.0031 

-.0001 

1.0031 

-.0001 

1.0031 

-.0001 

23 

.8238 

-.0010 

.8238 

-.0010 

.8238 

-.0010 

24 

.8941 

-.0007 

.8941 

-.0007 

.8941 

-.0007 

25 

.8435 

-.0012 

.8435 

-.0012 

.8435 

-.0012 

26 

.8030 

-.0010 

.8030 

-.0010 

.8030 

-.0010 

27 

.9238 

-.0011 

.9238 

-.0011 

.9238 

-.0011 

28 

.9509 

-.0009 

.9509 

-.0009 

.9509 

-.0009 

29 

.8235 

-.0011 

.8235 

-.0011 

.8235 

-.0011 

30 

.2721 

-.0004 

.6260 

-.0008 

.6260 

-.0008 

31 

.9830 

-.0003 

.9830 

-.0003 

.9830 

-.0003 

32 

.9810 

-.0007 

.9810 

-.0007 

.9810 

-.0007 

33 

.9830 

-.0004 

.9830 

-.0004 

.9830 

-.0004 

34 

.9116 

-.0014 

.9116 

-.0014 

.9116 

-.0014 

35 

.9970 

-.0001 

.9970 

-.0001 

.9970 

-.0001 

36 

.7840 

-.0012 

.7840 

-.0012 

.7840 

-.0012 

37 

.9558 

-.0012 

.9558 

-.0012 

.9558 

-.0012 

38 

.92  81 

-.0010 

.92  81 

-.0010 

.92  81 

-.0010 

39 

.8206 

-.0008 

.8206 

-.0008 

.8206 

-.0008 

40 

.9234 

-.0008 

.9234 

-.0008 

.9234 

-.0008 

41 

.9496 

-.0007 

.9496 

-.0007 

.9496 

-.0007 

42 

.8543 

-.0011 

.85  43 

-.0011 

.85  43 

-.0011 

43 

.3395 

-.0004 

.6146 

-.0007 

.6146 

-.0007 

44 

1.0121 

-.0012 

1.0121 

-.0012 

1.0121 

-.0012 

45 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

46 

.92  86 

-.0009 

.9286 

-.0009 

.92  86 

-.0009 

47 

.9628 

-.0006 

.9628 

-.0006 

.9628 

-.0006 

48 

1.0121 

-.0012 

1.0121 

-.0012 

1.0121 

-.0012 

49 

.9854 

-.0000 

.9854 

-.0000 

.9854 

-.0000 

50 

.6667 

-.0000 

.3333 

-.0000 

.6667 

-.0000 

51 

.6247 

-.0002 

.2946 

-.0003 

.6247 

-.0002 

52 

.6366 

-.0002 

.6366 

-.0002 

.2884 

-.0004 

53 

.9194 

-.0004 

.9194 

-.0004 

.9194 

-.0004 

54 

.9733 

-.0001 

.9733 

-.0001 

.9733 

-.0001 
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I 


I 


Amino 

Ancestral 

Human 

Rabbit 

Acid 

Intercept  Slope 

Intercept  Slope 

Intercept  Slope 

55 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

56 

.5990 

-.0007 

.5990 

-.0007 

.3438 

-.0001 

57 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

58 

.9063 

-.0010 

.9063 

-.0010 

.9063 

-.0010 

59 

.7966 

-.0012 

.7966 

-.0012 

.7966 

-.0012 

60 

.9718 

-.0008 

.9718 

-.0008 

.9718 

-.0008 

61 

.9725 

-.0009 

.9725 

-.0009 

.9725 

-.0009 

62 

.9177 

-.0011 

.9177 

-.0011 

.9177 

-.0011 

63 

.9915 

-.0012 

.9915 

-.0012 

.9915 

-.0012 

64 

.8265 

-.0012 

.8265 

-.0012 

.8265 

-.0012 

65 

.8558 

-.0008 

.8558 

-.0008 

.8558 

-.0008 

66 

.9701 

-.0006 

.9701 

-.0006 

.9701 

-.0006 

67 

.8903 

-.0012 

.8903 

-.0012 

.8903 

-.0012 

68 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

69 

.6173 

-.0008 

.6173 

-.0008 

.3351 

-.0001 

70 

1.0290 

-.0012 

1.0290 

-.0012 

1.0290 

-.0012 

71 

.8859 

-.0013 

.8859 

-.0013 

.8859 

-.0013 

72 

.3340 

-.0002 

.6655 

-.0001 

.6655 

-.0001 

73 

.5582 

-.0008 

.5582 

-.0008 

.3788 

-.0001 

74 

.9140 

-.0012 

.9140 

-.0012 

.9140 

-.0012 

75 

1.0112 

-.0010 

1.0112 

-.0010 

1.0112 

-.0010 

76 

.2940 

-.0040 

.2991 

-.0001 

.2907 

-.0004 

77 

.8775 

-.0013 

.8775 

-.0013 

.8775 

-.0013 

78 

.9657 

-.0004 

.9657 

-.0004 

.9657 

-.0004 

79 

.8073 

-.0012 

.8073 

-.0012 

.8073 

-.0012 

80 

.8931 

-.0012 

.8931 

-.0012 

.8931 

-.0012 

81 

.9540 

-.0006 

.9540 

-.0006 

.9540 

-.0006 

82 

.8131 

-.0012 

.8131 

-.0012 

.8131 

-.0012 

83 

.8376 

-.0008 

.8376 

-.0008 

.8376 

-.0008 

84 

.8001 

-.0009 

.8001 

-.0009 

.8001 

-.0009 

85 

.9737 

-.0008 

.9737 

-.0008 

.9737 

-.0008 

86 

.8804 

-.0005 

.8804 

-.0005 

.8804 

-.0005 

87 

.5838 

-.0007 

.3245 

-.0001 

.5838 

-.00.07 

88 

.9879 

-.0007 

.9879 

-.0007 

.9879 

-.0007 

89 

.9042 

-.0011 

.9042 

-.0011 

.90  42 

-.0011 

90 

.9147 

-.0012 

.9147 

-.0012 

.9147 

-.0012 

91 

1.0142 

-.0004 

1.0142 

-.0004 

1.0142 

-.0004 

92 

.8991 

-.0012 

.8991 

-.0012 

.8991 

-.0012 

93 

.93  96 

-.0013 

.93  96 

-.0013 

.93  96 

-.0013 

94 

.9773 

-.0002 

.9773 

-.0002 

.9773 

-.0002 

95 

.9402 

-.0013 

.9402 

-.0013 

.9402 

-.0013 

96 

.9761 

-.0006 

.9761 

-.0006 

.9761 

-.0006 

97 

.8365 

-.0012 

.8365 

-.0012 

.8365 

-.0012 

98 

1.0031 

-.0001 

1.0031 

-.0001 

1.0031 

-.0001 

99 

.9302 

-.0012 

.9301 

-.0012 

.9301 

-.0012 

100 

.8815 

-.0004 

.8815 

-.0004 

.8815 

-.0004 

101 

.8933 

-.0013 

.8933 

-.0013 

.8933 

-.0013 

102 

.9253 

-.0011 

.9253 

-.0011 

.9253 

-.0011 
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Table  4.11  Continued 


Amino 

Ancestral 

Human 

Rabbit 

Acid 

Intercept  Slope 

Intercept  Slope 

Intercept  Slope 

103 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

104 

.2590 

-.0004 

.6556 

-.0005 

.6556 

-.0005 

105 

1.0045 

-.0006 

1.0045 

-.0006 

1.0045 

-.0006 

106 

1.0142 

-.0012 

1.0142 

-.0012 

1.0142 

-.0012 

107 

.9279 

-.0012 

.9279 

-.0012 

.9279 

-.0012 

108 

.8790 

-.0013 

.8790 

-.0013 

.8790 

-.0013 

109 

.9319 

-.0013 

.9319 

-.0013 

.9319 

-.0013 

110 

.9910 

-.0007 

.9910 

-.0007 

.9910 

-.0007 

111 

1.0031 

-.0001 

1.0031 

-.0001 

1.0031 

-.0001 

112 

.6079 

-.0008 

.3344 

-.0000 

.6079 

-.0008 

113 

.9678 

-.0012 

.9678 

-.0012 

.9678 

-.0012 

114 

.8968 

-.0011 

.8968 

-.0011 

.8968 

-.0011 

115 

.6673 

-.0001 

.6673 

-.0001 

.3133 

-.0003 

116 

.3263 

-.0000 

.6667 

-.0000 

.6667 

-.0000 

117 

.9594 

-.0011 

.9594 

-.0011 

.9594 

-.0011 

118 

.9532 

-.0011 

.9532 

-.0011 

.9532 

-.0011 

119 

.9799 

-.0005 

.9799 

-.0005 

.9799 

-.0005 

120 

.9948 

-.0005 

.9948 

-.0005 

.9948 

-.0005 

121 

.9183 

-.0011 

.9183 

-.0011 

.9183 

-.0011 

122 

.9898 

-.0006 

.9898 

-.0006 

.9898 

-  .0006 

123 

.9224 

-.0012 

.9224 

-.0012 

.9224 

-.0012 

124 

.8221 

-.0012 

.8221 

-.0012 

.8221 

-.0012 

125 

.2888 

-.0004 

.2922 

-.0004 

.3389 

-.0001 

126 

1.0095 

-.0006 

1.0095 

-.0006 

1.0095 

-.0006 

127 

.8977 

-.0013 

.8977 

-.0013 

.8977 

-.0013 

128 

.80  82 

-.0010 

.8082 

-.0010 

.80  82 

-.0010 

129 

.9188 

-  .0013 

.9188 

-.0013 

.9188 

-.0013 

130 

.3354 

-.0000 

.6708 

-.0001 

.6708 

-.0001 

131 

.9113 

-.0011 

.9113 

-.0011 

.9113 

-.0011 

132 

.9964 

-.0009 

.9964 

-.0009 

.9964 

-.0009 

133 

.9437 

-.0007 

.9437 

-.0007 

.9437 

-.0007 

134 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

135 

.9389 

-.0011 

.9389 

-.0011 

.9389 

-.0011 

136 

.9448 

-.0007 

.9448 

-.0007 

.9448 

-.0007 

137 

1.0034 

-.0007 

1.0034 

-.0007 

1.0034 

-.0007 

138 

.8927 

-.0006 

.8927 

-.0006 

.8927 

-.0006 

139 

.9017 

-.0012 

.9017 

-.0012 

.9017 

-.0012 

140 

.9292 

-.0012 

.9292 

-.0012 

.9292 

-.0012 

141 

.9297 

-.0007 

.9297 

-.0007 

.9297 

-.0007 

142 

1.0006 

-.0001 

1.0006 

-.0001 

1.0006 

-.0001 

143 

.9778 

-.0002 

.9778 

-.0  02 

.9778 

-.0002 

144 

.8739 

-.0012 

.8739 

-.0012 

.8739 

-.0012 

145 

.9206 

-.0013 

.9206 

-.0013 

.9206 

-.0013 

146 

.9919 

-.0009 

.9919 

-.0009 

.9919 

-.0009 

CHAPTER  V 


CONCLUSIONS 

The  general izability  of  results  from  this  work  is 
restricted  due  to  the  limited  availability  of  sequence  data 
for  protein  generating  sequences.  The  DNA  sequencing  data 
for  human  alpha,  beta,  and  sigma  hemoglobin;  for  mouse 
alpha  and  beta  hemoglobin;  and  for  rabbit  alpha  and  sigma 
hemoglobin  is  currently  known.  Until  the  sequencing  data 
from  a  wide  variety  of  different  proteins  from  a  large 
cross-section  of  different  species,  the  full  implication  of 
the  use  of  these  techniques  will  be  difficult  to  completely 
evaluate. 

The  review  of  the  literature  revealed  a  great  deal  of 
literature  relative  to  the  possible  shemes  for  the  evolution 
of  DNA.  Literature  from  basic  fields  of  study  have  not  yet 
converged.  Literature  being  produced  in  the  field  of 
population  genetics  have  addressed  the  dynamics  of 
phenotypic  introduction  and  propogation  of  mutations  in  the 
genetic  material  of  the  organism.  The  general  results  from 
the  work  in  the  field  to  date  is  that  mutations  resulting  in 
neutral,  or  nearly  nei:.  ral  advantage  to  the  survival  of  an 
organism  in  the  environment  are  not  supportable  from  the 
frequency  of  mutations  observed  when  compared  to  the  best 


models  of  the  rate  of  acceptance  of  this  type  of  mutation. 
This  effort  found  in  almost  all  cases  an  order  of  magnitude 
higher  rejection  than  acceptance,  even  using  the  extremely 
lax  criteria  used  in  this  effort  for  acceptance. 

The  molecular  biologists  have  found  evidence  for 
several  effects  that  suggest  that  many  of  the  mutations  that 
are  expected  to  be  found  in  the  primary  sequence  of  amino 
acids  in  protein  should  result  in  selectively  neutral 
mutations.  Further,  studies  which  relate  the  differences  in 
the  primary  sequence  of  proteins  such  as  hemoglobin, 
cytochrome  c,  and  fibrinogin  suggest  a  constant  rate  of 
amino  acid  replacement  for  each  protein  which  provides 
evolutionary  relationships  between  species  which  agrees  well 
with  the  phylogenic  relationships  that  have  been  previously 
established.  This  strongly  suggests  that  neutral  mutations 
would  be  expected  to  be  the  primary  mechanism  envolved. 
Some  of  the  assumptions  used  for  estimating  the  transition 
probabilities  of  amino  acid  may  be  suspect.  Use  of  the 
nucleotides  sequences  instead  of  the  amino  acid  sequences 
as  the  methodology  adapted  in  this  effort  should  avoid  these 
same  difficulties  and  was  the  motivating  force  for  the 
selection  ofthe  modeling  techniquechosen  in  this  effort. 

Information  theory  and  the  application  of  information 
theory  have  been  applied  to  information  on  the  evolutionary 
system.  Involvement  of  this  field  to  the  problems  existing 
in  the  genetic  system  have  not  had  a  large  impact  to  da_e. 


The  genetic  code  and  the  reduncies  of  the  system  have  been 
successfully  evaluated.  Application  to  the  mutation  of 
amino  acid  sequences  have  not  provided  significant  insight 
into  the  behavior  of  the  system.  More  basic  information 
provided  from  nucleotide  mutation  behavior  should  allow 
further  analysis  in  this  area  with  a  better  prospect  for 
success.  Recent  information  that  has  been  obtained  in 
studies  of  the  chemical  and  physical  characteristics  of  the 
nucleic  acid  found  in  organisms  has  changed  many  ideas  about 
possible  interactions  that  may  be  occurring.  The 
observation  of  multiple  repeated  sequences  of  nucleotides  in 
the  DN A  of  eukaryotes  and  the  finding  of  intervening 
sequences  within  the  gene  that  codes  for  DNA  will  have  to  be 
understood  and  incorporated  into  future  work. 

Also,  the  field  of  chemical  physics  has  yielded 
information  on  chain  stability  caused  by  the  ratio  of 
guanine  and  cytosine  and  the  stabilization  of  mRNA  by  the 
formation  of  autocomplementary  regions  and  the  formation  of 
loops  offer  potential  for  understanding  the  selective 
advantage  of  nucleotide  changes  that  result  in  amino  acid 
point  mutations  which  do  not  appear  to  have  selective 
advantage  when  examined  on  the  basis  of  protein  function. 

Thefields  of  study  involved  in  trying  to 
understand  the  underlying  mechanisms  does  not  appear  to  havt 
been  incorporated  into  a  systematic  framework  necessary  for 
the  evaluation  of  system  performance.  This  thesis  effort  is 
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an  attempt  to  try  to  evaluate  the  system  by  examining  the 
output  produced  by  several  different  sets  of  inputs  with 
different  methods  of  modeling  system  performance,  since  the 
underlying  mechanisms  were  uncharacterized.  Estimates  of 
the  transition  probabilites  were  obtained  from  three 
different  sources  as  well  as  assumed  to  be  completely 
random.  The  output  of  the  system  was  characterized  at  the 
nucleotide  level  and  at  the  level  of  the  amino  acid.  The 
output  of  the  system  was  "locked"  in  place  by  constraining 
the  acceptable  amino  acids  at  each  position  of  the  chain  to 
those  residues  that  have  been  characterized  in  previous 
studies.  This  artificial  constraint  was  an  attempt  to 
restrict  the  model  to  a  domain  that  would  be  evolutionary 
significant.  In  other  words,  this  was  a  primative  attempt 
to  collectively  model  the  possibility  or  random  mutation 
resulting  in  selective  advantage  at  the  nucleotide  and  the 
amino  acid  level.  As  a  better  understanding  of  mechanisms 
that  have  physical  significance  is  formulated,  more  detailed 
modeling  of  the  process  will  be  possible. 

There  were  several  interesting  observations  that  can  be 
made  regarding  the  results  of  the  modeling  effort.  All  of 
the  transition  frequency  were  based  upon  estimates  of  an 
accepted  amino  acid  change  in  10,000,000  years  for  a  protein 
140  residues  long..  The  data  in  table  4.6  reflects  that 
using  the  Dayhoff  transition  frequencies  resulted  in  the 
highest  acceptance  of  new  amino  acid  residues  while  the 
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nearest-neighbor  transition  frequencies  resulted  in  the 
lowest  acceptance  rate.  The  random  data  was  nearly 
identical  with  the  nearest-neighbor  rate.  All  three  simple 
transition  methodologies,  Dayhoff,  Zuckerkandl,  and  Pitch 
were  very  close  to  having  the  same  rate  of  accepting  a 
mutation.  The  rate  at  which  a  previous  amino  acid  residue 
was  back  mutated  to  a  residue  that  was  already  at  the  same 
position  was  very  high  and  uniform  across  methods.  99%  of 
all  accepted  mutations  using  the  random  transition 
probabilities  resulted  in  a  back  mutation.  97%,  96%,  and 
95%  were  the  respective  back  mutation  rate  for  Dayhoff, 
Zuckerkandl,  and  Fitch.  Nearest-neighbor  back  mutations 
accounted  for  93%  of  the  accepted  mutations.  This  extremely 
high  percentage  of  back  mutations  at  the  amino  acid  level 
was  not  expected  before  this  analysis  was  begun.  Table  4.4 
for  the  nucleotide  transition  regression  sheds  some  light  on 
the  situation.  All  five  methods  resulted  in  the  acceptance 
of  only  15%  to  16%  of  the  nucleotide  changes.  Of  the 
nucleotide  changes  that  did  occur,  17%,  44%,  42%,  42%,  and 
19%  of  the  changes  were  backmutations.  The  rates  of 
transition  are  nearly  the  same  for  allowing  nucleotide 
changes  as  for  accepting  amino  acid  changes.  At  first 
thought,  it  would  appear  that  some  form  of  the  wobble 
hypothesis  of  Crick  might  provide  the  answer.  The  wobble 
hypothesis  has  been  applied  to  the  degenerate  pairing  of  the 
third  nucleotide  of  mRNA  to  tRNA  in  a  manner  which  preserves 


the  reliability  of  obtaining  the  same  amino  acid  in  the 
peptide  chain.  When  the  data  from  table  4.7  is  examined,  no 
significant  trend  of  either  allowing,  repeating,  or 
rejecting  nucleotide  based  upon  the  position  in  the  codon 
can  be  found  to  support  a  wobble  relationship  in  the 
mutations  at  the  DNA  level.  Upon  further  investigation,  the 
relationship  between  nucleotide  transition  frequency  and 
amino  acid  transition  frequency  are  tied  in  two  cycles.  The 
amino  acids  acceptable  at  any  particular  position  typically 
form  a  small  subset  of  the  20  residues  possible.  As  a 
result,  nucleotides  changes  are  allowed  at  any  of  the 
positions  in  the  codon  with  few  restrictions,  but  after  a 
very  few  resulting  changes  in  amino  acids,  they  begin 
repeating.  The  repetive  cycling  of  amino  acids  was  much 
higher  than  anticipated  by  any  of  the  literature.  Most  of 
the  literature  dealing  with  the  evolutionary  fixation  of 
mutations  in  amino  acids  did  not  reflect  the  high  degree  of 
repetition  that  was  found  to  exist  at  the  amino  acid  level. 

Interestingly  enough,  the  method  used  in  the  simulation 
proved  to  be  the  driver  for  the  type  of  repetition  seen. 
The  repetition  cycles  were  very  dependent  on  the  position. 
At  one  end  of  the  spectrum  was  the  use  of  the  random  method. 
Most  of  the  repeats  observed  were  simple  back  mutations  at 
the  third  position  in  the  DNA  codon.  As  the  transitions  of 
the  nucleotides  is  more  controlled  such  as  in  the  nearest- 


neighbor  approach,  repeats  of  amino  acids  occur  more  often 
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only  after  cycles  of  two,  three,  or  somtimes  even  more 
changes  in  the  amino  acid  residue.  These  cycles  are  found 
to  regularly  repeat.  The  structured  transition 
probabilities  of  the  three  methods  seemed  to  be  intermediate 
in  behavior.  Some  two  and  three  cycle  transition  patterns 
were  seen,  but  they  were  not  regularily  repeated  as  were  the 
transitions  observed  with  the  nearest-neighbor  method. 

This  difference  in  mutational  patterns  depending  on  the 
non-random  transition  probabilities  seems  to  paint  a 
different  picture  of  the  process  of  evolution.  The  fixation 
of  a  point  mutation  may  be  part  of  an  equil librium  searching 
mechanism  important  in  the  survival  of  the  organism.  The 
strength  of  the  Ising  model  used  with  the  nearest-neighbor 
transition  probabilities  and  the  allowable  amino  acid 
replacements  may  enable  the  organism  to  control  random 
genetic  drift  and  neutral  mutations  at  the  protein  level  to 
control  the  stability  of  the  DNA  and  the  RNA  chains.  The 
ability  to  control  the  guanine-cytosine  content  of  DNA  is 
important  in  the  information  potential  of  DNA.  It  is  also 
important  in  the  stablity  of  the  chain  to  denaturalization 
from  heat  as  well  as  ion  levels.  It  controls  the  rate  of 
chain  separation.  These  physical  properties  of  the  DNA 
chain  have  been  shown  to  follow  the  same  behavior  and 
performance  as  would  be  expected  from  an  Ising  model. 

Much  further  work  needs  to  be  done.  Further 
identification  of  more  DNA  from  more  species  and  for 
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different  proteins  needs  to  be  conducted.  With  a  larger 
database,  more  closure  on  the  process  of  evolution  may  be 
examined.  DNA  from  species  of  very  distant,  or  non- 
existant,  connection  to  the  environment  would  be 
instrumental  in  identifying  several  points.  An  opportunity 
to  use  DNA  sequences  from  species  found  around  the  hydrogen 
sulfide  vents  on  the  ocean  floor  in  the  Atlantic  and  the 
Pacific  oceans  has  high  potential  to  shed  light  on  many  of 
the  theories  existing  in  the  field  of  evolutionary  biology. 


PROGRAM  SOURCE  LISTING 


00005 

00010 

00020 

00030 

00040 

00050 

00060 

00061 

00070 

00080 

00081 

00090 

00110 

00120 

00130 

00190 

00191 

00192 

00193 

00194 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00360 

00370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 


PROGRAM  (COMPUTERANDOM  MUTATIONS  FOR  BETA  Hb) 
INTEGER  E,F,K 

DIMENSION  NA( 438,50) ,MA(146 ,50) 

DIMENSION  NULL(438) ,PCT(438) 

DIMENSION  NCOUNT (4,4) , MCOUNT (28 , 2 8) 

DIMENSION  MCRIT (146 ,14) 

DIMENSION  MGC (16 ,4) fN(450) ,M(150) ,NP(450 ,4)  , 

IMP (150 ,4) 

DIMENSION  T(64,3) 

DIMENSION  NR ( 43 8)  ,MR(150)  , MCH (150)  fMCR(150)  , 
1MCA (150) 

IXS=9636 
CALL  HB(NP) 

CALL  GC (MGC) 

CALL  INIT1 (NP,MP,MGC) 

CALL  AAAR (MCRIT) 

DO  1001  INTE=1,3 
IF  (INTE.EQ.l)  INTER=7  00 
IF  (INTE.EQ.2)  INTER=40 
IF  (INTE.EQ.3)  INTER=6 

CALL  INIT ( IX , IXS , F , E , N , INTER , NCOUNT  ,  MCOUNT , NULL 
1NA, MA , M  , MGC , NP , NIP , NR , MR , MP ,  IH  , MCH , MCR , MCA) 

DO  10  L=1 , INTER 
DO  20  1=2,438 

CALL  TRAN1 (NUC ,  I , INTER, N, IX) 

IF  (NUC.EQ.N(I) )  GO  TO  20 

CALL  MUTATE ( I , MCRIT , NULL , NUC , NCOUNT , MCOUNT , NA , 
1MA,N,M,MGC, E,F,NR,MR) 

20  CONTINUE 
10  CONTINUE 
PRINT  2 
METH=1 

CALL  OUT ( M , N , NA , MA , NULL , MCOUNT , NCOUNT , I H , IXS , 
1METH , INTER , NR , MR , MCH , MCR , MCA , MP ) 

DO  30  K=1 ,3 

CALL  INIT ( IX , IXS , F , E , N , INTER , NCOUNT, MCOUNT , NULL 
1NA,MA,M,MGC ,NP ,NIP ,NR,MR,MP, I H, MCH, MCR, MCA) 

DO  40  L=l, INTER 
DO  50  1=2,438 

CALL  TRAN2 ( NUC , I , INTER , N , K , IX ) 

IF  (NUC.EQ. (N (I) ) )  GO  TO  50 

CALL  MUTATE ( I , MCRIT , NULL , NUC , NCOUNT , MCOUNT , NA , 
1MA,N,M,MGC,E,F,NR,MR) 

50  CONTINUE 
40  CONTINUE 

IF  (K.EQ.l)  PRINT  3 


00460 
00470 
00480 
00490 
00500 
00510 
00520 
00530 
00540 
00550 
00560 
00570 
00580 
00590 
00600 
00610 
00620 
00630 
00640 
00650 
00660 
00670 
00680 
00690 
00700 
00710 
00715 
00750 
00760 
00761 
00770 
00780 
00781 
00790 
00  810 
00820 
10010 
10020 
10030 
10040 
10050 
10060 
10070 
10110 
10120 
10130 

10140 

10141 

10142 
10150 
10160 
10169 


IF 

(K.EQ.l) 

METH=2 

IF 

(K.EQ.2) 

PRINT  4 

IF 

(K.EQ.2) 

METH=3 

IF 

(R.EQ.3) 

PRINT  5 

IF 

(K.EQ.3) 

METH=4 

C  CALL  RESULT (M , N , NA , MA , NULL , MCOUNT , NCOUNT) 

CALL  OUT (M,N,NA,MA, NULL , MCOUNT , NCOUNT , I H , IXS , 

1 METH , INTER, NR, MR, MCH , MCR , MCA , MP ) 

30  CONTINUE 

CALL  INIT(IX,IXS ,F,E,N, INTER , NCOUNT , MCOUNT , NULL , 
1NA,MA,M,MGC,NP,NIP,NR,MR,MP, IH,MCH ,MCR,MCA) 

CALL  TP (T) 

DO  60  L=1 , INTER 
DO  70  1=2 ,438 

CALL  TRAN3 (NUC , I , INTER,N,T, IX) 

IF  (NUC.EQ. (N(I) ) )  GO  TO  70 

CALL  MUTATE ( I , MCR I T , NULL ,NUC, NCOUNT , MCOUNT , NA , 
1MA,N,M,MGC, E,F,NR,MR) 

70  CONTINUE 
60  CONTINUE 
PRINT  6 
METH=5 

C  CALL  RESULT ( M ,  N , NA , MA , NULL , MCOUNT , NCOUNT ) 

CALL  OUT ( M ,  N , NA , MA , NULL , MCOUNT , NCOUNT , IH, IXS , 
1METH , INTER , NR , MR , MCH , MCR , MCA , MP } 

1001  CONTINUE 
1000  CONTINUE 

2  FORMAT  (16H0RANDUM  MUTATION) 

3  FORMAT  (42HONUCLEOTIDE  DEPENDENT  MUTATION- 
1 Z UK ERRAND L ) 

4  FORMAT  (36HONUCLEOTIDE  DEPENDENT  MUTATION-FITCH) 

5  FORMAT  (38HONUCLEOTIDE  DEPENDENT  MUTATION- 
1DAYHOFF) 

6  FORMAT  (36H0NEAREST  NEIGHBOR  DEPENDENT  MUTATION) 
STOP 

END 

SUBROUTINE  HB(NP) 

DIMENSION  NP ( 450 / 4) 

DO  10  1=1,450 
NP  ( 1 , 1 )  =0 
NP ( I , 2 ) =0 
NP  (1 ,3)  =0 
10  CONTINUE 
NP ( 436 ,1) =2 
NP( 437,1) =1 
NP ( 43  8 ,1) =2 
NP  (433  ,1)  =3 
NP  ( 433 , 2)  =1 
NP  ( 433 ,3)  =1 
NP (434 ,1) =4 
NP (435,1) =3 
NP ( 430 ,1) =2 


81 


i 


10170 

NP(431,1)=1 

10180 

NP  ( 432,1)  =3 

10181 

NP(432 ,2) =2 

10182 

NP (432,3) =2 

10190 

NP  (427 ,1)  =1 

10200 

NP (428 ,1) =3 

10210 

NP ( 429 ,1) =4 

10220 

NP (424,1) =1 

10230 

NP  ( 425 ,1)  =3 

10231 

NP  (425 ,2)  =2 

10240 

NP  ( 426 , 1)  =3 

10241 

NP (426 ,2) =4 

10242 

NP  ( 426  ,3)  =2 

10250 

NP  (421 ,1)  =2 

10260 

NP ( 422 ,1) =4 

10270 

NP  (423 ,1) =2 

10280 

NP ( 418 ,1) =2 

10290 

NP (419,1) =4 

10300 

NP (420 ,1) =2 

10310 

NP(415,1) =2 

10320 

NP (416 ,1) =4 

10330 

NP (417,1) =4 

10340 

NP ( 412 ,1) =1 

10350 

NP  ( 413  ,1)  =3 

10360 

NP ( 414 ,1) =1 

10362 

NP ( 414 ,3) =2 

10370 

NP ( 409 ,1) -2 

10380 

NP  (410 ,1)  =3 

10390 

NP (411 ,1) -2 

10400 

NP  (406 ,1)  =1 

10401 

NP  ( 406 , 2)  =3 

10410 

NP  ( 407 ,1)  =1 

10420 

NP  (408,3)  =2 

10422 

NP  (40  8 ,3 )  =4 

10430 

NP  ( 403  ,1)  =1 

10440 

NP (404,1) =3 

10450 

NP  ( 405 , 1)  =4 

10460 

NP  (400 ,1)  =3 

10470 

NP  ( 401 , 1)  *3 

10480 

NP (402 ,1) =2 

10482 

NP  { 402  ,3)  «1 

10490 

NP (397 ,1) =2 

10500 

NP(398,1) =1 

10510 

NP (399,1) =3 

10520 

NP (3  94 ,1) =2 

10530 

NP (395 ,1) =2 

10540 

NP  (3  96  ,1)  =1 

10550 

NP  (3  91 ,1)  =3 

10560 

NP (392 ,1) =2 

10570 

NP (393 ,1) =2 

10580 

NP (388 ,1) =2 

10590 

NP (389 ,1) =4 
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10600 

10610 

10620 

10630 

10640 

10641 

10642 
10650 
10660 
10670 
10680 
10690 

10700 

10701 

10702 
10710 
10720 
10730 
10740 
10750 
10760 
10770 
10780 
10790 
10800 
10810 
10820 
10830 
10840 
10850 
10860 
10870 
10880 
10890 
10900 
10910 
10920 
10930 
10940 
10950 
10960 
10970 
10980 
10981 
10990 
11000 
11010 
11020 
11030 
11040 
11050 
11060 


NP(390  ,1)=4 
NP(385,1) =2 
NP(386,1)  =1 
NP (3  87  ,1) =2 
NP  (3  82  ,1)  =3 
NP  ( 3  82 , 2 )  =1 
NP (382 ,3) =1 
NP (383  fl) =4 
NP(384 ,1)  =4 
NP(379,1) =2 
NP(380,1)  =1 
NP  (3  81 ,1)  =2 
NP (376 ,1) =1 
NP(376,2)=4 
NP ( 376 ,3) =4 
NP (377  fl) =4 
NP(378 ,1) -2 
NP (373 ,1) =4 
NP(374,1)=4 
NP (375,1) =2 
NP (370 ,1) =1 
NP (371,1) =1 
NP(372 ,1) =2 
NP  (367 ,1) =1 
NP (368 ,1) =2 
NP (369,1) =2 
NP (364 ,1) *1 
NP (365 ,1) «2 
NP(366 ,1) =2 
NP ( 3  6 1 , 1 )  =2 
NP (362 ,1) =4 
NP (363 ,1) =2 
NP (358 ,1) =3 
NP  ,359 ,1) =3 
NP(360,1)=2 
NP (355 ,1) -2 
NP(356 ,1) =1 
NP (357 ,1) =3 
NP (352 ,1) =3 
NP (353 ,1) =2 
NP (354 ,1) =2 
NP ( 3  4  9 , 1 ) =2 
NP(350 ,1) =1 
NP (350 ,2) =2 
NP(351 ,1) =4 
NP  (346,1)  ==2 
NP(3  47 ,1)  *1 
NP ( 3  4  8 , 1 ) *3 
NP(343 ,1)  *2 
NP (344,1) =1 
NP (345 ,1) =3 
NP (340 ,1) =2 


11061 

NP(340,2)=1 

11062 

NP (340 ,3) =1 

11070 

NP (341 ,1) =1 

11080 

NP(342,1)=2 

11090 

NP(337,1)=3 

11100 

NP (338,1) =1 

11110 

NP (339 ,1) =2 

11120 

NP ( 3  3  4 , 1 ) =3 

11130 

NP(335,1)=4 

11140 

NP (336 ,1) =1 

11150 

NP(331,1)=1 

11160 

NP (332,1) =3 

11170 

NP(333 ,1) =3 

11180 

NP (328, 1) =2 

11190 

NP ( 329 ,1) =2 

11200 

NP (330 ,1) =1 

11210 

NP  ( 3  2  5 , 1 )  =3 

11220 

NP (326 ,1) =3 

11230 

NP (327 ,1) =4 

11240 

NP (322 ,1)  =2 

11250 

NP(323 ,1) =4 

11260 

NP (324 ,1) =3 

11270 

NP (319 ,1) =2 

11280 

NP (320 ,1) =2 

11290 

NP (321 ,1) =4 

11300 

NP (316 ,1) =3 

11310 

NP (317 ,1) =1 

11320 

NP (318,1) =1 

11330 

NP (313 ,1) ”1 

11331 

NP (313 -3 

113  40 

NP  (314  ,i,  =1 

11350 

NP (315 ,1) =1 

11360 

NP (310 ,1) =2 

11362 

NP (310 ,3) *3 

11370 

NP (311 ,1) =4 

113  80 

NP (312 ,1 ) =2 

113  90 

NP (307 ,1) =3 

11400 

NP (308,1) =3 

11410 

NP (309 ,1) =1 

11420 

NP (304 ,1) =1 

11430 

NP (305 ,1) =1 

11440 

NP  (306 ,1) =1 

11450 

NP (301 ,1) =2 

11460 

NP (302,1) =2 

11470 

NP (303 ,1) =2 

11480 

NP (29€ ,1) =1 

11481 

NP(298 ,2) =3 

11490 

NP (299,1) =4 

11500 

NP(300 ,1) =2 

11510 

NP (295 ,1) =2 

11520 

NP(2  96 ,1)  =1 

11530 

NP (297 ,1) =3 

11540 

11550 

11560 

11570 

11580 

11590 

11591 

11592 
11600 
11601 
11610 
11620 
11621 
11622 
11630 
11640 

11650 

11651 
11660 
11670 
11680 
11690 
11700 
11710 
11720 
11730 
11740 
11750 

11760 

11761 

11770 

11771 
11780 
117  81 
11782 
11790 
11800 
11810 
11820 
11830 
11840 
11850 
11860 
11870 
11880 
11890 
11900 
11910 
11920 
11930 
11940 
11950 


NP  (292 ,1)  =3 
SIP  (293,1)  =3 
NP(294,1) =1 
NP (289 ,1) =1 
NP (290 ,1) =3 
NP(291,1)=4 
NP  (291 ,2)  =1 
NP (291 ,3) =1 
NP  (2  86 ,1)  =1 
NP (286 ,2) =4 
NP  (2  87  ,1)  =3 
NP (288 ,1) =3 
NP (288 , 2) =2 
NP(288,3) =2 
NP  ( 2  83  , 1")  =1 
NP (284 ,1) =4 
NP(285,1)=2 
NP  (285 ,2)  =4 
NP  (2  80 ,1)  =4 
NP (281 ,1) =3 
NP (282 ,1) =2 
NP  (277  ,1)  =1 
NP ( 278 ,1) =1 
NP (279,1) =2 
NP (274 ,1) -2 
NP ( 275 ,1) =1 
NP (276 ,1) -A 
NP (271 ,1) -3 
NP (272 ,1) -2 
NP (272,2) =4 
NP ( 273 , 1) =2 
NP (273 ,2) =4 
NP (268 ,1) =3 
NP (268,2) =1 
NP  (268 ,3)  =1 
NP (269,1) =4 
NP (270 ,1) =4 
NP ( 265 ,1 ) =1 
NP (266 ,1) =3 
NP  (267 ,1)  =3 
NP (262 ,1) -2 
NP  (263 ,1)  =4 
NP ( 264 ,1) =4 
NP ( 259 ,1) =2 
NP  ( 260 ,1)  =1 
NP ( 261 , 1) =2 
NP ( 256 ,1) =2 
NP  (257  ,1)  =4 
NP (258 , 1) =4 
NP  (253  ,1)  =1 
NP (254 , 1) =3 
NP  (255 ,1)  =2 
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2/2 

NL 


11960 

11970 

11980 

11990 

12000 

12110 

12020 

12030 

12040 

12050 

12051 
12060 
12070 
12080 
12090 
12100 
12110 
12111 
12112 
12120 
12130 
12140 

12150 

12151 
12160 
12170 
12180 
12190 
12200 
12210 
12211 
12220 
12230 
12240 
12250 
12252 
12260 
12261 
12270 
12280 

12290 

12291 

12292 
12300 
12310 
12320 
12330 
12340 
123  50 
12352 

12360 

12361 


NP (250 ,1) =1 
NP  < 251 ,1) =4 
NP  (252  /l)  =3 
NP (247 ,1) =3 
NP (248  f 1) =2 
NP (249 ,1) =2 
NP(244 ,1) =2 
NP (245 ,1) =4 
NP (246 /l) =4 
NP (241 ,1)  =4 
NP (241 , 2} =2 
NP  (242 ,1)  =4 
NP (243  f 1) =4 
NP  ( 23  8 ,1)  =2 
NP (239  ,1) =1 
NP  (240 ,1)  =2 
NP (235 ,1) =3 
NP (235  #2) =2 
NP(235 ,3) =2 
NP (236  ,1)  =1 
NP  ( 237  ,1)  =3 
NP (232 , 1) =1 
NP (233 ,1) -2 
NP  (233 ,2)  =3 
NP(234,1)=2 
NP ( 229, 1) =3 
NP(230 ,1) =3 
NP  (231 ,1)  =2 
NP (226 , 1) *1 
NP (227,1) =1 
NP (227 , 2) =3 
NP (228,1) =1 
NP(223,1)=1 
NP (224 ,1) -2 
NP(225 ,1) =4 
NP (225,3) =2 
NP  (220  ,1)  =1 
NP  (220 ,2)  -2 
NP (221 ,1) =4 
NP (222,1) =2 
NP  ( 217 ,1)  =3 
NP  (217 ,2)  =1 
NP  ( 217  ,3 )  =1 
NP  ( 21 8,1)  =2 
NP ( 219 , 1) =2 
NP (214,1) =2 
NP  ( 215 , 1)  =1 
NP  (216 ,1)  =3 
NP ( 211 ,1) =1 
NP ( 211 ,3 ) =4 
NP  (  212 , 1)  =3 
NP(212,2) =2 
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12362 

12370 

12371 

12372 
123  80 
123  90 
12400 
12410 
12420 
12430 
12440 
12450 
12460 
12470 
12480 
12490 
12500 
12510 
12520 

12530 

12531 

12532 
12540 
12550 
12560 
12570 
125  80 
12590 
12600 
12610 
12620 
12630 
12640 

12650 

12651 

12652 
12660 
12670 
12680 
12681 

12690 

12691 

12692 
12700 
12710 
12720 
12730 
12740 
12750 
12760 

12770 

12771 


NP(212,3)=4 
NP(213,l)-2 
NP (213 ,2) =4 
NP (213  >3 )  =4 
NP (208 ,1) =3 
NP (209,1) =4 
NP  (210 ,1)  =3 
NP (205 ,1) =2 
NP(206,1)=1 
NP  (207 ,1)  =3 
NP (202 ,1) =3 
NP  (203 ,1)  =4 
NP(204 ,1) -2 
NP (199,1) =3 
NP(200 ,1) =4 
NP ( 201 ,1) =4 
NP(196 ,1)  =3 
NP (197 ,1 ) =1 
NP (198 ,1) “3 
NP (193 ,1) =2 
NP  (193  ,2)  =4 
NP (193 ,3) =4 
NP (194 ,1) =4 
NP (195 ,1) -4 
NP (190 ,1) =3 
NP (191,1) =2 
NP  (192 ,1)  =*2 
NP  ( 1 87 , 1 )  =3 
NP (188 , 1) =3 
NP (189,1) =4 
NP (184 ,1) =1 
NP (185 ,1) *1 
NP (186 ,1) =1 
NP (181,1) =3 
NP  (181 , 2)  =1 
NP (181 ,3 ) =1 
NP  (182,1)  =3 
NP (183 ,1) =2 
NP (178 ,1) =4 
NP(178,2) -2 
NP(179,1)=3 
NP (179,2) =4 
NP  (179 ,3)  *4 
NP (180 ,1) *4 
NP (175  ,1) =2 
NP (176 , 1) =1 
NP ( 177 , 1) =3 
NP (172,1) =1 
NP(173 ,1) -2 
NP (174,1) =4 
NP (169 ,1) »2 
NP  (16 9, 2)  =4 


87 


12772 
12780 
127  90 
12800 
12810 
12820 
12  830 
12840 
12850 
12860 
12870 
12880 
12890 
12900 
12910 
12920 
12930 
12940 
12950 
12960 
12970 
12980 
12990 
13000 
13010 
13020 
13030 
13040 
13050 
13060 
13070 
13080 
13090 
13100 
13110 
13120 
13130 
13140 
13150 
13160 
13170 
13180 
13190 
13200 
13202 
13210 
13220 
13230 
13240 
13250 
13260 
13270 


NP (169 
NP ( 170 
NP (171 
NP (166 
NP (167 
NP (168 
NP ( 163 
NP  (164 
NP ( 165 
NP (160 
NP ( 161 
NP (162 
NP ( 1 57 
NP ( 1 5  8 
NP (159 
NP ( 1 5  4 
NP (155 
NP (156 
NP (151 
NP (152 
NP (153 
NP (148 
NP (149 
NP (150 
NP (145 
NP (146 
NP ( 1 47 
NP (142 
NP  (143 
NP (144 
NP (139 
NP (140 
NP ( 1 41 
NP  (136 
NP (137 
NP (13  8 
NP (133 
NP (13  4 
NP(135 
NP (130 
NP(131 
NP (132 
NP(127 
NP  (128 
NP(128 
NP (129 
NP ( 124 
NP ( 1 2  5 
NP(126 
NP (121 
NP ( 1 2  2 
NP (123 


/3)  =4 
,1)=4 
,1)=2 
r  1)  *2 
,1)=1 
,1)=3 
,1)=3 
,1)=4 
,1)=3 
i-l)  =1 

f  1)  =2 

rl)=l 

,1)=3 

,1)=4 

f  1)  =2 
rl)=2 

/ 1)  =4 
f  1 )  =4 
i-l)  =2 
,1)=1 
,  1)  “3 
rl)  =3 
f  1)  -4 
f  1)  =3 
,1)=2 
/1)=1 
,1)=2 
rl)=l 
f  1)  =4 
#1)  -2 
,1)=1 
f  1)  =3 
a)  =3 
a)  =2 

a)=4 

a)=2 
a)  =3 
a)  =4 
a)=4 

rl)  =3 

a)=i 
a)=i 
a)=2 
a)  =2 

» 3)  =4 

a)  =4 
a)  =3 

rl)=l 

a)  =3 
a)  =2 
a)=i 

a)  =3 


13280 
13290 
13300 
13310 
13320 
13330 
13340 
13350 
13360 
13370 
133  80 
133  90 
13  400 

13401 

13402 
13410 
13  420 
13  430 

13440 

13441 

13442 

13450 

13451 

13452 
13  460 
13  470 
13480 
13490 
13500 
13510 

13520 

13521 

13522 

13530 

13531 
13540 
13550 
13  560 
13562 
13570 

13580 

13581 

13582 
13590 
13600 
13610 
13620 
13630 
13640 
13650 
13660 
13670 


NP(118,1)=3 
NP(119,1) *2 
NP(120,1)=2 
NP(115,1)=3 
NP(116,1)*4 
NP(117,1)=4 
NP (112 f 1) =2 
NP  (113  /l)  “1 
NP(114,1)=2 
NP (109fl) =2 
NP(110,1)=1 
NP  (111  ,1)  =3 
NP (106 ,1) =3 
NP  (106  f  2)  =1 
NP (106  f 3) =1 
NP  (107  r  1)  =1 
NP (108 /l) -2 
NP  (103 ,1)  =1 
NP (104 ,1) =2 
NP (104/2) =1 
NP(104,3)=1 
NP (105 , 1) =1 
NP (105 ,2) =4 
NP (105  f3 ) =4 
NP (100  rl) -2 
NP(101/1)=1 
NP (102 /l) =2 
NP  ( 97  / 1 )  =2 
NP  (98  /l)  =1 
NP  ( 99  /  l)  =3 
NP (94/1) =3 
NP(94,2)=1 
NP  (94  f3)  -1 
NP  (95  /  l)  =3 
NP  (95/2)  -2 
NP  (96  /  l)  =2 
NP  ( 91  /l)  =1 
NP  (92  rl)  =4 
NP  (92  f  3)  =2 
NP  ( 93  f  1 )  =3 
NP(88/1) *3 
NP  ( 8  8  /  2 )  =1 
NP(88,3)=1 
NP  ( 89  /l)  =4 
NP  (90  /l)  *3 
NP  ( 85  f  1)  *1 
NP ( 86 ,1) -1 
NP  (87  rl)  =1 
NP  (  82  / 1)  *3 
NP  ( 83  rl)  “2 
NP  ( 84  /  I)  *2 
NP (79/1) *4 
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m 

L-'- 

13680 

NP(80,1)=4 

13690 

NP  ( 81 ,1)  =4 

13700 

NP(76,1)=4 

( 

13710 

NP  (77 ,1)  =4 

C-V 

13720 

NP(78,1)=2 

^V 

■.*  • 

13730 

NP (73 ,1)  =3 

13740 

NP(74,1) =1 

13750 

NP (75 ,1) =1 

yJL 

13760 

NP(70,1)  =3 

13761 

NP  (70 ,2)  =1 

13762 

NP (70 ,3) =1 

.  v 

13770 

NP (71,1) =3 

•V 

l'  ■ 

13780 

NP(72,1)=4 

13790 

NP  (67  ,1)  =4 

h2 

13791 

NP(67,2)=1 

13792 

NP  (67  ,3)  =1 

13  800 

NP(68 ,1) =3 

13  810 

NP  (69/1)  =3 

•V 

13  820 

NP  (64  / 1)  =4 

13  821 

NP (64/2) =2 

-»4 

13830 

NP (65 ,1) =3 

13  831 

NP  ( 65 ,2)  =4 

13832 

NP  (65 ,3)  =4 

13840 

NP  (66  ,1)  =3 

13  842 

NP(66/3)=2 

13  850 

NP(61/1) =2 

m 

bn 

13860 

NP (62 /l) =1 

13  87  0 

NP  ( 63 , 1)  *2 

:y 

13880 

NP (58/1) =2 

13  890 

NP  (59,1)  =4 

13900 

NP (60  ,1)  =3 

■r. 

13  910 

NP  (55 ,1)  =1 

s 

13920 

NP  ( 56 ,1)  =3 

13  930 

NP  (57 ,1)  =2 

«• 

13940 

NP(52 ,1) =3 

13950 

NP  ( 53 , 1 )  =3 

13960 

NP (54 ,1) =2 

>  ' 

13  970 

NP  (49,1)  *1 

■—A 

13  972 

NP ( 4  9 , 3 ) =2 

13980 

NP  (50 ,1)  =4 

13982 

NP  (50 ,3)  -2 

f',\ 

13  990 

NP (51 ,1)  =1 

>  .* 

14000 

NP  (46,1)  -2 

14010 

NP  ( 47  ,1)  =4 

4 

14020 

NP(48 ,1) =3 

14030 

NP (43  ,1)  ®4 

**• 

14031 

NP  ( 43 , 2)  =2 

>v 

14040 

NP (44,1) =4 

14050 

NP  ( 45 ,1)  =4 

14060 

NP  ( 40 ,1 )  =2 

l‘d  14070 

NP  (41,1)  *1 

14080 

NP(42 ,1) =2 

£ 
.  V 

90 

i  .* 

£ 

*"-**•'* 

14090 

14100 

14110 

14120 

14121 
14130 
14140 
14150 
14160 
14170 
14180 
14190 
14200 
14210 
14220 
14230 
14240 
14250 
14260 
14270 

142  80 
14290 
14300 
14310 
14320 

14330 

14331 

14332 
14340 
14350 
14360 
14370 

143  80 
14390 
143  91 
14392 
14400 
14410 

14420 

14421 
14430 
14440 
14450 
14460 
14470 
14480 
14490 
14500 
14510 
14520 
14530 
23  440 


NP  (37  ,1) =2 
NP  ( 3  8 , 1 )  =1 
NP (39 ,1) =2 
NP (34,1) =1 
NP (34 ,2) =3 
NP (35 ,1) =3 
NP (36 ,1) =2 
NP (31,1) =1 
NP(32 ,1) =2 
NP  (33 ,1)  =2 
NP (28 ,1) =2 
NP ( 2  9 , 1 ) =1 
NP (30 ,1) =2 
NP ( 2  5 , 1 ) =1 
NP  (26 ,1)  =3 
NP  (27 ,1)  =2 
NP  (22  ,1)  =1 
NP  (23 ,1)  =4 
NP(24 ,1) =4 
NP ( 1 9 , 1 ) =3 
NP  (20  ,1)  =3 
NP(21 ,1) =2 
NP  ( 16 ,1)  =2 
NP  ( 17 ,1)  =1 
NP (18 ,1) =3 
NP  (13 ,1)  =3 
NP  (13  ,2)  =1 
NP  (13  ,3)  =1 
NP  ( 1 4 , 1 )  *3 
NP  ( 1 5 , 1 )  =2 
NP (10 ,1) =3 
NP  (11 ,1)  =4 
NP(12 ,1) =3 
NP  (7 ,1)  =2 
NP(7 ,2) =4 
NP (7 ,3 ) =4 
NP (8 ,1) =4 
NP ( 9 , 1 ) =4 
NP  (4 ,1)  =1 
NP  ( 4 , 2 )  =3 
NP ( 5 , 1 ) =  4 
NP  (6,1)  =1 
NP(1 ,1) =3 
NP ( 2 , 1 ) =4 
NP  (3 ,1)  =3 
DO  20  1=1,438 
IF (NP (1,2)  .EQ.  0)  NP ( I , 2 
IF (NP (1,3)  .EQ.  0)  NP  (1 ,3 
20  CONTINUE 
RETURN 
END 

SUBROUTINE  GC(MGC) 


) =NP (1,1) 
) =NP (1,1) 
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23  450 

DIMENSION  MGC (16,4) 

23  460 

MGC(1,1)=15 

23470 

MGC (2 ,1) =17 

23480 

MGC (3  ,1)  =4 

23490 

MGC (4 ,1) =16 

23500 

MGC (5,1) =3 

23510 

MGC (6 ,1) =20 

23520 

MGC (7,1) =4 

23530 

MGC (8,1) =25 

23540 

MGC ( 9,1) =15 

23550 

MGC (10 ,1) =17 

23560 

MGC (11 ,1) =4 

23570 

MGC (12,1) =16 

23580 

MGC (13,1) =3 

23590 

MGC (14 ,1) =25 

23600 

MGC (15 ,1) =4 

23610 

MGC (16,1) =25 

23620 

MGC (1,2) =5 

23630 

MGC  (2, 2)  =1 

23640 

MGC (3, 2) =2 

23650 

MGC (4, 2) =10 

23660 

MGC (5, 2) =5 

23670 

MGC (6 ,2) =1 

236  80 

MGC ( 7 , 2 ) =2 

23690 

MGC ( 8 , 2 ) =  8 

23700 

MGC (9,2) =5 

23710 

MGC  (10, 2)  =1 

23720 

MGC (11, 2) =2 

23730 

MGC(12 ,2) =10 

23740 

MGC (13 ,2) =5 

23750 

MGC (14, 2) =1 

23760 

MGC (15 ,2) =2 

23770 

MGC  (16, 2)  =8 

23780 

MGC  ( 1 , 3 )  =3 

23790 

MGC ( 2 , 3 ) =1 2 

23  800 

MGC (3 ,3 ) =9 

23810 

MGC ( 4 , 3 ) =1 8 

23  820 

MGC (5, 3) =3 

23830 

MGC ( 6 , 3 ) =1 2 

23  840 

MGC (7, 3) =9 

23850 

MGC ( 8 , 3 ) =1 4 

23  860 

MGC (9,3) =3 

23  870 

MGC (10 ,3) =12 

23  880 

MGC (11, 3) =9 

23890 

MGC (12 ,3) =18 

23900 

MGC (13 ,3) =3 

23910 

MGC (14,3) =12 

23  920 

MGC (15 ,3) =9 

23  930 

MGC (16 ,3) *14 

23  940 

MGC (I, 4) =11 

23950 

MGC (2 ,4) =4 

23  960 

MGC (3,4) =7 
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23  970 
23  980 
23990 
24000 
24010 
24020 
24230 
24040 
24050 
24060 
24070 
24080 
24090 
24100 
24110 
40010 
40020 
40025 
40030 
40040 
40050 
40060 
40070 
400  80 
40090 
40100 
40105 
40110 
40120 
50010 
50020 
50030 
50040 
50050 
50060 
50070 
500  80 
500  90 
50100 
50110 
50120 
60010 
60020 
60030 
60040 
60050 
60060 
60070 
60080 
60090 
60100 
60110 


MGC (4 ,4) =13 
MGC(5,4)=19 
MGC (6 ,4) =12 
MGC  (7  ,4)  =2 
MGC (8, 4) =6 
MGC (9, 4) =11 
MGC (10 ,4) =4 
MGC (11, 4) =7 
MGC (12  ,4) =13 
MGC  (13 , 4)  =11 
MGC  (14  /  4)  =12 
MGC (15, 4) =7 
MGC (16 ,4) =6 
RETURN 
END 

SUBROUTINE  INIT1 (NP , MP ,MGC) 

DIMENSION  NP ( 450 , 4) ,MP(150,4) , MGC (16, 4) 
DO  10  1=1,3 
J=1 46 
12  J3=J*3 
J31=J3-1 
J32=J3-2 

NIP=4*NP(J32,2) +NP (J31,2) -4 
MP ( J, 2) =MGC (NIP ,NP ( J3 , 2) ) 

J=J-1 

IF  (J.GE.l)  GO  TO  12 
10  CONTINUE 
RETURN 
END 

subroutine  init3 (np,mp,mgc) 

DIMENSION  NP ( 450 ,4) ,MP(150,4) ,MGC(16,4) 
J=1 46 
12  J3=J*3 
J31=J3-1 
J32=J3-2 

NIP=4*NP ( J3  2 ,3 ) +NP ( J3 1 ,3 ) -4 
MP ( J,3) =MGC (NIP,NP ( J3 , 3) ) 

J=J-1 

IF  (J.GE.l)  GO  TO  12 

RETURN 

END 

SUBROUTINE  AAAR(MCRIT) 

DIMENSION  MCRIT(146 ,14) 

MCRIT (1,3)=  7 
MCRIT(1 , 4)  =19 
MCRIT (1 ,5) =  9 
MCRIT (1 ,6) =  1 
MCRIT (1,7) =25 
MCRIT (2 ,2) =10 
MCRIT ( 2 ,3 ) =12 
MCRIT (2, 4) =18 
MCRIT (2,5)=  3 


93 


60120 

MCRIT{2,6)  =14 

60130 

MCRIT(2,7)=  8 

60140 

MCRIT (2 , 8) =25 

60150 

MCRIT(2,9) =13 

60160 

MCRIT (2 ,10) =16 

60170 

MCRIT (3 ,2)  =  3 

60180 

MCRIT(3,3)=  6 

60190 

MCRIT  (3,4)  =15 

60200 

MCRIT (3 ,5) =20 

60210 

MCRIT (3 ,6) =25 

60220 

MCRIT(4,2)=  7 

60230 

MCRIT  (4,3)=  4 

60240 

MCRIT  (4,4)=  2 

60250 

MCRIT (4, 5) =25 

60260 

MCRIT(5 ,2) =  9 

60270 

MCRIT (5, 3) =10 

60280 

MCRIT (9, 5) =13 

60290 

MCRIT (12,11)=  5 

60300 

MCRIT  (13,2)=  2 

60310 

MCRIT (12,7)=  3 

60320 

MCRIT (12 , 8) =10 

60330 

MCRIT (12, 9) =12 

60340 

MCRIT (12 ,10) =11 

60350 

MCRIT (13 ,3) =  1 

60360 

MCRIT  (13,4)=  7 

60370 

MCRIT  (13,5)=  4 

603  80 

MCRIT  (13, 6)  =17 

60390 

MCRIT (13, 7) =13 

60400 

MCRIT (14 ,2) =  3 

60410 

MCRIT (9, 6) =  7 

60420 

MCRIT (9 ,7) =  5 

60430 

MCRIT (9, 8) =10 

60440 

MCRIT (9, 9) =17 

60450 

MCRIT (9,10)=  1 

60  460 

MCRIT (9 ,11) =14 

60470 

MCRIT (10,2)=  2 

60480 

MCRIT (10 ,3)  =  7 

60490 

MCRIT (10 ,4) =18 

60500 

MCRIT  (10,5)=  3 

60510 

MCRIT (10 ,6 ) =1 4 

60520 

MCRIT (10 ,7) =  4 

6053C 

MCRIT (10, 8) =  8 

60540 

MCRIT (11 ,2) =  5 

60550 

MCRIT (11, 3) =  1 

60560 

MCRIT (11, 4) =11 

60570 

MCRIT (11,5) =  3 

60580 

MCRIT(11 ,6) =14 

60590 

MCRIT (11, 7) =  4 

60600 

MCRIT (11,8)-  8 

60610 

MCRIT (12,2) =  6 

60620 

MCRIT (12 ,3) =  7 

60630 

MCRIT (12,4)=  4 
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60640 
60650 
60660 
60670 
60680 
60690 
60700 
60710 
60720 
60730 
60740 
60750 
60760 
60770 
607  80 
60790 
60800 
60810 
60820 
60830 
60840 
60850 
60860 
60870 
60880 
60890 
60900 
60910 
60920 
60930 
60940 
60950 
60960 
60970 
60980 
60990 
61000 
61010 
61020 
61030 
61040 
61050 
61060 
61070 
61080 
61090 
61100 
61110 
61120 
61130 
61140 
61150 


MCRIT(12,5) =13 
MCRIT(12,6)=  2 
MCRIT (5,4)  =  4 
MCRIT(5,5)=  1 
MCRIT (5 ,  6)  =  8 
MCRIT (5,7) =  2 
MCRIT (5 ,8) =  6 
MCRIT (5, 9) =25 
MCRIT (6 ,2) =  4 
MCRIT (6,3)=  8 
MCRIT (6 ,4) =  2 
MCRIT (6,5) =10 
MCRIT (6 ,6) =  1 
MCRIT (6, 7) =  6 
MCRIT  (6, 8)  =14 
MCRIT (6,9)=  9 
MCRIT (6 ,10) =  5 
MCRIT (6, 11) =25 
MCRIT (7 ,2) =  8 
MCRIT (7,3)=  2 
MCRIT  (7,4)=  1 
MCRIT (7,5) =10 
MCRIT (7, 6) =21 
MCRIT (7, 7) =13 
MCRIT  (7, 8)  =14 
MCRIT (7,9)=  6 
MCRIT  (8,2)=  6 
MCRIT (8, 3) =10 
MCRIT ( 8 , 4) =  4 
MCRIT ( 9,2) =18 
MCRIT (9 ,3) =  4 
MCRIT (9,4)=  2 
MCRIT (14 ,3) =  4 
MCRIT (14, 4)  =15 
MCRIT (14 ,5)  =12 
MCRIT (15, 2) =20 
MCRIT (15 ,3) =15 
MCRIT (15, 4) =  1 
MCRIT(15 ,5) =  4 
MCRIT (16 ,2) =  1 
MCRIT (16 ,3) =  4 
MCRIT (16 ,4) =  2 
MCRIT (16, 5) =10 
MCRIT (16 ,6) =12 
MCRIT  (17,2)=  6 
MCRIT  (17, 3)  =11 
MCRIT (17,4)=  8 
MCRIT (17 ,5) =18 
MCRIT (18 ,2) =  5 
MCRIT (18, 3) =  6 
MCRIT (18, 4) =11 
MCRIT (18,5) =18 


61160 
61170 
61180 
61190 
61200 
61210 
61220 
61230 
61240 
61250 
61260 
61270 
61280 
61290 
61300 
61310 
61320 
61330 
61340 
61350 
61360 
61370 
61380 
613  90 
61400 
61410 
61420 
61430 
61440 
61450 
61460 
61470 
61480 
61490 
61500 
61510 
61520 
61530 
61540 
61550 
61560 
61570 
615  80 
61590 
61600 
61610 
61620 
61630 
61640 
61650 
61660 
61670 


MCRIT(19,2) =13 
MCRIT (19,3)=  6 
MCRIT(19,4) =18 
MCRIT ( 1 9 , 5 ) =21 
MCRIT(19,6)=  1 
MCRIT(19,7)=14 
MCRIT  (19 ,8)=  8 
MCRIT (19 ,9) =  2 
MCRIT(19 ,10) =10 
MCRIT (20/2) =  5 
MCRIT (20  ,3)  =  2 
MCRIT  (20 ,  4)  =  4 
MCRIT (20 ,5)=  8 
MCRIT (20,6)=  3 
MCRIT(20,7)=11 
MCRIT (20 / 8) =10 
MCRIT  (20,9)=  1 
MCRIT (20 ,10) =  9 
MCRIT (21, 2) =10 
MCRIT (21,3)=  8 
MCRIT (21, 4) =13 
MCRIT ( 21 , 5 ) =1 8 
MCRIT (21,6)=  2 
MCRIT (21 ,7) =  1 
MCRIT  (21, 8)  =14 
MCRIT (21, 9) =13 
MCRIT (21,10)=  8 
MCRIT (22, 2) =  8 
MCRIT (22, 3) =10 
MCRIT (22, 4) =  6 
MCRIT (22,5)=  2 
MCRIT (22, 6) =13 
MCRIT (22 ,7) =14 
MCRIT (23 ,2) =  5 
MCRIT (23, 3) =10 
MCRIT (23 ,4) =  2 
MCRIT  (23,5)=  7 
MCRIT (23, 6) =17 
MCRIT ( 23 , 7 ) =1 1 
MCRIT (23 ,8) =25 
MCRIT  (23,9)=  3 
MCRIT (24,2)=  1 
MCRIT (24 ,3) =  3 
MCRIT ( 2  4 , 4 ) =1 2 
MCRIT (24 ,5) =  5 
MCRIT (25, 2) =  1 
MCRIT(25 ,3) =  6 
MCRIT (25 ,4) =  2 
MCRIT (25 ,5) =12 
MCRIT (26 ,2) =  8 
MCRIT (26 ,3) =14 
MCRIT (26 ,4) =16 
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62200 
62210 
62220 
62230 
62240 
62250 
62260 
62270 
62280 
62290 
62300 
62310 
62320 
62330 
623  40 
62350 
62360 
62370 
623  80 
623  90 
62400 
62410 
62420 
62430 
62440 
62450 
62460 
62470 
62480 
62490 
62500 
62510 
62520 
62530 
62540 
62550 
62560 
62570 
62580 
62590 
62600 
62610 
62620 
62630 
62640 
62650 
62660 
62670 
62680 
626  90 
62700 
62710 


MCRIT (49,2)=  4 
MCRIT (49,3)=  1 
MCRIT ( 50 , 2)  =  7 
MCRIT(50,3)=  4 
MCRIT  (50 ,4)  =13 
MCRIT (50 ,5) =10 
MCRIT  (50 ,6)=  6 
MCRIT (51,2)=  9 
MCRIT (51 ,  3}  =  2 
MCRIT ( 52  f 2) =10 
MCRIT(52,3) =13 
MCRIT (52,4)=  4 
MCRIT  (52, 5)  =18 
MCRIT (52 ,6) =  1 
MCRIT (52 ,7) =  6 
MCRIT (52,8)=  8 
MCRIT (52,9)=  2 
MCRIT  (53,2)=  2 
MCRIT  (54,2)=  5 
MCRIT (54, 3) =11 
MCRIT  (55, 2)  =19 
MCRIT (55 ,3 ) =  3 
MCRIT (55, 4) =17 
MCRIT (56 ,2) =  1 
MCRIT (56, 3) =13 
MCRIT  (56,4)=  4 
MCRIT (56, 5) =10 
MCRIT (56 ,6) =  2 
MCRIT ( 56 , 7 ) =1 8 
MCRIT (56, 8) =25 
MCRIT (57, 2) =13 
MCRIT (57, 3) =10 
MCRIT  (57,4)=  2 
MCRIT (58,2)=  9 
MCRIT (58 ,3) =  2 
MCRIT (58,4) =12 
MCRIT (58 ,5) =25 
MCRIT (59, 2) =  6 
MCRIT (59 ,3) =  8 
MCRIT (59,4) =  7 
MCRIT  (59, 5)  =2  5 
MCRIT (59,6) =19 
MCRIT (59 ,7) =14 
MCRIT (60,2)=  5 
MCRIT (61,2)=  6 
MCRIT (6 1,3) =  3 
MCRIT (61 , 4) =13 
MCRIT (61,5) =  8 
MCRIT(61 ,6) =12 
MCRIT (62,2)=  2 
MCRIT (62 ,3) =  1 
MCRIT (63 ,2) =18 


62720 

MCRIT(63 ,3) =12 

62730 

MCRIT (63,4) =16 

62740 

MCRIT(64,2)=  1 

62750 

MCRIT (65/2) =  6 

62760 

MCRIT(65,3)=  4 

62770 

MCRIT (65/4) =  2 

62780 

MCRIT  (65/5)=  8 

627  90 

MCRIT (66/2)=  6 

62800 

MCRIT (66 / 3 ) =  8 

62810 

MCRIT (67 /2) =  5 

62820 

MCRIT (67 /3) =  2 

62830 

MCRIT (67 /4) =10 

62840 

MCRIT (67,5)=  8 

62850 

MCRIT (68 ,2) =  3 

62860 

MCRIT (68 ,3) =11 

62870 

MCRIT (69,2) =  8 

62880 

MCRIT (69 ,3) =12 

62890 

MCRIT (6 9, 4) =  9 

62900 

MCRIT (69 ,5) =  1 

62910 

MCRIT (69,6) =  7 

62920 

MCRIT (69 ,7) =  2 

62930 

MCRIT (69,8) =  4 

62940 

MCRIT (69, 9) =10 

62950 

MCRIT (69, 10) =13 

62960 

MCRIT ( 6  9 , 1 1 ) =1 8 

62970 

MCRIT ( 69 ,12) =  5 

62980 

MCRIT (69 ,13) =14 

62990 

MCRIT  (70,2)=  2 

63000 

MCRIT  (70,3)=  7 

63010 

MCRIT (70 ,4)  =  4 

63020 

MCRIT (71 ,2) =15 

63030 

MCRIT (71,3)=  3 

63040 

MCRIT ( 7 1 , 4 ) =1 1 

63050 

MCRIT (71 ,5) =  4 

63060 

MCRIT  (72,2)=  4 

63070 

MCRIT (7 2, 3) =17 

63080 

MCRIT (72 , 4) =  1 

63090 

MCRIT (72 ,5) =  2 

63100 

MCRIT (72 ,6) =10 

63110 

MCRIT (72 ,7) =13 

63120 

MCRIT (73 ,2) =10 

63130 

MCRIT (73 ,3 ) =  8 

63140 

MCRIT (73 ,4) =14 

63150 

MCRIT (73, 5) =13 

63160 

MCRIT (73 ,6) =  4 

63170 

MCRIT (74 ,2)  =  1 

63180 

MCRIT (74 ,3)  =  2 

63190 

MCRIT (74 ,4) =10 

63200 

MCRIT (75 ,2) =  3 

63210 

MCRIT (75 ,3 ) =  5 

63220 

MCRIT  (75 , 4)  =19 

63230 

MCRIT (75 ,5 ) =11 

63240 

MCRIT (75  ,  6)  =  7 

63250 

MCRIT (76  ,2)  =  2 

63260 

MCRIT(76,3)=13 

63270 

MCRIT (76 ,4) =  4 

63280 

MCRIT(76 ,5)  =  1 

63290 

MCRIT (76  ,6)  =18 

63300 

MCRIT (76 ,  7)  =  6 

63310 

MCRIT (76, 8) =  1 

63320 

MCRIT  (76, 9)  =14 

63330 

MCRIT (76 ,10)  =  8 

633  40 

MCRIT (77 ,2)  =18 

63350 

MCRIT (77  ,3 )  =13 

63360 

MCRIT  (77,4)=  8 

63370 

MCRIT ( 77 , 5 ) =1 4 

63380 

MCRIT (77, 6) =10 

63390 

MCRIT ( 77 , 7 ) =1 4 

63400 

MCRIT  ( 78 , 2)  =  3 

63410 

MCRIT (7 8, 3) =  9 

63420 

MCRIT (79 ,2) =10 

63430 

MCRIT (7 9, 3) =  8 

63440 

MCRIT (79, 4) =13 

63450 

MCRIT (80, 2) =13 

63460 

MCRIT ( 80 ,3) =10 

63470 

MCRIT  (80,4)=  6 

63480 

MCRIT  (80,5)=  4 

63490 

MCRIT (81,2)=  3 

63500 

MCRIT  (81, 3)  =11 

63510 

MCRIT (82,2)=  6 

63520 

MCRIT  (83,2)=  2 

63530 

MCRIT  (83,3)=  1 

63540 

MCRIT  (83, 4)  =13 

63550 

MCRIT (84,2)=  7 

63560 

MCRIT ( 84,3)=  2 

63570 

MCRIT (84, 4) =18 

635  80 

MCRIT (85, 2) =15 

63590 

MCRIT ( 85 ,3 ) =16 

63600 

MCRIT  (86,2)=  2 

63610 

MCRIT  (86,3)=  4 

63620 

MCRIT  (87,2)=  7 

63630 

MCRIT  (87,3)=  4 

63640 

MCRIT  (87,4)=  6 

63650 

MCRIT (  87 ,5)  =14 

63660 

MCRIT  (87,6)=  2 

63670 

MCRIT  (87,7)=  8 

63680 

MCRIT  (87, 8)  =18 

63690 

MCRIT (87, 9) =13 

63700 

MCRIT (88,2)=  3 

63710 

MCRIT (88, 3) =12 

63720 

MCRIT (88,4)=  9 

63730 

MCRIT (89,2)=  4 

63740 

MCRIT (89 ,3) =  7 

63750 

MCRIT (90,2)=  8 

63760 

MCRIT (90/3) =14 

63770 

MCRIT (  90 , 4)  =  6 

63780 

MCRIT  (91 ,2)=  3 

63790 

MCRIT (91 ,3) =16 

63  800 

MCRIT  (91,4)=  9 

63810 

MCRIT (91 ,5) =25 

63820 

MCRIT  (92, 2)  =18 

63830 

MCRIT (92, 3) =16 

63840 

MCRIT (92, 4) =25 

63850 

MCRIT (93, 2)  =17 

63860 

MCRIT  (93,3)=  4 

63  870 

MCRIT ( 93 , 4 ) =2  5 

63880 

MCRIT (94, 2) =10 

63890 

MCRIT (94,3)=  5 

63  900 

MCRIT (94, 4) =13 

63910 

MCRIT (95,2)=  6 

63920 

MCRIT (95,3)=  2 

63930 

MCRIT  (95, 4)  =25 

63940 

MCRIT (95 ,5) =  8 

63950 

MCRIT ( 95 , 6 ) =1 4 

63  960 

MCRIT  (96,2)=  3 

63970 

MCRIT (97, 2) =18 

63980 

MCRIT (97, 3) =14 

63990 

MCRIT (97, 4) =12 

64000 

MCRIT (98,2)=  5 

64010 

MCRIT (98, 3) =19 

64020 

MCRIT (99, 2) =10 

64030 

MCRIT (99, 3) =13 

64040 

MCRIT (99, 4) =18 

64050 

MCRIT (99,5) =16 

64060 

MCRIT (100 ,2) =  9 

64070 

MCRIT (101 ,2) =  8 

640  80 

MCRIT (101, 3) =14 

640  90 

MCRIT(101,4) =  2 

64100 

MCRIT (102 ,2)  =13 

64110 

MCRIT (102,3)  =10 

64120 

MCRIT(102,4)  =  6 

64130 

MCRIT (102,5)=  7 

64140 

MCRIT(103 ,2) =15 

64150 

MCRIT (104,2) =12 

64160 

MCRIT (104 ,3)  =  6 

64170 

MCRIT (10 4 ,4) =13 

641 80 

MCRIT(i05 ,2) =  3 

64190 

MCRIT (105,3) =12 

64200 

MCRIT(106 ,2) =  3 

64210 

MCRIT(107 ,2)  =  1 

64220 

MCRIT(108 ,2) =13 

64230 

MCRIT (108,3)  =10 

64240 

MCRIT (10 9 ,2)  =  5 

64250 

MCRIT (10 9, 3) =19 

64260 

MCRIT  (109 , 4)  =  4 

64270 

MCRIT (10 9, 5) =11 
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64280 

MCRIT(109,6) -  2 

64290 

MCRIT (110  f2) =  3 

64300 

MCRIT(110,3) =11 

64310 

MCRIT { 11 0 , 4 ) =1 5 

64320 

MCRIT (110 ,5)=  5 

64330 

MCRIT (111 ,2) =  5 

64340 

MCRIT(111,3)=21 

64350 

MCRIT ( 111 ,4) =  2 

64360 

MCRIT(111,5) =11 

64370 

MCRIT (112 ,2) =17 

643  80 

MCRIT(112  ,3) =11 

64390 

MCRIT (112,4)  =  4 

64400 

MCRIT{112,5)=  5 

64410 

MCRIT (112, 6) =18 

64420 

MCRIT(112 ,7)  =11 

64430 

MCRIT (113 ,2) =  5 

64440 

MCRIT (113 ,3) =17 

64450 

MCRIT (113, 4) =  8 

64460 

MCRIT(113 ,5) =  2 

64470 

MCRIT (11 4, 2) =  3 

64480 

MCRIT (115 , 2) =  2 

64490 

MCRIT (115 ,3 ) =  1 

64500 

MCRIT(115 , 4) =  4 

64510 

MCRIT (116, 2) =18 

64520 

MCRIT (116, 3) =10 

64530 

MCRIT (116, 4) =12 

64540 

MCRIT(116 ,5)  =  7 

64550 

MCRIT (116 ,6) =  8 

64560 

MCRIT (116 ,7) =13 

64570 

MCRIT (117, 2) =18 

64580 

MCRIT (117, 3) =13 

64590 

MCRIT (117, 4) =12 

64600 

MCRIT (11 8 ,2) =15 

64610 

MCRIT (11 8, 3) =  3 

64620 

MCRIT (118, 4) =18 

64630 

MCRIT (11 9, 2) =  1 

64640 

MCRIT (11 9 ,3) =14 

64650 

MCRIT(119,4) =15 

64660 

MCRIT  (119,5)=  4 

64670 

MCRIT (119,6)=  6 

6  4680 

MCRIT(119 ,7) =10 

6  4690 

MCRIT (120 ,2) =  6 

64700 

MCRIT (120, 3) =10 

64710 

MCRIT(120 ,4) =13 

64720 

MCRIT(120 ,5) =  4 

64730 

MCRIT (120, 6) =18 

64740 

MCRIT (120 ,7) =  8 

64750 

MCRIT (121,2) =  8 

64760 

MCRIT (121, 3) =10 

64770 

MCRIT ( 1 21 , 4 ) =1 4 

647  80 

MCRIT (121 ,5) =13 

64790 

MCRIT (121,6)=  6 
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64800 

64810 

64820 

64830 

64840 

64850 

64860 

64870 

64880 

64890 

64900 

64910 

64920 

64930 

64940 

64950 

64960 

64970 

64980 

64990 

65000 

65010 

65020 

65030 

65040 

65050 

65060 

65070 

650  80 

65090 

65100 

65110 

65120 

65130 

65140 

65150 

65160 

65170 

65180 

65190 

65200 

65210 

65220 

65230 

65240 

65250 

65260 

65270 

65280 

65290 

65300 

65310 


MCRIT (122, 2) =15 
MCRIT (123 ,  2 )  =  7 
MCRIT ( 1 23 , 3 )  =1 0 
MCRIT (123 ,4) *  4 
MCRIT(123 ,5) =13 
MCRIT(124,2) =  9 
MCRIT (124 ,3) =11 
MCRIT (124/4) =12 
MCRIT(125,2)=  9 
MCRIT (125, 3) =  2 
MCRIT (125 ,4) =14 
MCRIT (125 ,5) =12 
MCRIT (125 ,6)  =  8 
MCRIT(125 ,7)  =  5 
MCRIT(125 , 8)  =  3 
MCRIT(125 ,9)  =  1 
MCRIT (125 ,10)  =10 
MCRIT (125 ,11) =17 
MCRIT (126,2)=  5 
MCRIT (126 ,3 ) =  2 
MCRIT (126 , 4) =  3 
MCRIT (126 ,5)  =19 
MCRIT(126 ,6)  =  7 
MCRIT (126 ,7) =  8 
MCRIT (126 ,8) =17 
MCRIT (127,2) =14 
MCRIT (127 ,3)  =  8 
MCRIT  (127 ,4)  =  3 
MCRIT (12 8 ,2)  =  2 
MCRIT (12 8 ,3)  =  5 
MCRIT{128 ,4) =18 
MCRIT(128,5) =  4 
MCRIT (12 9 ,2) =  2 
MCRIT (129,3)=  4 
MCRIT(129 , 4) =10 
MCRIT (129,5) =14 
MCRIT (129 ,6) =  8 
MCRIT (12 9,7) =  1 
MCRIT(129,8) =13 
MCRIT (130 ,2) =16 
MCRIT (130 ,3) =20 
MCRIT (130,4) =18 
MCRIT(130,5) =  3 
MCRIT (130,6) =15 
MCRIT (130 ,7) =10 
MCRIT (131,2) =14 
MCRIT(131,3) =  8 
MCRIT (131 , 4 ) =  3 
MCRIT(132,2) =  6 
MCRIT (13 2,3 ) =  2 
MCRIT  (13 2 ,4)  =14 
MCRIT (13 2 ,5) =  8 
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65320 
65330 
65340 
65350 
65360 
65370 
653  80 
653  90 
65400 
65410 
65420 
65430 
65440 
65450 
65460 
65470 
65480 
65490 
65500 
65510 
65520 
65530 
65540 
65550 
65560 
65570 
65580 
65590 
65600 
65610 
65620 
65630 
65640 
65650 
65660 
65670 
65680 
65690 
65700 
65710 
65720 
65730 
65740 
65750 
65760 
70080 

70090 

70091 
70100 

70110 

70111 
70120 


MCRIT(133 ,2)  =  5 
MCRIT ( 133  f 3 ) =1 9 
MCRIT(133,4)  =  3 
MCRIT (133, 5) =18 
MCRIT (13 4 ,2) =  5 
MCRIT (134,3) =15 
MCRIT (13 5 ,2) =  2 
MCRIT (13 5 ,3) =  7 
MCRIT (13 5, 4) =17 
MCRIT (13 5, 5) =  4 
MCRIT (13 5, 6) =12 
MCRIT (13 6 ,2) =  1 
MCRIT (136 ,3) =  2 
MCRIT (136,4) =10 
MCRIT(136 ,5) =  5 
MCRIT (13 7 ,2) =  5 
MCRIT (13 8 ,2) =  1 
MCRIT (13 8,3 ) =  2 
MCRIT (13 8 ,4) =  4 
MCRIT (13 9, 2)  =13 
MCRIT (13 9 ,3) =  2 
MCRIT (13 9, 4) =10 
MCRIT(139 ,5) =  7 
MCRIT (13 9, 6) =  4  . 

MCRIT (139 ,7) =18 
MCRIT (140 ,2) =  2 
MCRIT(141 ,2) =  3 
MCRIT ( 1 41 , 3 ) =1 2 
MCRIT (142 ,2) =  2 
MCRIT (142 ,3) =  4 
MCRIT (142 ,4) =  1 
MCRIT (143 ,2) =18 
MCRIT (143 ,3) =  4 
MCRIT (143 ,4) =  6 
MCRIT (143 ,5)  =12 
MCRIT (144,2) =  6 
MCRIT(144 ,3) =12 
MCRIT (144, 4) =  2 
MCRIT (145,2) =16 
MCRIT ( 145 ,3 ) =17 
MCRIT  (145 ,4)  =18 
MCRIT (146, 2) =18 
MCRIT (146 ,3) =10 
RETURN 
END 

SUBROUTINE  INIT(IX,IXS,F,E,N, INTER , NCOUNT , 
1MCOUNT , NULL , NA , MA , M , MGC , NP , NI P , NR , MR , MP , IH , 

1MCH , MCR , MCA) 

INTEGER  E,F 

DIMENSION  NA ( 43  8 ,50 ) ,MA(146,50)  ,NULL(438)  , 

1PCT ( 43  8) , NCOUNT (4,4) 

DIMENSION  MCOUNT(28,28) ,MGC(16,4) ,N(450) ,M(150) , 


104 


70121  1NP(450,4) 

70130  DIMENSION  NR(43 8) ,MR (146) ,MP (150 ,4) ,MCH (150) , 

70131  1MCR(150) , MCA (150) 

70140  IX=IXS 

70150  E=1 

70160  F=1 

70170  MCH(147) =0 

70180  MCR ( 1 47 ) =0 

70190  MCA (147) =0 

70200  N(439)  =2 

70210  DO  13  1=1 ,4 

70220  DO  13  J=l,4 

70230  NCOUNT(I,J) =0 

70240  13  CONTINUE 

70250  DO  14  J=l,28 

70260  DO  14  1=1,28 

70270  MCOUNT ( I ,  J)  =0 

70280  14  CONTINUE 

70290  DO  15  1=1,438 

70300  NULL ( I ) =0 

70310  NA ( I , 1 ) =1 

70320  NR (I) =0 

70330  PCT(I) =0.000 

70340  DO  151  J=2 ,50 

70350  NA  (I ,  J)  =0 

70360  151  CONTINUE 

70370  15  CONTINUE 

70380  DO  141  1=1,146 

70390  MA(I  ,1)  =1 

70400  MR (I) =0 

70410  M ( I ) =MP ( I , IH) 

70420  MCH ( I ) =0 

70430  MCR (I) =0 

70440  MCA ( I ) =0 

70450  DO  1411  J=2 , 50 

70460  MA(I , J) =0 

70470  1411  CONTINUE 

70480  141  CONTINUE 

70490  DO  101  1=1,438 

70500  N ( I) =NP ( I , IH) 

70510  101  CONTINUE 

70520  DO  1011  1=1,146,7 

70530  1011  CONTINUE 

70540  RETURN 

70550  END 

80010  SUBROUTINE  TRANF ( A, B1 ,B2 ,B3 ,C1 ,C2 ,C3 ,D1 ,D2 ,D3 , El 

80011  1E1,E2,E3,K) 

80020  INTEGER  K 

80030  A-0.3971 

80040  IF  (K.EQ.2)  GO  TO  20 

80050  IF  (K.EQ.3)  GO  TO  30 

80060  C  *******FROM  ZUCKERKANDL(1971) *********** 


80070  C1=0 .061 

800  80  C2=0  .0  97 

80090  C3=0 .166 

80100  Dl=0.085 

80110  D2=0.183 

80120  D3=0 .342 

80130  El=0.038 

80140  E2=0.130 

80150  E3=0 .219 

80160  B1=0 .058 

80170  B2=0.179 

80180  B3=0 .273 

80190  GO  TO  40 

80200  20  CONTINUE 

80210  C  *******FROM  FITCH (1967) *********** 

80220  C1=0 .0113 

80230  C2=0 .0113 

80240  C3=0 .0715 

80250  01=0.0458 

80260  D2=0 .1447 

80270  D3=0 .2407 

80280  E1=0 .0277 

80290  E2=0 .0  912 

80300  E3=0 .2907 

80310  B1  =0.1140 

80320  B2=0 .2377 

80330  B3  =0 . 3  97 1 

80340  GO  TO  40 

80350  30  CONTINUE 

80360  C  *****************FROM  DAYHOFF-P.226*********** 

80370  C1=0 .0244 

80380  C2=0  .0  975 

80390  C3=0 .1463 

80400  D1=0 .0122 

80410  D2=0  .0 976 

80420  D3=0 .3  903 

80430  El =0.0610 

80440  E2=0 .1586 

80450  E3=0 .2318 

80460  B1=0 .0366 

80470  B2=0 .1464 

80480  B3=0 .2318 

80490  40  RETURN 

80500  END 

90010  SUBROUTINE  TRAN1 (NUC , I , INTER, N, IX) 

90020  DIMENSION  N(450) 

90030  NUC=N ( I ) 

90040  B=0. 066666667 

90050  C=0. 033333333 

90060  CALL  RANDOM  (IX,P) 

90070  IF  (P.GT.0.1)  GO  TO  33 

90080  IF  (P.LE.0.0)  GO  TO  33 


s.  *.  .  V.\ 
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900  90 

IF  (P.GE.B)  GO  TO  31 

90100 

IF  (P.LE.C)  GO  TO  32 

90110 

N(I)=N(I)+1 

90120 

GO  TO  33 

90130 

31 

N(I)=N(I)+2 

90140 

GO  TO  33 

90150 

32 

N(I)=N(I)+3 

90160 

33 

CONTINUE 

90170 

IF  (N (I)  .GE.5)  N(I)=N(I)-4 

90180 

RETURN 

90190 

END 

100010 

SUBROUTINE  MUTATE ( I , MCRIT , NULL , NUC , NCOUNT , 

100020 

1MC0UNT,NA,MA,N,M,MGC,E,F,NR,MR) 

100030 

INTEGER  E,F 

100040 

DIMENSION  NA ( 43  8 , 50)  ,MA(146,50)  ,NULL(43  8)  , 

100041 

INCOUNT (4/4) 

100050 

DIMENSION  MCOUNT(28,28)  , MCRIT (146 ,1 4) ,N(450 

10C060 

DIMENSION  MGC (16,4) ,NR(438) ,MR(146) ,M(150) 

100070 

JIN=l+((I-l)/3) 

100080 

J3=3*JIN 

100090 

J31=J3-1 

100100 

J32=J3-2 

100110 

NIP=4*N(J32) +N( J31) -4 

100120 

IF  (M(JIN) . EQ.MGC (NIP,N ( J3) ) )  GO  TO  30 

100130 

50 

ICRIT=JIN 

100140 

K=1 

100150 

DO  51  10=2,14 

100160 

IF  (MGC(NIP,N(J3) ) .EQ. MCRIT (ICRIT,IC) )  K=2 

100170 

51 

CONTINUE 

100180 

IF  (K.EQ.2)  GO  TO  34 

100190 

NULL ( I ) =NULL ( I ) +1 

100200 

N ( I ) =NUC 

100210 

GO  TO  30 

100220 

34 

CONTINUE 

100230 

NCOUNT ( NUC ,  N  ( I ) } =NCOUNT ( NUC ,  N  ( I )  )  +1 

100240 

MCOUNT (M ( JIN) ,MGC (NIP,N( J3) ) ) =MCOUNT (M ( JIN) 

100241 

1MGC(NIP,N(J3) )  )+l 

100250 

E=NA(I,1) +1 

100260 

IF  (E.GE.50)  E=50 

100270 

NA(IfE) =NUC 

100280 

IP1=0 

100290 

DO  130  IP=2 ,50 

100300 

IPIT=NA(I , IP) 

100310 

IF  (IPIT.EQ.N(I) )  IP1=1 

100320 

130 

CONTINUE 

100330 

NA (1,1) =NA ( I , 1) +1 

100340 

F=MA( JIN , 1) +1 

100350 

IF  (F.GE.50)  F=50 

100360 

1 

MA( JIN,F) =M (JIN) 

100370 

IP2=0 

100380 

DO  131  IP=2 ,50 

1003  90 

IPIT=MA ( JIN , IP) 

100400 

100410 

131 

IF  (IPIT.EQ.M ( JIN) )  IP2=1 

CONTINUE 

'  * 

100420 

MA(JIN,1)=MA(JIN,1)  +1 

( 

100430 

M ( JIN) =MGC (NIP, N ( J3 ) ) 

100440 

30 

CONTINUE 

100450 

IF  (IP1.EQ.1)  NR ( I ) =NR ( I ) +1 

v’; 

V.* 

100460 

IF  (IP2.EQ.1)  MR (JIN) =MR ( JIN) +1 

100470 

IP1=0 

» v 

100480 

IP2=0 

M 

100490 

RETURN 

p 

100500 

END 

k  *  . 

110010 

SUBROUTINE  OUT ( M , N , NA , MA , NULL , MCOUNT , NCOUNT , 

110020 

1 HB , IX  S , METH , INTER , NR , MR , MCH , MCR , MCA , MP ) 

110030 

DIMENSION  NA( 438,50) ,MA(146,50) ,NULL(438) , 

110031 

1PCT  ( 43  8) 

* 

110040 

DIMENSION  NCOUNT (4,4) , MCOUNT (2 8, 2 8)  ,N(450)  , 

Si 

110041 

1M (150) 

110050 

DIMENSION  NR ( 43  8)  ,MR(146)  ,MCH(150)  ,MCR(150)  , 

E-S 

110051 

1MCA(150) 

Sj] 

110060 

DIMENSION  MP (150 ,4) 

110070 

DO  1000  1=1,21 

a 

110080 

MCOUNT (I, 21) =MCOUNT(I ,25) 

110090 

MCOUNT (21,1) =MCOUNT (25,1) 

110100 

MCOUNT (21, 21) =MCOUNT(25 ,25) 

110110 

DO  1001  J=1 ,21 

110120 

IVAL=MCOUNT ( I , J) 

’•  ■' 

110130 

WRITE  (10,100)  (HB,METH, IX S, INTER, IVAL) 

110140 

1001 

CONTINUE 

*V 

110150 

1000 

CONTINUE 

%? 

110160 

DO  1002  1=1,4 

O 

110170 

DO  1003  J=1 ,4 

'_V 

110180 

IVAL=NCOUNT ( I , J) 

110190 

WRITE  (11,100)  (HB,METH,IXS, INTER, IVAL) 

1 

110200 

1003 

CONTINUE 

•V 

110210 

1002 

CONTINUE 

yv 

110220 

DO  1005  1=1,438 

j- 

110230 

NA  (1,1) =NA ( 1 , 1 ) -1 

110240 

IVAL=NA  (1,1) 

*  *  * 

110250 

WRITE  (13,100)  (HB,METH,IXS, INTER, IVAL) 

m 

110260 

1005 

CONTINUE 

k.  '*  • 

110270 

DO  1006  1=1,438 

110280 

IVAL=NR (I) 

,y, 

110290 

WRITE  (14,100)  (HB, METH, IXS, INTER, IVAL) 

1 

110300 

1006 

CONTINUE 

110310 

DO  1007  1=1,146 

110320 

IVAL=MR ( I ) 

110330 

WRITE  (15,100)  (HB, METH, IXS, INTER, IVAL) 

% 
i.»  < 

110340 

1007 

CONTINUE 

i ,  * 
\\ 

110350 

DO  1004  1=1,438 

!  *J* 

110360 

IDEM1=I 

'/m  * 
-  1 

110370 

IVAL=NULL ( I ) 

P 

1103  80 

WRITE  (12,100)  (HB, METH, IXS, INTER, IVAL) 
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r . 

i 

” 

f  1103  90 

1004 

CONTINUE 

110400 

DO  1  J-1,3 

i  110410 

DO  2  1=1,146 

110420 

IP  (MP(I,J) .EQ.M(I) )  GO  TO  3 

110430 

GO  TO  2 

110440 

3 

IP  (J.EQ.l)  MCH ( I ) =1 

110450 

IP  (J.EQ.2)  MCR (I) =1 

110460 

IF  (J.EQ.3)  MCA(I) =1 

i  110470 

2 

CONTINUE 

1  110480 

1 

CONTINUE 

S  t  110490 

MCH (147) =0 

110500 

MCR (147 ) =0 

\  110510 

DO  4  1=1,146 

3  -  110520 

MCH  ( 1 47 )  =MCH  (147)  +MCH  ( I ) 

i  110530 

MCR (147) =MCR (147) +MCR ( I ) 

110540 

MCA(147)=MCA(147)+MCA(I) 

110550 

4 

CONTINUE 

110560 

DO  5  1=1,147 

110570 

WRITE  (20,100)  (HB ,METH, IXS , INTER, MCH ( I) ) 

;  110580 

WRITE  (21,100)  (HB,METH,IXS, INTER, MCR(I) ) 

i  110590 

WRITE  (22,100)  (HB, METH, IXS, INTER, MCA(I) ) 

110600 

5 

CONTINUE 

110610 

100 

FORMAT  (10X, 12, IX, 12, IX, 18, IX, 16, IX, 18) 

110620 

RETURN 

110630 

END 

120010 

SUBROUTINE  RANDOM ( IX, P) 

120020 

IY=IX*1220703125 

120030 

IF (IY) 3,4,4 

‘  120040 

3 

IY=IY+21 47483647+1 

120050 

4 

RN-IY 

120060 

RN«RN*0. 46 56613 E-9 

120070 

P»24.3*RN-10. 

120080 

IX=IY 

!  120090 

RETURN 

120100 

END 

130010 

SUBROUTINE  RESULT ( M , N , NA , MA , NULL , MCOUNT , NCOUNT) 

i  130020 

DIMENSION  NA(  438 , 50)  ,MA(146  ,50)  ,NULL(438)  , 

;  130021 

1PCT(43  8) 

130030 

DIMENSION  NCOUNT (4, 4) , MCOUNT (2 8, 2 8) ,N(450) , 

130031 

1M  ( 1 5  0 ) 

130040 

IXS-3513 

;  130050 

PRINT  38, (M( J) , J*1 ,146) 

130060 

PRINT  45, ( (MCOUNT (I, J) ,1=1,25) ,J=1,25) 

.  130070 

PRINT  391 

130080 

PRINT  39,  (N ( I )  ,1=1,438) 

130090 

PRINT  44, ( (NCOUNT (I, J) ,1=1,4) ,J=1,4) 

l  130100 

PRINT  421 

130110 

DO  40  1=1,43  8 

•  130120 

NA1=NA (1,1) 

:  130130 

IP  (NA1.EQ.1)  GO  TO  40 

|  130140 

IP  (NA1.GE.50)  NA1=50 

;  130150 

PRINT  42,1, (NA(I,J) ,J=1,NA1) 

i» 
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,•  V  V  \*  V  V  .•  /  V 

130160 
130170 
130180 
130190 
130200 
130210 
130220 
130230 
130240 
130250 
130260 
130270 
130280 
130290 
130300 
130310 
130320 
130330 
130340 
130350 
130351 
130360 
130370 
1303  80 
1303  90 
1303  91 
130400 
130410 
130420 
130430 
130440 
130441 
130450 
130451 
130460 
130470 
130471 
130480 
130481 
130490 
130491 
130500 
130501 
130510 
130511 
130520 
130521 
130530 
130531 
130540 
130550 
130560 


40  CONTINUE 
PRINT  431 

DO  41  1=1,146 
MAI  =MA  (1,1) 

IP  (MA1.EQ.1)  GO  TO  41 
IF  (MA1.GE.25)  MAI =2 5 
PRINT  43,1, (MA(I,J) ,J=1, MAI) 

41  CONTINUE 
PRINT  551 

DO  56  1=1,438,6 

IF  (NULL (I) . EQ.0)  GO  TO  56 

PRINT  55 ,I,NULL(I) ,NULL(I+1) ,NULL(I+2) , 

1NULL ( 1+3) ,NULL (1+4) ,NULL(I+5) 

56  CONTINUE 
PRINT  531 
DO  54  1=1,438 

PCT  ( I ) =1 0  0 . *NA ( 1 , 1 ) / ( NA (1,1) +NULL ( I )  ) 

54  CONTINUE 

DO  52  1=1,438,6 

PRINT  53,1 ,PCT (I ) ,PCT (1+1) ,PCT(I+2) ,PCT(I+3) , 
1PCT (1+4) , PCT (1+5) 

52  CONTINUE 

36  FORMAT  (15HOON  INTERATION  ,I7,10H  POSITION  ,14, 
16H  FROM  ,I2,4H  TO  ,12) 

37  FORMAT  (4H  AA  ,I3,14H  CHANGED  FROM  ,I2,4H  TO  , 
112) 

38  FORMAT  (2013) 

39  FORMAT  (3012) 

42  FORMAT  (14 ,2X,I2 ,4913) 

43  FORMAT  (14, 2X, 12, 2415) 

44  FORMAT  (26H0NUCLEID  ACID  CHANGE  TABLE/15X , 2HTO/ 
16H  FROM  ,3X,1HA, 

13X,1HC,3X,1HG,3X,1 HT/ 3X , 1 HA , 2X , 4 I 4/ 3X , 1 HC , 2X , 
14I4/3X  ,1HG,2X  ,414/ 

13X,1HT,2X,4I4) 

45  FORMAT  (25H0  AMINO  ACID  CHANGE  TABLE/3 9X , 2HTO/ 
16H  FROM  , 

IIHG ,  2X ,1HA, 2X , 1HL, 2X , 1HS , 2X , 1HV,2X ,1HK , 2X , 1HT, 
12X ,1HE , 2X ,1HP , 

1 2X , 1 HD , 2X , 1HI , 2X , 1HR , 2X , 1 HN , 2X , 1 HQ , 2X , 1 HF , 2X , 
1H7 , 2X , 1HC ,  2X , 

IIHH, 2X,1HM,2X,1HW,12X,4HST0P/2H  G,2X,25I3/2H  A, 
12X ,25I3/2H  L, 

12X,25I3/2H  S, 2X , 25I3/2H  V  2X,25I3/2H  K,2X,25I3/ 
12H  T,2X ,2513/ 

12H  E,2X,25I3/2H  P,2X,25I3/2H  D,2X,25I3/2H  I,2X, 
125I3/2H  R,2X , 

125I3/2H  N, 2X , 25I3/2H  Q,2X,25I3/2H  F,2X,25I3/ 

12H  Y,2X,25I3/2H  C, 

12X,25I3/2H  H, 2X , 25I3/2H  M,2X,25I3/2H  W,2X,25I3/ 
14X, 25I3/4X, 2513/ 4X,25I3/4X,25I3/1X,4HST0P, 2513) 

53  FORMAT (I4,6F10.3) 


130570 

130580 

130590 

130600 

130610 

130611 

130620 

130621 

130630 

130640 

130641 

130650 

130660 

140010 

140020 

140030 

140040 

140050 

140060 

140070 

140071 

140080 

140090 

140100 

140110 

140120 

140130 

140140 

140150 

140160 

140170 

140180 

140190 

140200 

140210 

140220 

140230 

140240 

140250 

140260 

140270 

140280 

140290 

140300 

140310 

140320 

140330 

140340 

150010 

150020 

150030 

150040 


55  FORMAT  (14,615) 

381  FORMAT (21H0STARTING  AMINO  ACIDS/) 

382  FORMAT (19H0 ENDING  AMINO  ACIDS/) 

391  FORMAT (21 HO ENDING  NUCLEIC  ACIDS/) 

421  FORMAT  (34H0  NUCLEOTIDE  MUTATIONS  BY  POSITION// 
15H  NA  ,5HQUANT/) 

431  FORMAT  (34H0  AMINO  ACID  MUTATIONS  BY  POSITION// 
15H  AA  , 5HQUANT/) 

531  FORMAT  ( 29H0PERCENT  ACCEPTABLE  MUTATIONS/) 

551  FORMAT  (46H0  COUNT  OF  DISALLOWED  MUTATIONS  BY  N 
1A  POSITION) 

RETURN 

END 

SUBROUTINE  TRAN 2 ( NUC , I , INTER , N , K , IX ) 

INTEGER  K 
DIMENSION  N ( 450 ) 

NUC=N (I) 

CALL  RANDOM  (IX,P) 

IF  (P.LT.0.0)  GO  TO  341 

CALL  TRANF ( A , B1 , B2 , B3 , Cl ,C2 ,C3 ,D1 ,D2 ,D3 rEl ,E2 , 
1E3 ,  K) 

IF  (P.GT.A)  GO  TO  341 
NUC=N ( I ) 

IF  (NUC.EQ.l)  GO  TO  31 


IF 

(NUC.EQ.2) 

GO 

TO  3  2 

IF 

(NUC.EQ.3) 

GO 

TO  33 

IF 

(P.GE.B3) 

GO 

TO  341 

IF 

(P.LE.B1) 

N  (I 

)  =N  ( I )  -1 

IF 

(P.GE.B2) 

N  ( I 

) =N ( I ) +1 

N  ( ] 

[ ) =N( I ) -2 

GO 

TO  341 

31 

IF 

(P.GE.C3) 

GO 

TO  3  41 

N  () 

C)  =N  (I)  +2 

IF 

(P.LE.C1) 

N(I) 

=N(I)-1 

IF 

(P.GE.C2) 

N(I) 

=N (I) +1 

GO 

TO  3  41 

32 

IF 

(P.GE.D3) 

GO 

TO  3  41 

N  ( ] 

:)=N(I)+1 

IF 

(P.LE.D1) 

N(I) 

=N ( I ) -2 

IF 

(P.GE.D2) 

N  ( I ) 

=N ( I ) +1 

GO 

TO  3  41 

33 

IF 

(P.GE.E3) 

GO 

TO  3  41 

N  ( ] 

:)  =N(I)  -1 

IF 

(P.LE.E1) 

N  ( I ) 

=N(I)-1 

IF 

(P.GE.E2) 

N  ( I ) 

=N (I) +2 

341  CONTINUE 
RETURN 
END 

SUBROUTINE  TP(T) 
DIMENSION  T( 64 , 3 ) 
T(1 ,1)  =  .0540 
T(lr2) =.0016 


150570 

150580 

150590 

150600 

150610 

150620 

150630 

150640 

150650 

150660 

150670 

150680 

150690 

150700 

150710 

150720 

150730 

150740 

150750 

150760 

150770 

150780 

150790 

150800 

150810 

150820 

150830 

150840 

150850 

150860 

150  870 

150880 

150890 

150900 

150910 

150920 

150930 

150  940 

150950 

150960 

150  970 

150980 

150990 

151000 

151010 

151020 

151030 

151040 

151050 

151060 

151070 

151080 


T  (19 

.1)  = 

■  .0530 

T  (19 

,2)  = 

.0198 

T  (19 

,3)  = 

.0308 

T(20 

f  1 )  = 

.1260 

T(20 

f  2)  = 

.0342 

T(20 

,3)  = 

.0657 

T(21 

,1)  = 

.0240 

T(21 

f  2)  = 

.0072 

T  ( 21 

f  3)  = 

.0240 

T(22 

.1)  = 

.0650 

T(22 

,2)  = 

.0063 

T(22 

/  3 )  = 

.0503 

T  ( 23 

,1)  = 

.0870 

T(23 

/  2 )  = 

.0930 

T  ( 23 

r  3)  = 

.0217 

T(24 

.1)  = 

.0870 

T  (24 

f  2)  = 

.0032 

T(24 

/  3 )  = 

.0292 

T  (25 

/1)  = 

.1600 

T  ( 25 

/2)  = 

.0295 

T(25 

/  3 )  = 

.1558 

T(26 

r  1 )  ® 

.1330 

T  ( 26 

f  2)  = 

.0640 

T  ( 26 

r  3 )  = 

.0689 

T  ( 27  i 

rl)  = 

.0740 

T  (27  i 

r  2 )  = 

.0167 

T  ( 27  , 

r  3  )  = 

.0382 

T(28, 

r  1 )  ~ 

.1290 

T(28, 

r  2)  = 

.0165 

T  ( 2  8  ( 

r3)  = 

.0825 

T(29  ( 

r  1)  ~ 

.1600 

T  ( 29 1 

-  2 )  = 

.1333 

T  (29  i 

r3)  = 

.1333 

T(30, 

-1)  = 

.1350 

T (30  i 

'  2)  = 

.0053 

T(30 , 

3 )  = 

,0957 

T(31, 

1)  = 

.1240 

T(31, 

2)  = 

.0157 

T(31, 

3)  = 

.1033 

T(32, 

1)  = 

,0740 

T(32, 

2)  = 

,0072 

T(32  r 

3 )  = 

,0489 

T  (33  / 

1)  = 

,0540 

T(33f 

2)  = 

0135 

T(33  / 

3 )  = 

0270 

T(34, 

1)  = 

0990 

T<34, 

2)  = 

0198 

T(34, 

3 )  = 

0912 

T  (35  f 

1)  = 

0530 

T  ( 3  5  f 

2 )  = 

0217 

T(35  r 

3 )  = 

0337 

T  (36  f 

1)  = 

1280 
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151090 

T(36,2) 

=.0199 

151100 

T(36,3) 

=  .0995 

151110 

T(37,l) 

=  .0270 

151120 

T(37,2) 

=.0000 

151130 

T(37,3) 

=.0221 

151140 

T(38,l) 

=.0650 

151150 

T(38,2) 

=.0061 

151160 

1(38,3) 

=.0488 

151170 

T (39 ,1) 

=  .0870 

151180 

T(39,2) 

=  .0094 

151190 

T (39 ,3) 

=.0220 

151200 

T(40 ,1) 

=.1600 

151210 

T(40,2) 

=.0125 

151220 

T(40,3) 

=  .0972 

151230 

T(41 ,1) 

=.1600 

151240 

T( 41 ,2) 

=  .0205 

151250 

T  ( 41 , 3) 

=  .1486 

151260 

T(42,l) 

=.1330 

151270 

T  ( 42 ,2) 

=.0692 

151280 

T( 42 ,3 ) 

=.0745 

151290 

T  ( 43 , 1) 

=  .0740 

151300 

T(43 ,2) 

=  .0185 

151310 

T(43,3) 

=.0423 

151320 

T ( 44 , 1) 

=  .1290 

151330 

T(44 ,2) 

=  .0183 

1513  40 

T(44 , 3 ) 

=.1062 

151350 

T( 45 ,1) 

=.0550 

151360 

T(45,2) 

=.0467 

151370 

T (45  f 3) 

=.0517 

1513  80 

T(46 ,1) 

=  .1410 

1513  90 

T (46 ,2) 

=.0101 

151400 

T  (  46 ,3 ) 

=  .0831 

151410 

T(47,l) 

=  .0880 

151420 

T(47,2) 

=.0152 

151430 

T  ( 47 ,3) 

=.0684 

151440 

T  ( 48 ,1 ) 

=  .1780 

151450 

T (48 , 2) 

=.0123 

151460 

T(48,3) 

=.1666 

151470 

T(49 ,1) 

=  .0530 

151480 

T(49, 2) 

=.0035 

151490 

T(49 ,3) 

=  .0283 

151500 

T ( 50 ,1) 

=  .0990 

151510 

T (50 , 2) 

=  .0194 

151520 

T(50 ,3) 

=  .0893 

151530 

T  ( 51 ,1) 

=.0530 

151540 

T(51 ,2) 

=.0198 

151550 

T(51 ,3) 

=.0308 

151560 

T(52,l) 

=  .1280 

151570 

T ( 52 , 2) 

=.0259 

151580 

T(52 ,3) 

=  .0743 

151590 

T( 53 , 1) 

=  .0240 

151600 

T(53 ,2) 

=.0032 
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151610 

T(53  f 3)  =  .0240 

151620 

T<54,1) =.0650 

151630 

T(54  f  2)  = .0063 

151640 

T(54,3)=.0503 

151650 

T(55,l)  =  .0870 

151660 

T(55  ,2) =.0093 

151670 

T(55,3)  =  .0217 

1516  80 

T{56,1)=.0870 

151690 

T(56,2)=.0034 

151700 

T(56,3)=.0428 

151710 

T(57,l)=.1600 

151720 

T(57,2) =.0215 

151730 

T(57  ,3)  =  .1569 

1517  40 

T(58,l) =.1330 

151750 

T(58,2)=.0640 

151760 

T(58,3) =.0689 

151770 

T  (59  f  1)  =  .07 40 

1517  80 

T(59,2) =.0167 

151790 

T(59f3)=.0382 

151800 

T(60/l) =.1290 

151810 

T(60 ,2)  =  .0227 

151820 

T(60,3)=.0992 

151830 

T(61 f 1)  =  .1600 

151840 

T  ( 6 1 r  2 ) =.1333 

151850 

T(61  /3)  =  .1333 

151860 

T(62rl) =.1350 

151870 

T(62  f 2)  =  .0101 

151880 

T(62,3)=.0605 

151890 

T(63 fl) =.1240 

151900 

T(63,2) =.0231 

151910 

T(63f3)  =  .0  93  8 

151920 

T(64f 1) =.07  40 

151930 

T(64  f  2)  =  .0072 

151940 

T(64f3) =.0489 

151950 

RETURN 

151960 

END 

160080 

SUBROUTINE  TRAN3 (NUC f I , INTER, N, T, IX) 

160090 

DIMENSION  N(450) ,T(64,3) 

160100 

CALL  RANDOM  (IX fP) 

160110 

NUC=N (I) 

160120 

IP  (P.LT.0.0)  GO  TO  341 

160130 

IA=I+1 

160140 

IB=I-1 

160150 

IT=N(IA) +4* (N(I) -1) +16* (N(IB) -1) 

160160 

IF  (P.GE.T (IT, 1) )  GO  TO  341 

160170 

NUC=N { I ) 

160180 

IF  (NUC. EQ .1)  GO  TO  31 

160190 

IF  (NUC.EQ.2)  GO  TO  32 

160200 

IF  (NUC.EQ.3)  GO  TO  33 

160210 

IF  (P.LE.T(IT,2) )  N( I ) =N ( I ) -1 

160220 

IF  (P.GE.T(IT,3) )  N (I) =N (I) +1 

160230 

N( I ) =N( I ) -2 

160240 
160250 
160260 
160270 
160280 
160290 
160300 
160310 
160320 
160330 
160340 
160350 
160360 
160370 
1603  80 


GO  TO  341 

31  N ( I )  =N ( I ) +2 

IP  (P.LE.T ( IT, 2) )  N(I) =N(I) -1 
IF  (P.GE.T(IT,3) )  N(I)  =N(I) +1 
GO  TO  3  41 

32  N ( I ) =N ( I ) +1 

IP  (P.LE.T(IT,2))  N (I) =N(I) -2 
IP  (P.GE.T(IT,3) )  N ( I )  =N ( I ) +1 
GO  TO  3  41 

33  N ( I ) =N ( I ) -1 

IP  (P.LE.T(IT,2) )  N(I)-N(I)-1 
IF  (P .GE.T (IT, 3 ) )  N  ( I )  =N  ( I )  +2 
341  CONTINUE 
RETURN 
END 
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